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PREFACE 

The  Second  Plowshare  Symposium  was  held  in  Sap  Francisco,  California, 
on  May  13=15,  1959.  The  main  topics  covered  were:  Underground  Explosion 
Phenomenology,  Excavation  with  Nuclear  Explosives,  Power  and  Isotope  Re¬ 
covery  Studies,  Industrial  Applications,  and  Scientific  Experiments.  In  order 
to  publish  the  proceedings  in  volumes  of  a  reasonable  size,  each  of  the  major 
sessions  is  being  published  separately. 

As  is  indicated  by  the  topic  headings,  the  meeting  followed  a  planned  se= 
quence  from  the  past,  through  the  present,  and  ended  with  some  future  plans 
and  ideas.  The  order  of  presentation  was  (1)  the  data  on  past  underground  ex¬ 
plosions,  (2)  the  current  theories  on  phenomena  produced  in  underground  ex¬ 
plosions,  (3)  the  extrapolation  of  these  facts  and  theories  to  excavation  and 
power- recovery  projects  planned  for  the  near  future,  and  (4)  some  thoughts 
on  industrial  and  scientific  applications  envisioned  somewhat  farther  in  the 
future.  This  first  volume  contains  the  talks  of  the  first  day,  May  13,  con= 
cerning  the  underground  nuclear  explosions  that  took  place  in  Nevada  in  1957 
and  1958, 

In  addition  to  the  speakers,  many  people  contributed  to  the  successful 
carrying  out  of  the  symposium.  Marshall  Peterson  handled  the  meeting  ar¬ 
rangements;  Rodney  Southwick  was  responsible  for  the  news  coverage  and 
public  information  releases;  David  Rabb,  Nanette  Jaffe,  and  Harlan  Zodtner 
helped  in  editing  the  papers  and  session  recordings  for  publication;  the  staff 
of  the  Technical  Information  Division  was  responsible  for  the  preparation, 
publication  and  distribution  of  the  symposium  reports. 


W.  J.  Frank 
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A.  Opening  Talks 


1.  INTRODUCTORY  REMARKS 
Harold  Brown 

Lawrence  Radiation  Laboratory*  Livermore 

The  first  Plowshare  Symposium  was  held  somewhat  over  two  years  ago. 

It  would  seem  appropriate  at  this  time  to  take  a  look  at  where  we  now  stand, 
at  what  we*ve  done  since  that  first  Symposium.,  That  is  one  reason  for  this 
second  Symposium. 

Perhaps  an  even  more  important  reason  is  to  help  those  people  who  are 
interested  in  Plowshare,  or  have  heard  something  of  Plowshare,  learn  what*s 
going  on  in  the  subject,  by  hearing  the  people  who  are  working  in  it  say  what 
they  are  doing,  instead  of  by  reading  scattered  and  sometimes  overpopular¬ 
ized  articles.  In  particular,  we  hope  that  you  will  realize  that  this  program 
does  not  consist  of  a  group  of  physicists  in  a  back  room  drawing  conclusions 
about  geology,  chemical  processing,  or  economics,  by.t  that  qualified  people 
have  been  brought  in  as  participants  in  the  program  and  as  consultants  for 
each  specific  application  which  is  being  considered. 

I  don*t  think  it*s  my  job,  particularly  at  the  beginning  of  a  meeting  of  this 
kind,  to  evaluate  the  progress  of  the  last  two  years;  but  I  do  think  that  two 
points  should  be  made.  First,  there  have  now  been  a  number  of  underground 
nuclear  explosions.  These  have  not  been  Plowshare  experiments  in  the  sense 
that  they  were  entirely  motivated  by  the  Plowshare  program;  they  had  other 
motivations  as  well.  However,  they  have  cast  a  considerable  amount  of  light 
on  the  characteristic  effects  of  underground  explosions  and  on  the  specific 
Plowshare  applications  which  use,  or  might  use,  those  effects. 

Second,  in  evaluating  the  Plowshare  program,  one  important  criterion 
(it's  not  the  only  one)  is  whether  we  are  learning  more.  That  is,  a  criterion 
of  the  success  of  the  program  is  whether  we  are  learning  more  and  whether 
what  we  have  learned  is  favorable  or  unfavorable  to  the  specific  applications 
that  we  have  in  mind.  Although  I  think  the  accomplishment  of  specific  projects 
is  in  the  long  run  the  most  important  criterion  of  progress,  it  is  by  no  means 
the  only  test  of  progress,  and  I  think  we  should  keep  that  in  mind  as  we  listen 
to  the  subsequent  papers. 
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2.  WELCOMING  ADDRESS 
Ellison  Co  Shute 

San  Francisco  Operations  Office,  USAEC 

On  behalf  of  the  University  of  California  Lawrence  Radiation  Laboratory, 
and  the  Atomic  Energy  Commission,  I  extend  a  wholehearted  welcome  to  all 
of  you  attending  the  Second  Plowshare  Symposium,  It  is  our  hope  that  these 
discussions  will  stimulate  the  exchange  of  ideas  on  the  peaceful  applications 
of  nuclear  explosives  and  will  further  contribute  to  our  knowledge  in  this  field0 
The  intriguing  possibility  of  using  nuclear  explosives  for  peacetime  purposes 
was  considered  as  early  as  fourteen  years  ago  at  the  inception  of  the  atomic 
energy  program.  In  more  recent  years,  scientists  and  engineers  in  our  lab¬ 
oratories  and  in  industrial  and  private  institutions  have  continued  to  consider 
ways  and  means  in  which  this  new  explosive  could  be  put  to  work.  The 
Lawrence  Radiation  Laboratory  at  Livermore  provided  the  initiative,  starting 
a  formal  program  of  investigation  which  has  become  known  as  the  Plowshare 
program.  A  first  major  step  in  this  direction,  as  Dr.  Brown  indicated,  was 
the  First  Plowshare  Symposium  held  in  February  1957,  Although  it  was  con¬ 
ducted  as  a  classified  meeting,  most  of  the  papers  presented  there  were  sub¬ 
sequently  declassified  and  published.  Additional  emphasis  was  given  to  the 
program  by  the  underground  weapons  tests  conducted  at  the  Nevada  Test  Site 
in  1957  and  1958.  Those  tests  have  proved  to  be  a  valuable  source  of  refer¬ 
ence  data  in  considering  possible  uses  of  nuclear  detonations  for  entirely 
peaceful  purposes  such  as  mining,  excavation,  and  earth-moving  jobs  previ¬ 
ously  entirely  impractical  or  impossible  before  the  advent  of  this  new  dimen¬ 
sion  in  explosives.  Other  possibilities  worth  investigating  are  the  production 
of  power  and  isotopes,  and  the  replenishment  of  ground  water.  Nuclear  ex¬ 
plosives  have  been  known  to  the  public  only  as  nuclear  weapons  and  as  such 
have  gained  a  somewhat  fearsome  reputation.  History  indicates,  however, 
that  this  preoccupation  of  wartime  application  is  not  entirely  justified.  Many 
other  explosives  have  been  invented  for  wartime  use  and  were  considered 
horrible;  in  the  long  run,  however,  peacetime  uses  have  proved  to  be  far 
more  numerous  than  the  military  uses. 

Two  formal  projects  are  currently  under  way  in  the  Plowshare  program. 
They  are;  Project  Gnome  which  is  a  proposed  detonation  of  a  10-kiloton 
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device  in  the  salt  beds  near  Carlsbad,  New  Mexico,  for  the  purpose  of  inves¬ 
tigating  the  feasibility  of  recovering  power  and  isotopes,  and  Project  Chariot 
which  is  a  proposed  experiment  to  use  two  200-kiloton  and  three  20-kiloton 
devices  to  demonstrate  the  feasibility  of  excavation  for  such  purposes  as 
harbors  and  canals.  In  addition,  there  have  been  discussions  with  the  petro¬ 
leum  industry  regarding  its  interest  in  supporting  an  AEC-U,  S.  Bureau  of 
Mines  project  to  investigate  the  practicability  of  using  nuclear  devices  to 
crush  oil  shales  thereby  permitting  the  recovery  of  oil  by  in- situ  retorting. 
Let  me  re-emphasize  that  the  Gnome  and  Chariot  Projects  are  under  study. 
Thorough:  consideration  is  being  given  to  every  phase  of  these  projects  and 
the  Commission  requires  that  every  precaution  be  taken  to  assure  the  public 
health  and  safety.  Extensive  studies  and  surveys  are  being  and  will  be  con¬ 
ducted  to  determine  all  possible  effects  of  blast,  radiation,  and  heat  on  any 
environment  being  considered.  The  Commission  will  give  us  approval  to  pro¬ 
ceed  only  after  such  studies  have  established  that  the  nuclear  explosives  can 
be  used  with  complete  safety.  We  all  look  forward  to  hearing  the  ideas  and 
papers  being  presented  and  the  discussions  that  X  am  sure  will  follow.  I  hope 
you  will  find  the  sessions  both  informative  and  interesting. 
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3.  THE  PLOWSHARE  PROGRAM 
Edward  Teller 

Lawrence  Radiation  Laboratory,  Livermore 

It  is  a  very  great  pleasure  to  see  all  of  you  here.  It8s  rather  a  bigger 
gathering  than  we  had  two  years  ago,  and  I  think  that  what  we  know  about  the 
subject  is  also  a  little  bit  larger.  I  will  try  to  give  you,  in  a  very  few  minutes, 
some  sort  of  a  perspective  on  what  we  are  after. 

Conventional  explosives  have  been  used  in  peacetime  as  much,  perhaps 
more,  than  they  have  been  used  in  wartime.  They  have  been  used  in  the  main 
for  two  purposes.  They  have  been  used  for  transportation  in  the  sense  that 
obstacles  have  been  removed  in  building  roads,  waterways  and  other  similar 
things.  And  they  have  been  used  in  mining.  You  will  see  that,  according  to 
our  present  ideas,  these  are  also  the  two  main  uses  of  nuclear  explosives, 
although  you  will  hear  about  other  classes  of  possible  applications. 

The  reason  that  nuclear  explosives  have  been  so  slow  in  entering  our  ideas 
for  peaceful  purposes  is  twofold.  One  reason  is  that  nuclear  explosives  have 
sort  of  biggish  effects  with  the  result  that  one  tends  to  think  twice  or  perhaps 
even  more  than  twice  before  one  actually  tries  to  undertake  something.  The 
other  reason  is  that  nuclear  explosives  produce  radioactivity.  This  radio¬ 
activity  in  the  immediate  neighborhood  of  the  nuclear  explosive  can,  indeed, 
be  quite  dangerous;  one  has  to  take,  and  we  are  going  to  take,  as  Mr.  Shute 
has  just  told  you,  the  most  conscientious  precautions  that  this  radioactivity 
will  do  no  harm.  Of  course,  one  should  also  realize  that  radioactivity  in  ex¬ 
ceedingly  small  doses  (a  subject  that  often  enters  the  public  discussipn)  is 
really  far  from  dangerous;  one  cannot  even  determine  sometimes  whether  a 
danger  exists. 

I  should  say  right  away  that  many  of  the  applications  are  underground  ex¬ 
plosions  where  the  radioactivity  is  completely  contained  and  where,  provided 
one  pays  serious  attention  to  the  question  of  possible  ground  water  contamina¬ 
tion,  the  danger  will  be  quite  negligible.  Furthermore,  I  should  tell  you  that 
the  radioactivity  we  are  likely  to  release  underground  is  very  much  smaller 
than  the  radioactivity  that  one  buries  when  one  tries  to  dispose  of  reactor  pro¬ 
ducts.  The  difficulties  due  to  possible  contamination  merit  serious  and  con¬ 
scientious  consideration.  But  they  are  by  no  means  insurmountable. 
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It  has  been  an  exceedingly  wise  act  of  the  Atomic  Energy  Commission  to 
have  declassified  the  whole  subject  of  Plowshare  with  one  exception,  with  the 
exception  of  the  black  box  which  contains  the  main  actor  in  the  drama.  This 
item  we  cannot  at  the  present  time  declassify*  I  would  like  to  say  that  all  of 
us  who  are  involved  in  this  business  are  very  anxious  to  conduct  it  with  the 
greatest  possible  openness  and  to  give  you  information  about  the  nuclear  ex¬ 
plosive  itself  insofar  as  this  should  prove  compatible  with  the  information 
which  we  have  to  guard  for  military  reasons*  Actually,  almost  everything 
that  needs  to  be  said  for  an  effective  application  can  be  said  in  an  unclassified 
manner  and  will  be  said  here. 

I  have  told  you  that  two  great  difficulties  connected  with  Plowshare  are 
the  fact  that  nuclear  explosives  are  big  and  that  they  produce  radioactivity. 

I  would  like  to  assure  you  that  we  are  working  very  hard  on  these  two  topics. 

It  stands  to  reason  that  the  full  effect  of  nuclear  explosives  can  be  exploited 
only  if  the  biggest  kind  of  nuclear  explosions  are  included  in  our  program  and 
in  the  long  run  we  hope  to  include  them.  At  the  same  time  it  will  add  a  very 
great  deal  to  the  flexibility  and  the  scope  of  the  program  if  smaller  nuclear 
explosions  can  be  carried  out  and  particularly  if  they  can  be  carried  out 
cheaply.  This  is  one  of  the  aims  of  continued  research  into  nuclear  explosives 
and  I  only  wish  I  could  tell  you  more  about  this  exceedingly  interesting  topic. 
Perhaps  someday  I  shall  be  able,  or  somebody  shall  be  able,  to  do  so. 

Second,  I  should  like  to  tell  you  that  the  difficulty  of  radioactive  contami¬ 
nation  is  also  something  which  we  can  attack.  Even  at  the  present  time  we 
know  how  to  make  big  nuclear  explosions  fairly  clean.  We  are  hard  at  work 
to  make  big  and  small  nuclear  explosions  as  clean  as  possible  and,  hopefully, 
we  may  make  them  completely  clean.  It  will  require  invention,  it  will  require 
work,  but  if  we  could  get  rid  of  radioactive  contamination,  the  applicability  of 
this  new  tool  would  be  very,  very  greatly  increased,  and  this  is  one  of  the 
most  important  aims  toward  which  our  experiments  are,  and  must  be,  directed. 

Now,  as  to  the  program  itself,  we  need  help.  We  have  knowledge  —  good 
knowledge,  reliable  knowledge  —  about  nuclear  explosives:  how  to  make  them, 
how  to  handle  them,  how  to  apply  them..  We  acquired  some  knowledge  of  the 
behavior  of  nuclear  explosives  in  the  air  quite  some  time  ago,  and  underground 
more  recently.  This  knowledge  comes  in  very  handy,  but  we  have  to  have  more 
knowledge,  particularly  about  underground  nuclear  explosions. 
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However,  there  are  very,  very  big  areas  where  we  are  in  dire  need  of 
advice  and  cooperation;  in  fact,  this  program  cannot  prosper  without  a  great 
number  of  people,  earth  scientists,  chemists  of  every  description,  scientific 
people  of  every  kind.  If  these  people  will  participate,  and  I  know  many  of  you 
are  in  these  various  fields,  then  the  program  will  prosper;  otherwise  it  will 
not. 

Let  me  sketch  briefly  the  state  of  our  prospects  and  ideas.  First  of  all, 
we  have  the  earth-moving  projects  —  building  harbors,  building  canals.  To 
move  earth,  and  to  move  earth  cheaply,  is  something  that  we  know  we  can  do 
today e  We  do  not  know  the  best  way  to  do  it,  we  do  not  know  in  every  detail 
how  the  radioactivity  can  be  handled,  and  we  can  only  hope  that  the  radioac¬ 
tivity  might  be  avoided  altogether,  as  I  have  told  you.  But,  the  question  does 
not  end  here.  Today  we  can  dig  a  harbor  for  the  price,  or  for  less  than  the 
price,  that  it  costs  to  equip  that  harbor.  But,  when  you  want  to  dig  a  harbor, 
the  question  immediately  arises  "where?"  Where  is  it  worthwhile  to  dig  a 
harbor?  Will  that  harbor  last  or  will  it  silt  up?  Again,  for  a  canal,  the  same 
kind  of  economic  questions  arise.  These  are  tasks  which  we  can  perform, 
which  we  hope  to  perform  wisely  with  the  help  of  a  wide  and  thorough  discus¬ 
sion  both  with  respect  to  safety  and  with  respect  to  the  question  of  whether  the 
projects  are  worthwhile. 

Then,  there  are  many  projects  which  do  not  depend  on  the  moving  of  earth 
but  which  depend  on  the  breaking  up  of  rock.  These  include  possible  regula¬ 
tion  of  underground  flow  of  water,  storage  ,of  water,  and  all  manner  of  mining 
including  oil  production.  These  two  words,  water  and  oil,  are  magic.  Oil  is 
our  most  important  fuel  and  the  lifeblood  of  modern  industry.  It  may  be  a 
strange  reversal  when  nuclear  energy  will  be  used  to  make  the  main  source 
of  conventional  energy  cheaper  and  more  plentiful.  We  seem  to  be  putting  some 
cart  before  some  horse. 

Water  is  even  more  important;  but  whether  these  things  can  be  done  in  a 
practical  sense  is  an  exceedingly  difficult  question.  In  fact,  I  think  that  any 
thing  or  project  of  this  description  is  likely  to  turn  out  to  be  not  feasible,  at 
least  for  the  time  being.  We  have  found  that  we  can  crush  rock;  but,  we  also 
have  found  out  that  this  rock  may  solidify  again.  We  have  found  out  that  we 
can  break  up  rock,  but  only  one  kind  of  rock  —  volcanic  ash;  we  do  not  know 
whether  we  can  break  up  granite,  limestone,  or  oil  shale,  and  in  what  precise 
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manner.  Yet,  for  these  applications  of  water  storage,,  underground  regula¬ 
tion  of  water  flow,  mining,  oil  production,  these  questions  of  breaking  up  rock 
are  of  basic  significance. 

And  here  arises  the  point  that  Harold  Brown  has  so  correctly  mentioned  — 
we  must  for  the  next  few  years  consider  ourselves  successful  if  we  have 
gathered  information  about  these  important  things.  We  may  have  a  chance,  a 
small  chance,  a  10%  chance,  or  a  smaller  chance  even,  greatly  to  enrich  the 
oil  reserves  of  the  world  by  making  available  oil  sands  or  oil  shales.  The 
odds  may  be  long,  but  the  stakes  are  high  and  our  investment  need  not  be  very 
great. 

There  is  a  big  area  of,  in  my  opinion,  even  more  adventurous  thought. 

We  can  produce  heat  and  we  can  mine  the  heat.  If  we  can  make  the  capital 
investment  of  piping  and  other  things  serve  for  many  shots  in  the  same  place, 
we  can  hope  to  do  this  in  a  really  economical  way  and  to  have  really  cheap  and 
abundant  thermonuclear  power.  We  can  make  isotopes  of  all  descriptions 
underground;  some  of  these  isotopes  will  be  very  valuable,  provided  we  can 
extract  them  in  a  cheap  and  effective  manner. 

Now  I  would  like  to  mention  to  you  two  fields  which,  I  believe,  have  been 
considered  much  too  little.  They  have  been  considered,  but  the  consideration 
has  not  been  sufficient.  One  is  the  use  of  nuclear  explosions  for  purely  scien¬ 
tific  purposes.  A  nuclear  explosion  underground  is  a  source  of  geophysical 
disturbance,  a  smallish  earthquake  which  can  be  used  to  find  out  more  about 
the  crust  and  interior  of  our  earth.  A  nuclear  explosion  is  a  source  of  radia¬ 
tion  and  neutrons;  we  can  use  nuclear  explosives  to  find  out  more  about  neutrons 
and  nuclear  physics  in  general.  One  can  use  a  nuclear  explosive  as  it  was  used 
unintentionally  in  the  first  big  nuclear  explosion  in  the  Pacifier  to  capture 
many  successive  neutrons  in  the  same  nucleus  and  thereby  get  to  a  new  species 
which  then  can  be  studied.  One  can  use  nuclear  explosives  for -various  strange 
purposes  to  find  out  the  properties  of  the  upper  atmosphere  and  the  earth*s 
magnetic  field,  as  it  has  been  done  by  the  well-publicized  Argus  experiment. 

I  am  sure  that  this  does  not  begin  to  exhaust  the  list.  We  haven jt  thought  . 
enough  about  the  subject  and  we  shall  hear  about  it  during  this  symposium.,  I 
also  should  like  to  tell  you  that,  in  order  to  discuss  this  matter  more  fully, 
with  complete  knowledge  of  what  exactly  goes  into  and  comes  out  of  the  black 
box,  a  classified  conference  is.being  planned  during  the  summer  at  Los  Alamos. 
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I  hope  it  will  be  very  successful  and  produce  new  scientific  ideas  and  new  sci¬ 
entific  approaches. 

There  is  another  and  very  practical  field  which  also  has  been  neglected  — 
at  least  neglected  in  comparison  with  its  enormous  potential.  One  of  our  good 
friends,  John  Grebe,  will  talk  to  you  about  this  field:  the  use  of  a  nuclear 
explosion  underground  to  produce  the  heat  and  pressure  to  induce  chemical 
reactions. 

I  should  like  to  start  this  topic  in  a  somewhat  lighter  vein  by  mentioning 
to  you  a  relatively  small  and  interesting  object  which  actually  commands  more 
interest  than  is  justified.  I  have  diamonds  in  mind.  Now,  I  don8t  think  that 
diamonds  are  important,  except  in  some  industrial  processes,  and  if  we  ever 
succeed  in  making  diamonds,  all  we  are  going  to  accomplish  is  to  ruin  the 
diamond  market.  The  importance  of  diamonds  cannot  compare  with  the  im¬ 
portance  of  oil  and  water.  But,  just  for  the  fun  of  it,  let  me  talk  for  a  moment 
about  diamonds. 

We  happen  to  know  diamonds  can  be  produced  if  we  have  three  things  — 
high  pressure,  low  temperature,  and  time;  then  diamonds  will  grow  out  of  pure 
carbon.  Now,  high  pressure  from  nuclear  explosives  we  have,  but  low  tem¬ 
peratures  we  do  not.  Perhaps  if  we  are  very  clever  about  it  we  might  be  able 
to  arrange  conditions  so  as  to  give  in  certain  regions  high  pressure  without 
excessive  temperature.  This  much  I  think  could  be  done.  What  we  do  not 
have  in  nuclear  explosions  is  time.  We  cannot  maintain  the  high  pressure  for 
a  long  period.  This  should  not  necessarily  be  discouraging.  Maybe  time  can 
be  replaced  by  exceedingly  high  pressure,  and  perhaps  we  have  a  chance  to 
produce  this  exceedingly  high  pressure.  Even  so,  I  don’t  have  a  hope  for  the 
Hope  diamond.  What  we  are  much  more  likely  to  produce,  lest  the  diamond 
industry  react  to  my  remarks,  is  very  fine  and  very  imperfect  diamond  powder. 

All  this  I  am  using  only  as  an  example.  Maybe  we  can  initiate  geochemical 
and  geophysical  transformations  which  usually  go  on  deep  inside  the  earth  and 
reproduce  types  of  materials  never  seen  near  the  surface  of  the  earth. 

We  can  produce  a  number  of  chemical  reactions  in  the  strange  reactors 
which  we  will  have  underground.  Provided  we  can  make  the  products  either 
volatile  or  fluid,  and  some  of  them  will  be,  it  will  not  be  very  difficult  to  get 
them  to  the  surface,  and  our  chemical  industry  might  be  seriously  interested. 
Now  this  kind  of  thing  cannot  be  done  without  the  cooperation  of  chemists, 
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chemical  engineers,  and  economists.  In  a  laboratory  whose  aims  are  really 
a  little  bit  different,  we  cannot  do  much  more  than  give  you  the  means,  draw 
the  opportunities  to  your  attention,  and  then  hope  for  the  best. 

This  group  here,  if  you'll  pardon  me  for  saying  so,  seems  to  be  similar 
to  a  bunch  of  starved  and  blind  people  stumbling  around  in  an  extremely  rich 
garden.  The  food,  I'm  sure,  is  there  to  pick;  the  only  question  you  have  to 
decide  is  in  which  direction  you  want  to  extend  your  hand;  then  you  have  to 
rely  on  the  best  methods  of  touch,  smell,  and  intelligence  lest  you  bite  into 
a  bitter  fruit.  I  am  sure  that  the  opportunity  is  here.  I'm  sure  that  some¬ 
where  among  you  there  is  the  idea,  the  initiative,  the  perseverance  and  the 
genius  which  will  bring  this  Plowshare  program  to  fruitation. 
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B.  The  Data  from  Underground 
Nuclear  Explosions 

1.  SOME  BASIC  DATA  ON  UNDERGROUND 
NUCLEAR  EXPLOSIONS 

Gerald  W.  Johnson 

Lawrence  Radiation  Laboratory,  Livermore 

General  Description  of  an  Underground  Explosion 

In  discussing  the  phenomenology  of  underground  explosions,  it  is  useful 

to  consider  several  different  time  scales,  with  their  associated  phenomena, 

12 

The  Rainier  event,  in  which  10  7  kilotons  of  energy  (la  7  X  10  calories, 

6.  7  X  10^  Btu,  or  2.  0  X  10^  kw-hr)  were  released  in  a  room  6  feet  square, 

7  feet  high,  and  800  feet  below  the  surface  of  the  earth,  will  be  used  as  an 
example. 

a.  Nuclear  Phase  (Microseconds).  The  energy  of  the  nuclear  explosive 
is  generated  in  a  few  tenths  of  a  microsecond,  vaporizing  the  assembly  ma¬ 
terials  and  forming  a  rapidly  growing  fireball.  The  material  in  the  Rainier 
room  reached  a  temperature  of  ten  million  degrees  Fahrenheit  and  a  pressure 
of  seven  million  atmospheres. 

b.  Hydrodynamic  Phase  (Milliseconds),  A  shock  wave  proceeds  outward 
vaporizing,  melting,  and  crushing  the  surrounding  medium.  A  cavity  is 
formed  by  the  outward  motion  of  the  rock  near  the  center  of  the  explosion. 

In  Rainier,  the  shock  vaporized  a  3-foot  thickness  of  rock  (about  100  tons), 
melted  7  more  feet  (about  600  tons),  and  crushed  rock  out  to  130  feet.  The 
cavity  in  Rainier  reached  its  final  radius  of  55  feet  in  about  a  hundred  milli¬ 
seconds. 

c.  Quasi- static  Phase  (Seconds  to  Minutes).  The  cavity  cools  and  col¬ 
lapses  when  the  internal  pressure  can  no  longer  support  it.  In  Rainier,  the 
pressure  dropped  to  forty  atmospheres  in  about  a  minute.  A  fissure  must 
then  have  opened,  allowing  the  gases  to  escape  and  releasing  the  internal 
pressure;  the  roof  caved  in  and  a  chimney  of  rubble  was  formed.  About  two 
hundred  thousand  tons  of  broken  rock  were  produced. 

d.  Long-term  Phase  (Days  to  Years).  A  slow  transport  of  heat  through 
the  rock  takes  place,  and  the  radioactive  products  decay.  In  Rainier,  the 
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water  present  in  the  rock  (20%  by  weight)  vaporized  and  distributed  the  heat 
over  a  large  volume;  the  highest  temperature  now  present  is  at  the  boiling 
point  of  water. 

Containment  and  Cratering  Data 

Nine  past  explosions  will  be  of  particular  interest  in  this  symposium. 
Several  were  near  the  surface  and  were  thus  not  contained;  they  gave  useful 
data  for  cratering  studies.  Several  were  completely  contained,  producing  no 
craters  and  releasing  no  radioactivity  to  the  atmosphere.  Table  I  summarizes 
the  yield,  cratering,  and  released  radioactivity  data.  The  cratering  informa¬ 
tion  will  be  dealt  with  in  more  detail  in  subsequent  papers.  The  containment 
data  may  be  summarized  by  saying  that  for  depths  of  burial  greater  than 
D  -  450  W  ^  (where  W  is  the  yield  in  kilotons)  the  radioactivity  will  be  com¬ 
pletely  confined. 

Geological  Data 

The  six  explosions  which  were  contained  or  nearly  contained  took  place 
in  a  geological  formation  known  as  Oak  Springs  tuff.  This  formation,  which 
is  some  1900  feet  in  thickness  in  the  tunnel  area,  has  been  divided  into  eight 
lithologic  units,  TOS^  through  TOSg.  The  rocks  composing  these  units  are 
fine  to  coarse,  grained-bedded  to  massive  tuffs,  welded  tuff  tuffaceous  sand¬ 
stones,  tuff  breccias,  and  agglomerates.  The  rocks  of  the  Oak  Springs  for¬ 
mation  are  siliceous  and  contain  from  65  to  75%  SiO^  Vesicles  average  6% 
of  the  volume  of  the  nonwelded  tuffs,  but  only  3.  5%  of  the  volume  of  welded 
tuffs. 

The  tuffs  can  be  divided  into  three  groups:  welded  tuff,  units  TOS^  and 
TOSg;  friable  tuff,  unit  TOS?;  and  bedded  tuff,  units  TOSj  through  TOSg. 

Welded  tuff  forms  the  cap  rock  of  the  mesa  (unit  TOSQ).  The  color,  of 

o 

welded  tuff  ranges  from  reddish  purple  to  gray.  Phenocrysts  of  quartz, 
feldspar,  and  biotite  are  common.  The  rock  is  hard  and  dense. 

The  rocks  of  the  friable  tuff,  unit  TOSy,  are  weak.  Samples  can  easily 
be  broken  by  hand,  and  drill  core  recovery  is  poor.  Most  beds  are  light 
gray  or  brown  and  consist  of  granules  of  lapilli  and  fragments  of  quartz  and 
feldspar.  Biotite  is  relatively  conspicuous;  fragments  of  quartzite  and 
welded  tuff  are  locally  abundant.  It  was  in  this  zone  that  the  Rainier  event 
took  place. 
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eakthrough  to  surface  but  stemming  failed,  releasing  radioactivity  into  the  tunnel. 
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Bedded  tuff  comprises  units  TOS^  through  TOS^  They  are  not  welded 
but  are  composed  of  zeolitized  ash  shards  and  pumice,  fragments  of  welded 
tuff  and  quartzite,  and  phenocrysts  of  quartz,  feldspar,  and  biotite,  all  ce¬ 
mented  together  by  zeolite  minerals.  The  unit  is  well  bedded  with  individual 
distinctive  beds  ranging  from  a  few  inches  to  a  few  tens  of  feet  in  thickness. 
Most  of  the  beds  are  nearly  white  or  reddish  purple.  Pumice  fragments  are 
abundant  in  most  beds,  and  in  parts  of  unit  TOSj-  they  comprise  a  large  part 
of  the  rock.  The  Blanca  event  occurred  in  unit  TOS^  of  this  group. 

The  "basement”  rock  of  the  formation  is  carbonate,  with  limestone  as 
the  dominating  type.  Figure  1.  1  shows  the  bedding  structure  for  the  Rainier 
and  Blanca  shots. 
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2.  THE  EARLY  CAVITY  HISTORY  OF  RAINIER 
AND  COMMENTS  ON  ENERGY  STORAGE 

George  C.  Kennedy 

Institute  of  Geophysics,  University  of  California*  Los  Angeles 

There  are  three  major  chemical  environments  in  which  underground 
shots  can  be  made.  The  knowledge  of  explosions  in  these  environments  is 
important  if  they  are  to  be  used  for  the  underground  storage  of  heat.  The 
behavior  of  these  three  systems,  which  are  (1)  silicate  (lava,  granite,  slate), 
(2)  limestone  (CaCOs),  and  (3)  salt  (NaCI),  is  grossly  contrasting.  Since  Dr. 
Higgins  and  I  have  worked  out  in  some  detail  the  phenomenology  of  a  shot  in 
silicate1  (the  tuff  in  the  Rainier  event),  I  would  like  to  review  that  first.  A 
surprisingly  detailed  history  of  temperatures  and  pressures  in  the  under¬ 
ground  cavity  produced  by  the  Rainier  Event  can  be  worked  out  by  physical 
examination  of  the  sarhples  obtained  by  drilling  and  tunneling  in  the  shot  vi¬ 
cinity,  and  by  radiochemical  analysis  of  these  samples. 

Several  milliseconds  after  detonation  of  the  device,  a  cavity  was  formed, 
lined  with  approximately  10  cm  of  fused  rock  at  high  temperatures.  The  total 
quantity  of  fused  rock  has  been  estimated  from  the  concentrations  of  radio¬ 
active  products  contained  in  various  samples  to  be  approximately  700  metric 
tons.  The  gradient  of  the  temperature  between  the  shell  of  fused  rock  and  the 
surrounding  unfused  rock  was  very  steep.  Rhyolite  glass  is  exceedingly  vis¬ 
cous  up  to  temperatures  as  high  as  1400°C  and  flows  very  slowly  at  lower 
temperatures.  However,  numerous  thin,  small  fragments  showing  contact 
between  unfused  rock  and  molten  rock  were  recovered  in  drill  cores.  These 
fragments  show  a  maximum  of  l/4  inch  of  glass  adhering  to  unfused  rock. 
Through  that  l/4  inch  the  temperature  difference  must  have  been  on  the  order 
of  300~400°C.  Unaltered  tuff  fuses  at  about  1100°C,  which  may  serve  as  an 
approximate  boundary  temperature  between  the  glass  lining  and  the  unfused 
rhyolite.  However,  l/4  inch  in  toward  the  shot  chamber  from  the  fused  mar¬ 
gin,  the  rhyolite  was  sufficiently  hot  to  be  fluid  and  fall  as  minute  droplets, 
which  suggests  that  temperatures  on  the  inside  of  this  l/4  inch  rim  of  glass 
may  have  been  on  the  order  of  1500°C. 

1G.  C.  Kennedy  and  G.  H.  Higgins,  Temperatures  and  Pressures  Asso¬ 
ciated  with  the  Cavity  Produced  by  the  Rainier  Event,  UCRL-5281,  July  1958. 
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Numerous  small  droplets,  beads,  and  shot-like  particles  of  glass  were 
found  mixed  with  the  collapsed  rubble  recovered  from  the  bottom  of  the  cav¬ 
ity  (see  Fig.  2  0  1).  These  beads  apparently  are  droplets  of  glass  that  dripped 
from  the  roof  and  walls  of  the  chamber  early  in  the  sequence  of  events  follow- 
ing  the  explosion.  An  interesting  phenomenon  appears  upon  close  examina¬ 
tion  of  these  droplets  of  glass.  Although  they  were  obviously  once  quite  fluid, 
there  is  no  evidence  of  indentation  as  would  be  expected  if  they  fell  from  the 
roof  of  the  chamber  and  reached  the  floor  while  still  fluid.  Presumably  dur- 
ing  the  several  seconds  of  fall  from  the  roof  of  the  chamber  to  the  floor  of  the 
chamber,  these  droplets  changed  in  temperature  very  rapidly  and  on  impact 
with  the  bottom  were  cool  enough  so  that  the  glass  was  quite  stiff  if  not  solid. 
Exceedingly  rapid  and  early  cooling  apparently  occurred  within  the  shot 
chamber. 

An  estimate  of  the  vapor  pressure  of  water  in  the  chamber  at  the  time 
these  glass  droplets  froze  may  be  obtained.  The  glass  beads  vesiculate 
(that  is,  become  filled  with  small  cavities)  excessively  and  fuse  at  around 
920°C.  They  vesiculate  due  to  evolution  of  water  from  the  glass  and  from 
the  very  minute  bubbles  within  the  glass.  A  few  of  these  beads  were  weighed, 
heated  to  outgas  them,  and  reweighed.  A  sample  of  vesiculated  material 
recovered  from  the  site  and  a  sample  of  the  vesiculated  material  produced 
by  heating  the  bead  are  shown  in  Fig.  2.2.  About  40  atmospheres  of  vapor 

pressure  at  1200-1500°C  in  the  chamber  is  required  to  account  for  all  the 
water  in  the  glass  and  the  bubbles,  assuming  of  course  that  the  glass  beads 
were  in  equilibrium  with  the  vapor  inside  the  chamber. 

It  is  clear,  however,  that  though  these  droplets  were  quenched  during 
fall  at  a  pressure  of  approximately  40  atmospheres,  the  pressure  of  water 
vapor  in  the  shot  chamber  was  dropping  rapidly.  A  puddle  of  glass  (several 
inches  thick,  as  sampled  by  core  drilling)  had  earlier  accumulated  on  the 
floor  of  the  chamber  by  slumping  of  the  10  cm  glass  lining  of  the  chamber. 

This  puddle  was  saturated  with  water;  but  before  it  could  freeze,  the  water 
pressure  in  the  shot  chamber  had  dropped  sufficiently  to  make  the  puddle 
vesiculate  extensively  as  it  outgassed  (see  Fig.  2.  2).  Thus  the  evidence 
strongly  points  to  the  fact  that  though  the  vapor  pressure  was  40  atmospheres 
inside  the  chamber  immediately  after  the  shot  while  the  glass  droplets  were 
still  falling  from  the  roof  of  the  chamber,  a  minute  or  two  later  the  pressure 
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had  dropped  to  only  a  few  atmospheres.  A  heating  experiment  on  the  vesic- 
ulated  glass  from  the  bottom  of  the  chamber  indicated  that  it  had  been  nearly 
completely  outgassed.  This  particular  material  does  not  vesiculate  on  heat¬ 
ing  as  do: the  glass  droplets. 

Two  more  observations  concerning  the  small  droplets  of  glass  which 
solidified  as  beads  early  in  the  history  of  the  shot  throw  some  light  on  the 
sequence  of  events  inside  the  shot  chamber.  Most  of  the  droplets  of  quenched 
glass  show  minute  bumps  and  pimples  all  over  their  surfaces.  These  pimples 
are  apparently  the  result  of  other  exceedingly  small  droplets  of  glass  collect¬ 
ing  and  adhering  to  the  larger  droplets  during  cooling.  They  suggest  that 
the  cavity  contained  numerous  very  minute  dust  particles  which  probably  con¬ 
densed  from  the  vapor  phase  and,  because  of  their  small  size,  cooled  very 
quickly  and  adhered  to  various  glass  beads  and  other  free  surfaces.  Two 
glass  beads,  l/2  and  3/4  of  an  inch  in  diameter,  do  show  slight  flattening  on 
the  bottom  as  though  they  had  not  entirely  chilled  during  the  time  of  fall  from 
the  roof  to  the  bottom.  All  those  smaller  than  l/2  inch  show  no  sign  of  flat¬ 
tening  or  no  sign  of  being  molten  on  reaching  the  bottom  of  the  chamber. 

Thus  we  can  reconstruct  the  pressure-temperature  history  of  the  cham¬ 
ber,  starting  with  an  initial  water-vapor  pressure  based  on  the  assumption 

that  approximately  700  metric  tons  of  tuff  containing  something  on  the  order 

8 

of  15^c  water  were  melted,  vaporizing  approximately  10  grams  of  water. 

This  water  vapor  filled  the  cavity,  whose  volume  was  approximately  2  X  10^ 

3 

cm  .  Thus  the  specific  volume  of  water  in  the  chamber  was  approximately 
3 

200  cm  /g,  corresponding  to  a  pressure  of  almost  exactly  40  atmospheres  at 
a  temperature  of  1500°C.  The  agreement  between  this  computed  pressure, 
based  on  the  estimated  amount  of  rock  fused  and  its  initial  water  concentra¬ 
tion,  with  the  pressure  based  on  measurements  of  the  amount  of  water  in  the 
quenched  glass  droplets  is  remarkable.  Because  the  thin  puddle  of  glass  that 
formed  early  on  the  bottom  of  the  chamber  soon  vesiculated,  we  can  deduce 
that  the  pressure  dropped  rapidly  after  the  shot.  Coincident  with  the  pressure 
drop  was  a  very  sharp  and  marked  cooling  as  evidenced  by  the  glass  beads 
quenched  during  fall.  It  seems  quite  likely  that  the  extremely  rapid  cooling 
was  due  to  the  expansion  of  the  gas  in  the  chamber. 

Radiometric  analyses  of  the  unfused  rock  immediately  adjacent  to  the 
molten  shell  strongly  suggest  that  the  water  vapor  and  radioactive  gases  in 
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the  chamber  did  not  escape  by  a  general  outward  diffusion.  Levels  of  radio¬ 
activity  in  the  unfused  rock  adjacent  to  the  shell  were  not  detectable  above 
background.  Since  an  appreciable  quantity  of  gaseous  radioactive  products 
resulted  from  the  explosion,  it  is  clear  that  the  highly  compacted  tuff  adjacent 
to  the  walls  of  the  chamber  was  essentially  impermeable  to  gases. 

Evidence  of  the  vesiculated  melt  at  the  bottom  of  the  chamber,  the 
quenched  beads  of  glass,  the  marked  loss  of  such  short-lived  volatile  radio¬ 
active  products  as  krypton  and  xenon  from  the  chamber,  and  the  marked 
asymmetry  of  the  post- shot  thermal-energy  distribution  strongly  suggests 
the  chamber  vented  through  one  or  at  the  most  a  few  major  fractures.  Appar¬ 
ently  none  of  these  fractures  communicated  with  the  surface  as  no  activity 
was  observed  on  the  surface  of  the  ground  after  the  shot.  One  of  these  escape 
channels  was  probably  intersected  by  a  core  drill  hole  75  to  80  feet  beyond 
the  shot  site  and  diametrically  opposite  to  the  access  tunnel.  The  nonvola¬ 
tile  decay  products  of  the  short-lived  radioactive  gases  were  highly  concen¬ 
trated  in  the  fused  rock  along  this  channel. 

Since  the  decay  of  the  gaseous  products  into  nonvolatile  elements  is  known 
to  be  exponential  with  time,  the  percentage  loss  of  volatile  precursers  is  a 
direct  indication  of  the  time  of  escape.  Times  of  escape  have  been  computed 
for  a  number  of  elements  and  a  number  of  analyses.  The  range  of  times  of 
escape  computed  by  this  method  is  remarkably  narrow  and  suggests  that  the 
chamber  vented  between  30  seconds  and  two  minutes  after  the  shot.  The  time 
of  venting  derived  from  examination  of  the  glasses  in  the  chamber  and  the 
time  derived  from  radiochemical  sampling  are  in  close  agreement. 

Coincident  with  the  accumulation  of  the  puddle  of  glass  on  the  floor  of  the 
chamber  and  with  the  venting  during  the  first  minute  or  two  after  the  shot, 
other  phenomena  were  taking  place  in  the  cavity.  The  dense,  highly  com¬ 
pressed  tuff  of  the  walls,  with  a  water  content  of  approximately  15%,  was 
steadily  increasing  in  temperature  owing  to  conduction  of  heat  outward  from 
the  adjacent  thin,  molten  shell.  A  series  of  implosions  of  the  shell  into  the 
cavity  was  probably  due  in  large  part  to  expansions  of  water  in  the  dense  tuff. 
Nowhere  in  drill  holes  or  tunnel  opening  has  the  original  glass-lined  wall 
been  seen.  This  wall  material  was  all  imploded  into  the  cavity.  These  im¬ 
plosions  continued  over  a  considerable  time.  They  began  before  all  the  melt 
had  accumulated  on  the  bottom  of  the  cavity  and  continued  until  sometime  after 
the  cavity  had  vented  and  the  general  collapse  of  the  roof  had  begun. 
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It  is  interesting  to  compare  the  initial  gas  pressure  in  the  cavity  with  the 
weight  of  the  column  of  rock  overlying  the  cavity.  If  an  overburden  of  850 
feet  of  tuff  with  a  density  of  approximately  2,  0  g/cc  is  assumed,  the  rock 
pressure  on  the  cavity  is  50  to  55  atmospheres,  which  is  counterbalanced  by 
the  initial  gas  pressure  in  the  cavity  of  approximately  40  atmospheres.  Thus 
it  seems  likely  that  very  shortly  after  the  gas  escaped  from  the  cavity  and 
the  pressure  dropped  to  a  few  atmospheres,  the  roof  would  collapse.  Most 
of  this  roof  collapse  did  take  place  after  the  cavity  had  cooled  and  the  gas 
had  escaped.  The  total  absence  of  fragments  of  rock  in  the  puddle  of  molten 
lava  at  the  base  of  the  cavity  shows  that  all  the  lava  had  vesiculated  and 
cooled  before  collapse  of  the  cavity  took  place, 

NaCl  as  a  chemical  environment  for  the  storage  of  heat  will  be  studied 
in  the  Gnome  experiment.  Surprisingly  enough  the  NaChH^O  system  is  air- 
most  unexplored,  but  recently  we  have  made  a  large  number  of  measurements 
on  its  chemistry,  I  will  describe  these  measurements  later  in  this  symposium 
(see  p,  4  of  UCRL-5678,  Part  IV  of  the  Symposium  proceedings),  and  will 
only  summarize  a  few  of  the  details  now, 

A  most  unexpected  result  was  the  discovery  that  the  concentration  of  salt 
in  steam  is  almost  independent  of  temperature  up  to  600°C.  Further,  at  high 
pressures,  the  concentration  of  salt  in  steam  decreases  as  the  temperature 
rises  toward  600°C;  there  is  a  highly  optimum  temperature  in  the  system  at 
which  there  is  a  low  salt  content  in  the  steam,  amounting  to  around  200-300 
parts  per  million. 

The  least  explored  system  today  is  that  of  CaCO^0  H^O,  and  I  think  par¬ 
ticularly  serious  consideration  should  be  given  to  its  detailed  examination, 

A  casual  glance  suggests  that  it  is  one  of  the  most  interesting  of  all  systems 
in  which  to  store  energy  because  of  its  self- sealing  nature  and  because  of  the 
high  fluidity  of  the  molten  phase.  Solubility  of  CaCO^  in  gases  is  small,  and 
limestone  environments  are  geologically  abundant. 

Suppose  we  examine  the  details  of  an  explosion  in  limestone.  Immediately 
after  the  shot,  CaCO^  is  melted  at  a  temperature  of  approximately  1220°C 
and  a  pressure  of  40  atmospheres  of  CC^,  However,  in  the  CaCO^0  CaO 
system,  there  is  a  minimum  in  the  melting  point  curve. 
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1200°C 


Consequently  as  soon  as  a  little  CO2  is  formed,  CaO  will  form  in  equili¬ 
brium  with  CaCO^.  As  some  of  the  CO2  gas  streams  off,  the  melting  tempera¬ 
ture  will  drop  sharply  to  around  900°C;  there  is  the  possibility  of  having  a  large 
fluid  body  of  a  low  temperature  melt  in  the  system  CaCO^"  CaO  with  relatively 
low  CO^  vapor  pressure,  because  of  the  shape  of  the  vapor  pressure  vs  tem¬ 
perature  curve  in  a  reaction  of  this  kind. 


Furthermore,  the  addition  of  water  at  1000  atmospheres  will  probably 
lower  the  temperature  approximately  another  200 °C.  Therefore,  there  is  a 
good  probability  of  obtaining  highly  fluid  melts  at  temperatures  of  700-800°C. 

Limestone  beds  thus  may  make  effective  reservoirs  for  the  storage  of  heat. 
In  addition,  they  have  certain  properties  that  NaCl  beds  do  not  have.  For  in¬ 
stance,  they  are  self-sealing  in  case  of  a  crack.  If  we  have  a  ball  of  melt 


-27- 


UCRL-5675 


underground  in  the  system  CaCO^'H^O,  and  a  crack  opens  up  in  the  surround¬ 
ing  rock  and  melt  starts  leaking  away,  the  lowering  of  the  water  pressure  due 
to  escape  of  gas  will  raise  the  melting  point  of  the  carbonate.  The  melt  in  the 
crack  will  becomes  solid  and  seal  off  any  further  escape  of  the  melt.  The  sys¬ 
tem  thus  tends  to  be  self- sealing  like  a  bullet  proof  gas  tank. 

To  summarize  this  talk,  I  wish  to  emphasize  that  the  phenomenology  of 
underground  nuclear  explosions  in  the  three  most  likely  chemical  environments 
is  greatly  contrasting.  Particularly  serious  attention,  I  think,  should  be  given 
to  the  detailed  examination  of  CaCO^H^O  systems,  as  this  is  the  system  which 
is  the  most  interesting  and  which  is  least  explored  today  in  terms  of  physical 


behavior. 
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3.  CLOSE-IN  SHOCK  TIME-OF-ARRIVAL  MEASUREMENTS 
AND  HYDRODYNAMIC  YIELD1 

Fo  B.  Porzel 

Senior  Scientific  Advisor 

Armour  Research  Foundation  of  the  Illinois  Institute  of  Technology 

Part  A 


This  paper  briefly  describes  the  principal  ideas  in  the  theory,  experi¬ 
mental  techniques  and  results  for  the  programs  entitled,  nClose-In  Shock 
Time-of- Arrival  Measurements, 11  which  were  performed  by  Armour  Research 
Foundation  (ARF)  on  Operations  Rainier  (1957)  and  Hardtack  II  (1958)  under 
contract  to  the  University  of  California  Lawrence  Radiation  Laboratory.  The 
end  purpose  of  the  ARF  programs  was  to  determine  the  hydrodynamic  yield; 
they  exploited  methods  previously  used  by  ARF  for  similar  purposes  for  the 
Navy  in  Operations  Wigwam  (1955)  and  Hardtack  I  (1958),  and  are  based  in 
turn  on  similar  measurements  by  Group  J-10  of  the  Los  Alamos  Scientific 
Laboratory  on  Operation  Jangle  (1951)  at  the  Nevada  Test  Site.  Another  im¬ 
mediate  and  more  basic  purpose  of  these  measurements  is  to  provide  better 
understanding  of  the  phenomenology  of  underground  explosions. 

This  paper  covers  four  topics:  experimental  techniques,  theory,  Rainier 
results,  and  Hardtack  II  results. 


This  paper  is  in  two  parts.  Part  A  is  primarily  concerned  with  close-in 
shock  time-of- arrival  measurements  for  the  determination  of  hydrodynamic 
yield,  on  operations  Rainier  and  Hardtack  II.  This  part  is  essentially  the  ver¬ 
bal  presentation  made  by  the  author  at  the  Plowshare  Symposium  in  May,  1959 
at  San  Francisco.  The  test  report  for  Operation  Rainier  covers  this  section 
in  detail:  WT-1495,  Project  25.  1,  "Close-In  Time-of- Arrival  Measurements 
for  Yield  of  Underground  Rainier  Shot,n  University  of  California,  Lawrence 
Radiation  Laboratory. 

Part  B  is  primarily  concerned  with  the  theoretical  concepts  used  in  Part 
A,  and  in  particular,  discusses  some  implications  of  the  theory  for  these 
measurements  which  apply  to  control  and  use  of  explosions,  aside  from  the 
specific  test  operations  described  in  Part  A.  [Editor's  Note:  Part  B  proved 
too  long  to  publish  as  part  of  the  Plowshare  Symposium.  Sections  I,  II,  III, 

VII*  and  VIII,  which  are  directly  referred  to  in  Part.  A,  have  been  included. 

The  other  sections  can  be  obtained  in  Armour  Research  Foundation  Technical 
Memorandum  421,  nNew  Concepts  for  Control  and  Use  of  Nuclear  Explosions,11 
F.  B,  Porzel.  ] 


-29- 


UCRL-5675 


I.  Experimental  Techniques 

Two  principal  experimental  techniques  have  been  used:  time -of- arrival 
switches  and  Doppler.  By  far  the  greatest  reliance  so  far  has  been  placed  on 
the  switches  and  the  experimental  arrangement  for  Rainier  is  shown  schemat¬ 
ically  in  Fig.  3.  1.  The  switches  are  placed  at  essentially  equal  logarithmic 
spacings  from  the  nuclear  device  out  to  the  expected  limit  of  crushing,  beyond 
which  distance  the  shock  velocity  becomes  sonic  and  hence  ambiguous  for  di¬ 
agnosis  of  hydrodynamic  effects.  The  blast  switch  itself  is  essentially  an  open 
terminal  of  a  coaxial  cable  whose  inner  conductor  is  closed  by  crushing  of  an 
outer  conductor  upon  it  when  the  main  shock  pulse  arrives.  The  switch  is 
designed  not  to  close  from  weak  precursor  sound  signals.  The  signals  from 
a  single  string  of  switches  are  electronically  coded  for  pulse  length  and  pulse 
width  in  a  "mixing  box,M  then  transmitted  on  a  single  cable  to  a  shelter  out¬ 
side  the  tunnel,  where  they  are  displayed  on  an  oscilloscope  at  a  20,  000  cps 
raster  sweep  and  recorded  on  a  GR  No,  651-AH  streak  camera.  Time  is 
measured  with  almost  arbitrary  accuracy  because  the  abrupt  pulse  rise  on 
the  raster  can  be  read  to  one  microsecond  whereas  the  hydrodynamic  events 
occur  on  a  millisecond  time  scale.  Also,  because  the  velocities  are  deter¬ 
mined  by  differentiation  of  the  time-of-arrival  curve,  good  relative  positions 
of  the  gauges  to  its  neighbors  is  sufficient  for  velocity  and  to  check  on  the 
absolute  position  through  R  =  ^XJ  dt. 

The  essential  idea  in  the  Doppler  system  is  to  place  a  hollow  coaxial  cable 
on  a  radial  line  toward  the  device.  An  rf  signal  of  1580  megacycles  is  input 
and  forms  a  standing  wave.  As  the  shock  advances  and  crushes  the  outer  con¬ 
ductor  upon  the  inner  conductor,  the  standing  wave  shortens;  the  reflected 
wave  is  heterodyned  with  the  input,  and  the  Doppler  frequency  of  the  mixed 
signal  is  displayed  on  an  oscilloscope.  The  Doppler  frequency  is^ directly 
proportional  to  the  shock  velocity.  The  Doppler  technique  is  greatly  compli¬ 
cated  by  a  "pinch  jet,"  described  in  Section  VII  of  Part  B  of  this  paper:  the 
input  signal  is  reflected  from  a  shock  in  pure  air  induced  ahead  of  the  crushed 
cable.  But  even  without  special  provision  for  jetting,  the  Doppler  can  provide 
an  accurate  value  of  shock  velocity  near  the  walls  of  the  device  chamber. 

The  blast  switches  and  Doppler  system  offer  an  ideal  combination  for 
mutual  check  and  support  because  the  switches  yield  the  integrated  time-of- 
arrival  for  gross  phenomenology  and  the  Doppler  system  provides  great 
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Fig.  3.  1.  Rainier  equipment  layout. 
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detail  in  the  differentiated  counterpart  as  shock  velocity  versus  time. 

II.  Theory 

As  is  well  knoyn,  given  the  Rankin e- Hug oniot  adiabat,  the  measurement 
of  the  shock  velocity  permits  the  direct  calculation  of  peak  pressure,  density, 
material  velocity,  and  other  hydrodynamic  variables  of  interest  at  the  peak 
of  the  main  pressure  pulse.  It  is  important  for  the  case  of  soil  and  rock  to 
recall  that  these  relations  require  only  gross  conservation  of  mass,  momen¬ 
tum,  and  energy  across  a  narrow  domain  and  are  not  sensitive  to  the  exist¬ 
ence  of  a  truly  sharp  shock  front.  By  means  of  a  derivation  the  author  calls 
the  "matrix  formulation,"  it  is  then  possible  to  derive  the  space  gradients  and 
time  rate  of  decay  on  the  interior  of  the  wave  from  the  measured  peak  values 
at  the  front  and  from  their  measured  decay  with  distance  and  time.  Finally 
by  means  of  the  analytic  solution  —  being  the  repeated  use  of  conservation  of 
mass,  momentum,  and  energy  on  the  inferior  of  the  wave  —  the  complete 
waveforms  for  all  the  hydrodynamic  variables  on  the  interior  can  be  calcu¬ 
lated  as  a  function  of  space  and  time,  and  the  hydrodynamic  energy  can  be 
evaluated  by  integrating  the  energy  waveforms. 

The  waste  heat  concept  plays  a  controlling  role  in  the  description  of  ex¬ 
plosions  in  solids  and  liquids;  this  concept  is  discussed  in  Section  I  of  Part  B, 
(Also  in  Part  B  are  discussed:  in  Section  II,  the  general  characteristics  of 
wave  propagation;  and  in  Section  III,  some  implications  regarding  underground 
explosions.  )  The  waste  heat  concept  greatly  simplifies  the  description  for  the 
behavior  of  the  material  on  the  interior  of  the  wave  during  its  expansion  behind 
the  shock,  and  the  validity  of  the  concept  rests  in  good  part  on  considerations 
of  "grain  size"  of  the  mixture.  Characteristically,  strong  shocks  produce 
high  temperatures  in  a  gas  phase  but  usually  produce  low  temperature  in  the 
solid  phase;  simple  considerations  of  the  solutions  to  the  heat  conduction  equa¬ 
tion  will  show,  as  an  example,  that  for  particles  approximately  0.  1  milli¬ 
meter  in  diameter,  the  heat  flows  from  the  gas  phase  into  the  solid  phase  in 
times  like  0.  1  millisecond.  Such  considerations  justify  the  assumption  of 
temperature  equilibrium  for  the  time  scale  of  a  nuclear  explosion  in  soil,  but 
it  is  to  be  expected  that  this  assumption  is  not  justified  for  the  shorter  time 
scale  of  explosions  induced  in  soil  from  small  charges. 

The  equation  of  state  for  Nevada  soil  along  the  Rankine-Hugoniot  adiabat 
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has  been  experimentally  determined  by  LRL,  furnished  to  Armour  Research 
Foundation,  and  will  not  be  described  here.  Oddly  enough,  despite  the  com¬ 
plexity  of  the  thermodynamics  for  rock,  especially  soils  with  high  water  and 
air  content,  elaborate  thermodynamic  detail  is  not  required  in  the  equation 
of  state.  The  dynamic  compressibility  is  largely  controlled  by  the  air  con¬ 
tent. and  controlled  at  somewhat  higher  pressures  by  the  compressibility  of 
the  water  fraction.  A  quite  adequate  equation  of  state  can  be  constructed 
from  the  knowledge  of  the  rock  density  in  situ,  the  void  fractions  (essentially 
from  the  comparison  of  actual  density  to  crystal  density),  the  water  content, 
and  the  gross  mineralogical  composition.  We  require,  in  addition,  a  speci¬ 
fication  of  the  waste  heat  up  to  the  boiling  point  of  water  and  up  to  effective 
vaporization  of  the  soil.  The  latter  quantities  are  rather  easily  estimated 
from  the  boiling  temperature  of  water,  the  water  content  and  its  heat  of  va¬ 
porization.  Surprisingly  enough,  Nevada  soil  is  close  to  saturation  in  water 
and,  except  on  the  deep  interior  of  the  cavity,  temperatures  exceeding  the 
local  boiling  point  of  water  are  quite  unlikely.  The  ARF  studies  of  these 
phenomena  led  to  postulating  the  concept  "decrepitation";  when  the  waste  heat 
is  sufficiently  high  to  vaporize  irreversibly  the  water  fraction  and  waters  of 
crystalization,  the  final  state  of  the  rock  could  only  appear  as  exploded  ma¬ 
terial,  similar  to  the  process  of  making  puffed  cereals  and  like  the  familiar 
phenomena  of  decrepitation  in  certain  minerals.  Because  the  rock  samples 
obtained  to  determine  the  equation  of  state  for  Rainier  had  evidently  dried  out 
before  testing,  the  actual  void  ratio  for  Rainier  was  calculated  using  the  ob¬ 
served  crushing  and  cavity  radii.  We  assume  (1)  the  gross  hydrodynamic  dis 
placements  are  negligible  beyond  the  crushing  radius  as  far  as  they  affect  the 
cavity  radius,  and  (2)  the  solid  fraction  of  material  on  the  interior  was  com¬ 
pacted  in  a  spherical  shell  of  inner  radius  essentially  the  cavity  radius  x, 
when  the  shock  reached  the  crushing  radius  R.  Conservation  of  mass  then 
requires  that  the  average  density  compression  ratio  (essentially  the  solid 
fraction  only)  is  given  by 
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With  x  -  55  ft  and  R  —  130  ft  for  Rainier,  the  result  r|  =  1 .  08  indicates  the  small 
void  ratio  fraction  of  0.  08  in  comparison  with  void  fractions  like  20  or  30%  for 
dried- out  samples. 

III.  Rainier  Re  stilts 


The  Rainier  explosion  is  characterized  by  good  containment,  with  essen¬ 
tially  spherical  symmetry,  in  which  the  waste  heat  concept  provides  the  con¬ 
trolling  mechanism. 

The  idealized  equation  of  state  for  Nevada  soil  and  the  observed  shock  ve¬ 
locity  versus  distance  curve  for  Rainier  are  shown  qualitatively  in  Figs.  3.4a 
and  3.4b  of  Part  B,  which  also  describes  the  interrelations  of  various  regions 
on  those  two  figures.  The  curves  shown  here  represent  the  author's  views  at 
Armour  Research  Foundation  for  what  should  have  happened  and  what  did 
happen  on  Rainier,  but  the  reader  is  cautioned  that  this  interpretation,  includ¬ 
ing  the  "decrepitation"  concept,  do  not  yet  necessarily  represent  the  views  of 
the  sponsor. 

Doppler  measurements  fixed  the  shock  velocity  in  soil  as  70  meters/msec 
at  the  walls  of  the  nuclear  device  chamber,  an  effective  radius  like  1.  6  meters. 
The  switches  showed  that  the  shock  velocity  fell  quickly  to  a  value  of  3.  2 
meters/msec  at  4.  9  meters,  and  ran  constant  at  this  value  out  to  18.4  meters. 
This  region  corresponds  to  a  peak  pressure  of  40,  000  bars,  and  the  flat  por¬ 
tion  of  the  curve  represents  the  transition  from  an  explosion  in  a  gas  (made 
by  vaporization  of  rock  for  shock  radii  inside  4.  9  meters)  to  the  characteris¬ 
tics  typical  of  an  explosion  in  soil  (beyond  18.4  meters).  Ambient  sound  ve¬ 
locity  in  Nevada  is  about  8,  000  ft/sec;  shock  velocity  was  measured  as  low  as 
7,  400.  ft/sec,  and  was  still  falling  rapidly.  As  expected  and  as  the  gauges 
were  so  designed,  the  switches  beyond  the  crushed  zone  limit  at  130  feet  did 
not  close. 

It  is  emphasized  that  the  waste  heat  so  strongly  controls  the  phenomenol- 
ogy  that  a  nuclear  explosion  in  soil  is  best  regarded  as  a  calorimeter  experi¬ 
ment:  virtually  all  the  heat  was  expected  to  be  distributed  within  the  cavity 
and  the  relatively  thin  shell  of  crushed  material  out  to  the  crushing  radius. 

One  equation  for  yield  —  used  for  Rainier  —  is  instructive  on  this  point.  The 
energy  balance  in  Eq.  (1)  of  Part  B  is  readily  put  in  the  calorimetric  form 
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4  3  •  • 

The  term  ir  R_,  is  the  volume  out  to  the  point  where  any  vaporization  occurs 

in  soil,  which  gives  R2  =  18.4  meters  for  Rainier;  p0  Q2  is  the  energy  density 
deposited  by  waste  heat  up  to  the  boiling  or  decrepitation  point.  The  term 
(Q./Q-,  -  1)(R,R:?)  is  a  correction  factor  to  account  for  waste  heat  in  excess 
of  simply  raising  to  the  boiling  point;  on  Rainier  this  term  has  the  numerical 
value  0.  07  in  Eq.  (2).  The  term  Pf/p0Q2  represents  the  available  hydrodynamic 
energy  when  the  shock  reached  18.4  meters  and  had  the  value  0.  23  for  Rainier. 
The  fraction  l/l+O.  07+0.  23  =0.  77  means  that  77%  of  the  energy  was  dissipated 
within  18.4  meters  or  60. feet  from  the  device  merely  from  raising  the  rock  to 
the  boiling  point  of  water.  The  available  hydrodynamic  energy  already  repre¬ 
sented  only  a  relatively  small  fraction  of  total  energy  by  then,  and  the  total 
yield  is  not  sensitive  to  an  accurate  specification  of  hydrodynamic  waveforms. 

With  a  crushing  pressure  of  700  bars  (an  upper  estimate  in  UCRL-5124), 
Eq.  (3)  of  Part  B  indicates  that  a  maximum  for  the  available  energy  left  at  the 
cessation  of  the  crushing  stage  was 


W' 

W 


=  0.23 


700  y/2 

40, 000 J 


0.03  . 


In  short,  only  3%  of  the  energy  was  left,  even  at  130  feet,  to  go  into  the  seis¬ 
mic  wave;  97%  of  the  energy  had  already  been  converted  to  heat. 

The  analytic  solution  provides  a  similar  evaluation  of  the  yield  at  any 
other  radius  and  time  by  more  elaborate  counterparts  of  Eq.  (2)  above.  The 
constancy  of  W  =W(t)  so  determined  is  an  internal  check  on  the  validity  of  the 
theory,  provided  the  data  and  the  shot  geometry  are  sufficient  to  insure  rea¬ 
sonable,  spherical  symmetry.  The  probable  error  of  the  Rainier  hydrodynam¬ 
ic  yield,  about  1.  7  +  0.  15  kt,  is  largely  due  to  the  uncertainty  from  the  neces¬ 
sity  to  "construct"  the  equation  of  state  for  Nevada  soil. 

The  principal  damage  to  the  Rainier  tunnel,  beyond  the  crushed  zone,  was 
due  to  the  "pinch  jet"  similar  to  that  described  in  Section  VII  of  Part  B.  The 
temperatures  in  such  an  induced  jet  are  significantly  high.  With  a  shock  ve¬ 
locity  in  soil  of  7900  ft/sec,  the  model  of  Fig.  5.  6b.,  Part  B,  would  require 
air  shock  pressures -of  the  order  of  100  bars,  i.  e.  ,  "fireball"  temperatures. 
Such  temperatures  easily  account  for  the  observed  carbon  monoxide  in  the 
Rainier  tunnel.  The  carbon  monoxide  may  have  been  produced  by  air  shock 
action  on  the  mine  timbers. 
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IV.  Hardtack  II  Results 

The  Hardtack  II  events  were  characterized  by  the  channeling  concept,  to 
about  the  same  degree  that  Rainier  was  characterized  by  the  waste  heat  con¬ 
cept. 

The  channeling  concept  and  its  implications  are  described  in  Sections  VII 
and  VIII  of  Part  B.  The  shot  geometry  for  each  of  the  Hardtack  II  events  in¬ 
volved  line-of- sight  pipes  of  relatively  large  diameters,  such  as  1  or  2  feet. 
Also,  the  accelerated  time  schedule  of  Hardtack  II,  due  to  the  impending  mor¬ 
atorium  on  weapons  tests  on  1  November  1958,  required  the  positioning  of 
blast  switches  in  the  floor  of  tunnels  for  some  shots,  precisely  where  the 
asymmetries  due  to  jetting  were  most  pronounced.  As  such,  the  time-of- 
arrival  data  are  of  excellent  value  as  diagnostic  hydrodynamics,  whereas  it 
would  be  naive  to  expect  a  reliable  yield  determination  with  such  experimental 
improvisations  and  shot  geometry.  Here  again,  the  reader  is  cautioned  that 
the  results  given  here  represent  the  interpretation  of  the  author  at  ARF  for 
what  should  have  happened  and  what  did  happen  on  Hardtack  II,  but  this  pre¬ 
sentation  does  not  necessarily  imply  the  technical  concurrence  of  the  sponsor. 

Tamalpais  was  instrumented  like  Rainier  in  rock,  but  not  along  the  axis 
of  any  conjectured  jets.  The  results  confirmed  the  existence  of  subsonic  ve¬ 
locities  in  the  surrounding  rock,  even  close  to  the  device  chamber,  presum¬ 
ably  due  to  the  large  fraction  of  energy  which  jetted  elsewhere. 

The  results  for  Logan  are  shown  qualitatively  in  Fig.  3.  2.  The  line  =  of- 
sight  pipe  shown  was  evacuated,  had  a  maximum  diameter  of  24  inches,  and 
passed  through  three  intermittent  zones  of  sandbagging  in  the  tunnel.  The 
blast  switches  were  placed  on  the  floor  of  the  tunnel,  in  both  the  open  tunnel 
and  the  sandbagged  sections.  The  measurements  generally  showed  an  abrupt 
increase  in  shock  velocity  at  each  of  the  open  spaces,  and  a  much  slower  ve¬ 
locity  within  each  sandbagged  layer.  To  Armour  Research  Foundation,  this 
result  appears  to  be  a  vivid  demonstration  of  jetting  down  the  pipe,  flooding 
of  energy  into  the  open  tunnel  spaces,  partial  closure  by  sandbags  followed  by 
a  new  jet  through  the  sandbag  zone.  The  jet  velocity  showed  a  marked  de¬ 
crease  in  the  third  open  space,  as  if  the  third  set  of  sandbags  held,  but  one 
switch  at  187.  5  feet,  in  this  open  section  of  the  tunnel,  closed  earlier  than  a 
switch  in  the  sandbags  at  183.  5  feet,  indicating  that  the  shock  energy  was  flow- 
ing  backward  through  the  sandbags  at  this  interface.  It  is  not  clear  from  the 
data  +hat  even  the  last  set  of  sandbags  (beyond  the  line-of- sight  pipe)  actually 
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Fig.  3.  2.  Logan  —  time  of  arrival. 
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held  the  jet  at  the  far  end  of  the  tunnel;  the  last  switches  indicate  the  shock 
was  possibly  accelerating  beyond  that  point.  The  Logan  results  indicate  an 
apparent  energy  release  as  low  as  5  kt  where  the  explosion  was  partially 
stemmed  at  each  sandbag  barricade,  but  an  apparent  yield  like  30  kt  where 
jetting  opened  the  tunnel.  In  the  opinion  of  ARF,  the  jetting  effect  in  the 
Logan  tunnel  is  simply  the  explanation  for  the  damage  which  occurred  in  the 
main  tunnel. 

Evans  event  again  confirmed  the  subsonic  velocities  in  soil,  due  in  part 
(again  in  the  author  ?s  opinion)  to  jetting.  The  Evans  shot  involved  large  air- 
filled  line-of- sight  channels,  the  first  being  about  130  feet  long  which  opened 
into  a  large  room.  Significantly,  a  blast  switch  in  that  room,  placed  well  off 
the  line~of=  sight,  closed  at  0.  7  msec,  at  a  time  when  the  shock  had  advanced 
only  about  4  feet  into  the  rock  surrounding  the  device  chamber.  The  ratio 
130/4  is  in  excellent  agreement  with  the  ratio  of  30;  1  (air  shock  velocity  to 
ground  shock  velocity)  suggested  in  Part  B  for  jets  at  high  pressures. 

F or  Blanca,  the  tunnel  damage  from  Logan  precluded  the  use  of  the  orig¬ 
inal  nuclear  device  location  and  switch  layouts,  but  attempts  were  made  to 
improvise  a  layout  as  on  Logan.  The  close-in  string  of  switches  failed  for 
one  of  many  possible  but  unknown  reasons.  The  last  string  indicates  that  the 
last  set  of  sandbags  successfully  stemmed  the  explosion  at  that  point.  The 
overall  results  for  Blanca  are  sufficient  to  confirm  the  expectation  that  about 
10  times  as  much  energy  from  the  nuclear  device  was  available  for  crushing 
as  would  have  resulted  from  a  stemmed  charge,  and  a  correspondingly  larger 
fraction  of  energy  went  into  the  seismic  wave. 

V.  Summary 

In  summary  of  these  Armour  Research  Foundation  studies,  it  appears 
entirely  feasible  to  measure  the  hydrodynamic  yield  by  time-of-arrival  methods, 
provided  the  explosion  is  either  stemmed  to  permit  spherical  symmetry,  or, 
in  the  presence  of  open  line-of- sight  channels,  if  sufficient  strings  of  switches 
are  placed  along  different  radii  to  verify  the  jetting  and  measure  its  extent.  A 
principal  advantage  of  the  switch  method  is  the  simplicity  and  reliability  with 
which  many  hundreds  of  pertinent  points  of  information  can  be  obtained  at  pres¬ 
sure  levels  far  in  excess  of  the  capabilities  of  us'ual  pressure-instrumentation. 

Some  consequences  of  the  waste  heat  concept  on  wave  propagation  are  listed 
briefly  in  Section  II  of  Part  B,  and  their  implications  for  underground  explosions 
in  Section  III.  Some  implications  of  the  channeling  concept  (Section  VII)  are 
discussed  in  Section  VIII,  "Controlled  Excavation." 
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Part  B 


New  Concepts  for  Control  and  Use  of  Nuclear  Explosions 

This  section  describes  some  new  and  basic  concepts  for  control  and  use  of 
nuclear  explosions  in  peaceful  uses.  These  applications  involve  large  scale 
public  works  for  which  the  designation  "civil  engineering"  is  particularly  ap¬ 
propriate. 

Because  some  of  the  ideas  presented  here  may  seem  so  novel  as  to  oppose 
first  intuition  about  them,  it  probably  helps  to  state  at  the  outset  that  the  field 
"Fluid  Dynamics  of  Strong  Shocks"  is  literally  a  new  branch  of  physics.  His¬ 
torically,  this  field  grew  out  of  the  military  necessity  to  understand  the  blast 
and  other  hydrodynamic  effects  of  nuclear  weapons;  it  is  only  fair  to  add  that 
the  progress  would  not  have  been  possible  without  the  enormous  experience 
invested  in  that  field  during  the  past  fifteen  years.  The  military  applications 
clearly  could  not  be  paralleled  by  a  similar  growth  in  the  nonclassified  litera¬ 
ture.  But  the  basic  concepts  are  not  classified,  they  are  readily  derived  from 
conservation  of  mass,  momentum,  and  energy.  The  field  is  characterized  by 
the  nonlinearity  in  the  hydrodynamic  equations  of  motion  and  in  the  equation 
of  state  of  materials,  due  to  the  high  pressures  and  stresses  characteristic 
of  such  explosions.  The  end  result  is  an  emphasis  on  different  mathematical 
disciples  and  theoretical  techniques,  new  but  no  more  difficult  than  their  count¬ 
erparts  in  "acoustic"  theory.  Unfortunately,  the  linearity  assumption  pervades 
so  many  fields  of  science  and  engineering  —  aerodynamics,  elasticity,  vibra¬ 
tion*  electronics  —  that  we  develop  a  corresponding  intuition  about  the  behav¬ 
ior  of  natural  processes  as  if  linearity  were  a  universal  characteristic  of 
nature.  But  in  the  realm  of  very  large  displacements  an4  unbalance  of  forces 
the  concepts  of  linearity  are  more  of  a  liability  than  an  asset.  In  some  cases, 
the  linearity  assumption  literally  annihilates  the  most  significant  terms  in  the 
equations  of  motion  and  denies  the  essence  of  the  problem  we  are  trying  to 
understand.  A  corresponding  confusion  would  result  if  we  were  trying  to  under¬ 
stand  a  radio  or  television  set  which  was  built  with  nonlinear  resistances  so  that 
electrical  resistance  was  not  a  constant  property  of  material,  independent  of 
voltage  and  current,  but  proportional  to  some  high  power  of  the  current  itself. 

This  section  is  restricted  to  a  discussion  of  two  principal  concepts  and 
their  applications:  waste  heat  and  channeling.  The  waste  heat  concept  pro¬ 
vides  a  categorical  method  for  the  control  and  use  of  an  explosion  by  control 
of  the  properties  of  the  medium  surrounding  the  device  or  energy  source;  the 
applications  are  possible  because  of  nonlinearity  of  the  equation  of  state  for 
compression.  The  channeling  concept  provides  a  categorical  method  for  con¬ 
trol  and  use  of  an  explosion  by  control  of  the  geometry  of  the  material  surround¬ 
ing  the  energy  source;  the  applications  are  po s sibl e  mo stly  because  of  the  non¬ 
linearity  in  the  shock  velocity  and  in  the  material  velocity  behind  the  shock. 

I.  Waste  Heat 

Irreversible  heating  in  a  material,  caused  by  the  entropy  change  across 
a  shock  which  traverses  it,  is  well  known.  The  effect  is  proportional  to  the 
cube  of  the  overpressure  at  low  pressures  and  is,  therefore,  vanishingly  small 
in  the  acoustic  approximation.  There  is  a  tendency  to  overlook  the  correspond¬ 
ing  fact  then  that  the  entropy  change  is  necessarily  a  controlling  feature  of  the 
hydrodynamics  at  the  high  pressures  which  exist  close  to  any  rapid  energy  re~ 
lease.  Thus,  the  fireball  from  a  nuclear  device  in  air  is  not  due  to  the  device 
parts  directly  but  due  to  the  passage  of  the  shock  wave  over  the  surrounding 
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air;  if  the  air  were  adiabatic  in  compression  (let  alone  its  nonlinearity)  there 
would  be  no  fireball  from  a  nuclear  device.  Also,  the  "hole"  represented  by 
the  bubble  of  an  underwater  explosion  and  the  cavity  in  an  underground  explo¬ 
sion  would  exist  for  other  hydrodynamic  reasons,  but  they  would  be  different 
from  their  present  size* 

We  define  the  waste  heat  Q(P)  to  be  the  net  energy  per  unit  mass  (above 
initial  energy  content)  which  remains  in  the  material  (originally  shocked  to 
pressure  P)  after  it  has  finally  expanded  adiabatically  down  to  the  ambient 
pressure  P0  that  existed  before  the  explosion.  The  nwaste  heat  concept"  it¬ 
self  has  a  special  meaining  with  regard  to  explosions  in  solids,  liquids,  and 
their  mixtures  with  gases;  it  is  a  convenient,  useful,  and  accurate  engineering 
assumption  which  justifies  the  subsequent  bookkeeping  of  internal  energy  by 
dividing  it  into  "waste  heat"  and  "available  energy."  The  waste  heat  concept 
implies  important  but  deceptively  simple  differences  from  entropy  change: 
there  would  be  no  waste  heat  if  the  material  returned  to  zero  pressure  instead 
of  ambient  pressure.  This  bookkeeping  approximation  is  justified  by  the  fact 
that  the  bulk  modulus  for  compressibility  is  quite  insensitive  to  temperature 
and  for  a  variety  of  other  compensating  effects  both  in  the  process  of  shock 
and  in  the  subsequent  expansion. 

The  relative  magnitude  of  several  energy  quantities  and  the  method  of 
calculating  waste  heat  is  indicated  in  Fig.  3.  3.  The  figure  is  a  conventional 
pressure-volume  diagram;  energy  is  portrayed  as  an  area  on  such  a  diagram. 
The  full  curve  is  called  the  "Rankine-Hugoniot  curve"  which  gives  the  relation 
between  pressure  P  and  specific  volume  V  at  various  shock  strengths.  By 
conservation  of  energy  alone,  the  total  energy  added  to  the  material  by  the 
shock  is  given  by  E»p  -  E0  =  P(V0  -  V).  This  is  the  area  of  the  rectangle  be¬ 
low  the  pressure  P  and  between  V0  and  V.  The  kinetic  energy  E^  of  the  ma¬ 
terial,  with  velocity  u  behind  the  shock,  is  given  by  conservation  of  mass  and 
momentum  alone  as:  1/2  u^  =  1/2  (P  -  Po  )(V0  -  V).  This  is  the  triangle  in  the 
upper  right  of  the  total  energy  rectangle:  whatever  the  actual  path  during  shock 
compression,  the  latter  equation  shows  that  the  compres sion  is  effectively  as 
if  the  material  were  shocked  along  the  straight  line  joining  P0  with  P.  It  im¬ 
mediately  follows  that  the  internal  energy  in  the  material  is  the  quantity  of 
E^  -  E0  =  1/2  (P  +  P0  )(Y0  -  V)  which  is  clearly  the  total  trapezoid  in  the  total 
energy  rectangle  lying  below  the  straight  line  joining  P  and  P0  .  Although  these 
energy  expressions  and  the  Rankine-Hugoniot  relations  may  be  derived  in  a 
different  manner  with  more  restrictive  assumptions,  the  simple  derivation 
shown  here  is  sufficient  to  indicate  that  these  relations  are  quite  general  and 
apply  without  regard  to  the  equation  of  state  of  the  material  or  the  details  of 
the  shock  process. 

For  solids  and  liquids  the  subsequent  expansion  of  the  material  behind  the 
shock  can  be  described  in  a  simple  manner:  the  material  will  expand  by  essen¬ 
tially  adiabatic  processes  back  to  its  original  pressure  PQ  and  volume  V0  ,  and 
for  most  practical  purposes,  will  expand  along  the  Rankine-Hugoniot  curve. 
Since  the  area  under  that  curve  is  the  fP  dV  energy  subsequently  given  up 
during  the  expansion,  that  area  is  defined  as  the  available  internal  energy.  It 
follows  that  the  lens- shaped  area  lying  between  the  straight  line  and  the  curve 
is  by  definition  the  waste  heat  itself.  The  energy  represented  by  this  area 
appears  as  a  final  temperature  rise  in  the  material  left  long  after  the  shock 
has  passed.  Heat  conduction  processes  occur  on  time  scales  much  longer 
than  the  rapid  process  of  compression  and  expansion  described  here,  so  the 
deposition  of  waste  heat  is  literally  forever  so  far  as  the  shock  hydrodynamics 
is  concerned. 
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For  all  material s,  there  exists  some  shock  pressure  high  enough  that  the 
waste  heat  is  sufficient  to  vaporize  the  material  irreversibly,  even  after  it  has 
expanded  and  coded  back  to  ambient  pressure.  We  define  that  shock  pressure 
to  be  the  1  'vaporization M  pressure,  and  the  material  inside  the  corresponding 
radius  becomes  part  of  the  device  cavity.  The  energy  represented  by  the 
’’waste  heat  area"  in  excess  of  the  total  heat  to  vaporization  appears  as  super¬ 
heat  in  a  gas,  is  available  by  expansion,  and  is  not  counted  as  waste  heat. 

Very  close  to  the  device,  well  in  excess  of  vaporization  pressure,  all  mate¬ 
rials  are  initially  good  transmitters. 

It  is  easy  to  find  or  to  fabricate  materials  and  mixtures  which  "crush 
down"  at  some  relatively  low  pressure.  In  that  case  the  Rankine-Hugoniot 
curve  is  much  more  sharply  bowed,  as  shown  in  the  diagram  at  the  upper 
right  of  Fig.  3.  3,  and  nearly  all  the  internal  energy  appears  as  waste  heat. 

As  a  result  of  the  Rankine-Hugoniot  relations,  all  materials  compress  to 
some  limiting  density  ratio,  like  10  to  1  or  4  to  1  even  at  infinite  shock  pres- 
sureso  For  this  reason  all  materials  display  a  vertical  Rankine-Hugoniot 
curve  at  pressures  high  enough  and  all  materials  become  ideal  absorbers  at 
such  high  pressures  provided  they  are  not  greatly  superheated. 

At  low  pressures  the  curvature  in  the  Rankine-Hugoniot  curve  is  very 
slight;  in  the  acoustic  approximation  the  Rankine-Hugoniot  curve  is  a  straight 
line  passing  through  (P0  ,  V0).  This  is  the  domain  of  ordinary  elastic  theory. 
Strictly  speaking,  the  straight  line  approximation  is  never  really  reached,  but 
the  waste  heat  area  becomes  very  small  and  we  speak  of  material  as  an  ideal 
transmitter.  All  materials,  at  least  below  their  crushing  strength,  are  ideal 
transmitters  at  low  pressures,  as  shown  at  the  lower  right  of  Fig.  3.  3. 

II.  Wave  Propagation 

We  discuss  briefly  the  fraction  of  energy  which  is  hydrodynamically  avail¬ 
able  after  the  shock  has  run  to  some  appreciable  distance  R  from  an  instan¬ 
taneous  energy  release  in  a  liquid  or  a  solid.  The  damage,  and  hence  the  dif¬ 
ficulty  of  controlling  an  explosion,  can  often  be  specified  in  terms  of  the  peak 
pressure  P  at  the  shock  front. 

At  high  shock  pressures,  the  usual  theory  in  the  current  literature  is 
based  on  certain  similarity  as  sump  cions,  which  are  not  really  justified  in  their 
own  right,  but  the  similarity  condition  would  require  that  the  decay  of  pressure 
with  distance  be  as  P  ~R“3,  the  quantity  PR3~  W,  and  (d  In  P)/(d  In  R)  =  -3. 

At  low  pressures  the  usual  description  is  based,  on  the  acoustic  approximation 
and  gives  (P  -  PQ )  R  ~W  and  [d  in  (P  -  PQ)]  /(d  In  R)  =  -l.  In  both  cases,  the 
total  energy  W  in  the  wave  is  taken  to  be  constant  and  equal  to  the  original 
energy  released.  These  results,  both  at  high  and  low  pressures,  are  mislead¬ 
ing  oversimplifications  and  the  reasons  can  be  summarized  briefly  by  means 
of  the  waste  heat  concept. 

With  the  energy  partitioned  according  to  the  waste  heat  concept,  the  en¬ 
ergy  balance  at  any  time  is: 


Effective 

energy 


W» 


Original  _  Dissipated  __  Available  internal  plus 

energy  energy  kinetic  energy 

R 

w  -  4tt  y  p0QR2dR  =  ^  IT  R3(P  -  Po)  f. 
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The  validity  of  this  equation  is  wholly  a  matter  of  definition  and  conservation 
of  energy.  The  effective  energy  is  W*,  the  energy  left  to  propagate  the  shock 
and  to  do  damage  at  the  distance  R.  The  extreme  right-hand  term  contains 
the  volume  of  the  shock  (4/3)ttR^  and  the  peak  pressure  P  at  the  shock  front, 
pressure  having  dimensions  of  energy  per  unit  volume.  The  factor  f  is  de¬ 
fined  to  relate  the  average  energy  density  per  unit  volume  on  the  interior  of 
the  wave  (available  internal  energy  plus  kinetic  energy)  to  the  peak  pressure 
at  the  shock  front.  Because  the  sum  of  compressional  energy  and  kinetic  en¬ 
ergy  decreases  monotonically  on  the  interior  of  the  wave,  the  factor  f  will 
be  less  than  unity,  in  terms  of  the  peak  energy  density  at  the  shock  front  it¬ 
self;  f  will  be  a  number  like  1/2  for  strong  shocks  in  gases  but  be  a  much 
smaller  number  when  the  shock  becomes  weak.  The  waste  heat  concept  makes 
it  possible  to  calculate  the  dissipated  energy  by  integrating  over  the  past  his¬ 
tory  of  the  shock  front,  since  the  functions  P(R)  and  Q[P(R)]  are  known  and  the 
heat  is  ''stored  forever"  in  the  material.  Because  the  waste  heat  Q  is  always 
a  positive  quantity  the  integral  is  monotonically  increasing.  Equation  (1)  shows 
that  the  effective  energy  W1  (or  its  counterpart  on  the  extreme  right)  is  not 
proportional  to  a  constant  W  but  decreases  continually  as  thd  shock  progresses 
and  the  waste  heat  accumulates  in  the  material  left  behind. 

Assuming  that  the  original  energy  W  is  constant  (it  is  not  necessary  to 
assume  that  the  energy  release  was  instantaneous),  Eq.  (1)  may  be  differen¬ 
tiated  with  respect  to  time  or  distance  R  and  after  suitable  manipulations  the 
rigorous  result  is  obtained  that 


At  high  pressures,  f  =  constant  =  l/2  or  1/3,  (d  In  f)/[d  In  (?  -  P0)]  =  0  and 
at  the  same  time  the  quantity  p0  Q/(P  -  P0)  is  itself  close  to  1/2.  The  result 
is  that  the  logarithmic  decay  of  pressure  with  distance  is  two  or  three  times 
greater  than  the  similarity  assumption  would  indicate  and  means  that  the  peak 
pressure  decay  is  like  the  sixth  or  perhaps  the  ninth  power  with  distance  from 
the  source.  This  rapid  decay  of  pressure,  says  P  ~R~  ,  with  distance  instead 
of  P  ~R“^  is  due  to  dissipation  of  energy  through  the  waste  heat  process;  it 
results  in  dissipation  of  energy  fractions  like  99%  of  the  original  energy  W 
before  the  shock  wave  reaches  the  crushing  pressure  or  other  pressures  of 
practical  interest. 

It  can  be  shown  that  at  low  pressure^  the  factor  f  is  proportional  to  the 
square  of  the  overpressure,  f  ~(P  -  P0 )  ,  so  that  (d  In  f)/[d  In  (P  -  P0 )]  -  2 . 

It  follows  from  Eq.  (2)  that  [  d  In  (P  -  P0  )5  / (d  h  R)  =»  {l  +  [  p0  Q/(P  -  P0‘)f]  }  .  3 

But  in  the  acoustic  approximation  Q  ~(P  -  P0)^  and  because  (P  -  P0)f  ~(P  -  P0  ) 
it  follows  that  the  correction  term  p0Q/(P  -  P0  )f  (representing  the  waste  heat) 
does  not  vanish  but  approaches  a  constant  limiting  value  even  in  the  acoustic 
approximation. 

Despite  the  variations  in  the  factor  f  with  shock  strength  and  the  variation 
in  the  behavior  of  the  waste  heat  with  pressure,  it  is  possible  to  specify  the 
available  energy  W1  by  a  simple  relationship  with  pressure  as: 
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Pi  is  the  1  Vaporization  pressure11;  the  term  Wi  is  the  available  energy  at  the 
vaporization  pressure;  the  W*  will  usually  be  about  l/2  the  original  energy 
release  W.  The  pressure  Pj_  will  be  of  the  order  of  tens  of  kilobars  for  most 
soils  but  may  be  of  the  order  of  megabars  for  dense  homogeneous  materials 
like  concrete  or  steel  if  free  of  voids.  In  any  case,  Eq.  (3)  shows  why  the 
available  energy  W*  will  be  only  a  small  fraction  of  the  original  energy  re¬ 
lease  W,  and  will  decrease  with  distance. 

Some  consequence  of  the  waste  heat  concept  on  wave  propagation  follow 
from  such  considerations  as  these. 


1.  Most  of  the  energy  from  the  explosion,  fractions  like  99%,  is  deposited 
as  heat  very  close  to  the  source,  in  the  cavity  and  in  a  relatively  thin  shell  of 
compressed  material  immediately  surrounding  the  nuclear  device  cavity. 

2.  The  energy  available  to  produce  damage  or  seismic  waves  bears  no 
set  direct  relation  with  the  original  energy  release;  the  residual  energy  frac¬ 
tion  depends  upon  the  initial  configuration  of  the  material  and  the  room  im¬ 
mediately  surrounding  the  device,  upon  the  vaporization  pressure,  and  the 
dynamic  crushing  strength  of  the  rock  or  soil.  The  apparent  energy  at  long 
distances  will  be  a  small  fraction  of  the  original  energy  release,  like  0.  001. 

3.  The  concept  of  energy  partition,  so  much  going  into  the  bubble,  so 
much  going  into  blast,  is  essentially  a  misnomer;  the  available  energy  is  a 
constantly  decreasing  fraction  with  distance  and  the  accumulated  waste  energy 
rapidly  approaches  100%  of  the  original  energy. 

4.  It  is  doubtful  if  the  acoustic  approximation  ever  applies,  in  the  sense 
of  a  constant  energy  in  the  wave.  For  example,  Cole’s  "Underwater  Explo¬ 
sions"  gives  the  universal  experimental  for  TNT  explosions  in  water  in  the 


form 


(P  -  P0)  =  constant 


(4) 


From  the  exponent  1.13  and  Eq.  (2)  we  recognize  that  whatever  the  value  of  Q 
and  f  at  low  pressures,  the  limiting  value  is  such  that  Q/P  -  P0)f  =0.  13. 

5.  If  the  attempt  is  made  to  measure  the  energy  release  from  an  explo¬ 
sion  by  measuring  the  kinetic  and  compressional  energy  in  the  shock  at  ap¬ 
preciable  distances  from  the  source,  the  results  are  bound  to  indicate  very 
small  energy  releases  perhaps  in  error  by  a  factor  like  100  from  the  original 
energy  release.  Such  procedures  have  been  used  by  a  number  of  agencies, 
e.  g.  ,  in  attempting  to  evaluate  the  energy  release  from  metal-water  reactions 
suspected  to  occur  in  reactor  explosions;  the  results  indicated  very  small  energy 
releases.  The  point  is  that  a  similar  measurement  from  a  TNT  explosion  would 
also  lead  to  similarly  small  values  of  hydrodynamic  energy,  at  large  distances, 
entirely  inconsistent  with  the  known  energy  release  of  TNT,  as  used  for  scaling. 
Damage  is  usually  predicted  on  the  basis  of  scaling  from  TNT  results,  so 

that  direct  comparisons  are  valid,  and  this  "apparent  energy"  grossly  under¬ 
estimates  the  damage. 

Ill .  Underground  Explosions 

A  typic ah  history  from  an  explosion  in  soil  displays  striking  features,  in¬ 
dicated  in  Figs.  3.4a  and  3.4b,  which  would  not  be  suspected  on  the  basis  of 
either  linearized  or  conventional  shock  theory.  This  description  is  the  author’s 
interpretation  for  what  should  have  happened  and  what  did  happen  in  the  Rainier 
underground  shot  at  Nevada  Test  Site  in  September  1957;  the  reader  is  re¬ 
minded  that  the  interpretation  is  not  yet  universally  accepted. 
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Fig.  3.4a.  ,  Pressure- velocity  diagram  for  soil. 
Fig,  3.  4b.  Shock  velocity  versus  distance  in  soil. 


The  plot  of  shock  velocity  versus  distance  as  in  Fig,  3,4b  is  convenient  to 
describe  the  various  phases,  in  the  history  of  the  shock  growth.  At  high  enough 
pressures,  conservation  of  mass  and  momentum  alone  requires  the  shock  ve¬ 
locity  to  behave  as  U  ~pl/2  for  any  medium.  The  similarity  condition 
P  ~R~3  would  require  d  In  U/d  In  P  =  -3/2,  which  is  a  straight  line  of  slope 
-3/2  on  these  logarithmic  coordinates.  At  low  pressures  the  shock  velocity 
asymptotically  approaches  a  constant,  which  is  sound  velocity;  the  general 
behavior  according  to  conventional  theory  is  shown  as  the  dashed  line  in  Fig. 
3.4b.  For  the  purpose  of  exposition  now,  we  assume  the  Rankine-Hugoniot 
curve  has  the  shape  shown  in  the  pressure-volume  diagram  of  Fig.  3.4a.  The 
corresponding  regions  of  growth  in  Figs.  3,  4a  and  3.  4b  are  shown  by  labeling 
portions  of  the  curve  A,  B,  C,  D,  and  E.  The  region  A  is  shocked  well  above 
the  vaporization  pressure  of  the  material.  The  material  behaves  like  a  gas 
because  the  material  behind  it  is  irreversibly  vaporized  and  superheated.  But 
the  shock  velocity  curve  decays  more  rapidly  with  distance,  than  U  ~R~  / 

(like  U  ~pl/2),  because  of  the  enormous  energy  expended  on  latent  heat  of 
vaporization  of  the  material.  The  region  B  corresponds  to  the  range  of  pres¬ 
sures  between  the  pressure  at  which  100%‘of  the  material  is  irreversibly  va¬ 
porized  and  some  lower  pressure  at  which  all  of  the  material  can  be  heated 
to  the  boiling  point  by  waste  heat  but  too  weak  to  vaporize  any  of  it.  B  is  a 
region  of  more  or  less  constant  temperatures.  For  other  reasons  also,  the 
peak  pressure  is  probably  constant,  and  in  any  case,  the  shock  velocity  runs 
constant  for  some  distance.  In  Rainier,  this  .constant  velocity  region  extended 
for  all  material  which  originally  lay  between  5  meters  out  to  17  or  18  meters 
from  the  device.  The  regions  C  and  D  shown  are  best  understood  by  noting 
that  the  Rankine-Hugoniot  equation  for  shock  velocity  requires0. 


<P  -  P0) 
Vo'  -T 


(5) 
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This  equation  means  that  the  straight  line  joining  P  and  PQ  ,  as  shown  in  Fig. 

3C  3,  represents  the  square  of  the  shock  velocity,  U^.  Ambient  sound  velocity 
has  been  indicated  by  the  dashed  line  on  Fig,  3.  4a,  as  an  extension  of  a  straight 
line  portion  of  the  Rankine-Hugoniot  curve  at  low  pressures  in  the  region  E0 
The  region  C  corresponds  to  those  pressures  below  vaporization  but  at  which 
the  shock  is  still  sufficiently  strong  to  be  supersonic  with  respect  to  ambient 
sound  velocity.  The  region  D  is  particularly  interesting  because  it  represents 
a  region  of  pressures  well  above  the  crushing  pressure  but  in  whicK\the  shock 
velocity  is  below  ambient  sound  velocity.  In  this  region  of  pressures,  precur¬ 
sor  sound  signals  may  precede  the  main  shock  pulse,  but  they  are  necessarily 
insignificant  in  strength  compared  with  the  main  shock  pressure.  The  region 
E  represents  the  range  of  pressures  below  the  crushing  strength  of  the  rock. 

By  virtue  of  the  simplification  used  to  portray  compressibility  in  the  region 
E  in  Fig.  3.4b  and  in  the  corresponding  region  E  in  Fig0  3.  4a,  the  crushing 
pressure  is  shown  as  a  horizontal  straight  line. 

We  note  that  acoustic  theory  would  require  that  the  entire  history  of  the 
explosion  would  be  represented  by  straight  lines,  both  for  the  compressibility 
in  Fig.  30  4a  (dashed  line)  and  a  constant  value  (sound  velocity)  represented  by 
E  in  Fig.  3.  4b0 

Some  comments  on  the  history  of  an  underground  explosion  are  given 
below. 

1.  There  is  little  meaning  to  a  ^constant  seismic  energy^  in  the  acoustic 
region  Es  which  is  by  definition  beyond  the  crushed  zone,  a  distance  like  100 
feet  for  one  kiloton  in  a  medium  like  Nevada  soil.  Surprisingly  enough, 

Nevada  rock  is  a  relatively  good  transmitter  of  shock  energy  because  of  its 
high  water  content  at  depth.  Energy  fractions  like  1%  of  the  original  energy 
may  still  be  available  in  the  seismic  wave  at  this  distance,  but  this  would  be 
significant  only  if  the  material  beyond  the  crushed  zone  were  completely  homo¬ 
geneous  and  linearly  elastic. 

2.  The  existence  of  faults,  joints,  fissures,  and  free- surfaces  all  con¬ 
stitute  mechanisms  whereby  the  compression  (due  to  the  shock  transmitted 
beyond  the  crushed  zone)  is  readily  relieved  (at  the  expense  of  the  net  energy 
transmitted  further  out  in  a  seismic  wave).  As  such,  the  !!earthquakes n  gen¬ 
erated  by  these  explosions  are  more  or  less  coincidental  to  the  shock  process 
and  their  relationship  to  the  energy  release  is  ambiguous.  Both  the  Rainier 
event  (1957)  and  the  Hardtack  II  events  (1958)  resulted  in  displacements  of 
about  two  feet  along  joints  or  planes  of  structural  weakness  originally  in  the 
rock,  much  like  the  slippage  along  an  earthquake  fault.  Very  large  amounts 
of  energy  may  be  represented  as  potential  energy  from  a  fault  block  lying 
above  some  plane  of  structural  weakness,  and  the  passage  of  the  shock  may 
literally  trigger  this  energy.  In  any  case,  the  subsequent  movement  of  the 
block  (which  may  relieve  the  shock  stress  by  moving  the  block  as  a  whole) 
generates  a  disturbance  of  its  own  along  the  induced  fault  entirely  comparable 
to  the  mechanism  in  an  ordinary  earthquake,  but  of  a  much  smaller  magnitude 
of  energy.  It  is  suggested  that  the  seismic  signals  alleged  to  be  observed  as 

a  result  of  the  Rainier  event  —  if  they  had  any  significance  at  all  —  were  prob¬ 
ably  due  to  these  minor  induced  earthquakes.  It  is  of  considerable  reassur¬ 
ance,  therefore,  in  worrying  about  earthquakes  from  underground  explosions, 
to  note  that  earthquakes  did  occur  but  were  very  minor  events.  It  is  entirely 
probable  that  once  a  te st  site  has  been  tffired  in,n  by  relieving  the  eons-old 
accumulation  of  stresses  along  planes  of  structural  weakness,  no  further  se¬ 
rious  disturbances  would  result  from  subsequent  shots,  evenslike  these  very 
minor  events. 
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3e  Finally,  it  is  noted  that  the  typical  disturbance  produced  by  the  shock 
wave  itself  is  of  short  duration  and  of  an  entirely  different  nature  from  the 
train  of  waves  associated  with  slippage  along  a  fault  block,  which  have  a  rel¬ 
atively  long  time  duration.  There  is  no  strong  reason  to  suppose  that  the  in- 
itial  compression  wave  in  rock  is  different  from  any  other  shocks  a  sharp 
rise  followed  by  a  more  or  less  rapid  and  asymptotic  decay  back  to  ambient 
pressure.  In  the  low  pressure  approximations,  the  time  decay  of  pressure 
just  behind  the  front  in  such  a  wave  can  be  derived  explicitly  in  terms  of  the 
peak  pressure,  the  time  t  from  the  source,  and  the  properties  of  the  mate¬ 
rial  according  to  ~ 


9p  Po  co 

~w  =  ~ 


(6) 


Here  p0  is  the  density  and  Cq  is  the  sound  velocity.  Even  gross  considera¬ 
tions  with  this  equation  will  show  that  the  npositive  duration’1  for  pressure  in 
the  shock  wave  is  likely  to  be  of  the  order  of  milliseconds;  this  duration  is 
probably  too  short  for  the  average  seismic  instrument  to  resolve;  what  the 
instrument  saw  was  the  mountain  settling  slowly  after  its  initial  jolt. 

VII.  Channeling 

By  channeling  we  mean  that  the  tendency  for  shock  energy  to  propagate 
preferentially  into  rarified  media  instead  of  dense  media  is  so  pronounced  that 
the  energy  can  be  directed  almost  at  will  by  means  of  open  channels. 

The  basic  hydrodynamic  reasons  for  channeling  can  be  simply  shown. 
Compare  the  rate  of  work  W  per  unit  area  and  time  done  by  a  shock  front  on 
the  gas  compared  with  that  done  on  a  solid  ^medium,  assuming  the  configura¬ 
tion  shown  in  Fig.  3.  5.  The  rate  of  work  W  in  any  medium  is  given  by  the 
product  W  =Pu,  a  piston  moving  with  the  force  P  at  the  material  velocity  u. 

If  the  peak  pressure  at  the  shock  front  at  a  given  time  is  the  same  in  the  two 
media  —  and  the  hydrodynamics  will  nearly  always  guarantee  that  it  is  —  then 
the  rate  of  work  in  either  medium  is  described  by  the  equation: 

W  =  Pu  =  P  (P  -  P0)(V0  -  V). 

Under  the  condition  that  the  shock  pressure  is  equal  across  the  interface,  the 
ratio  of  the  rate  of  work  in  each  medium  is  readily  obtained  as 


W 


gas  _ 


W 


solid 


A V\ 

/solid 


(8) 


3  4 

A  factor  of  10  to  1  0  or  so  occurs  in  the  density  ratio  betwen  the  two  media 
and  is  alone  sufficient  to  introduce  a  factor  of  30  to  100  for  energy  propaga¬ 
tion  in  favor  of  the  gaseous  phase,  even  at  high  pressures  where  the  compres¬ 
sion  in  the  solid  is  comparable  to  that  in  the  gas.  However,  over  a  large 
range  of  pressures,  V/V0  in  the  gas  is  so  substantially  less  than  1  so  that  the 
corresponding  term  in  the  numerator  is  practically  1.  At  the  same  time  the 
density  change  aV/V0  in  the  solid  is  very  small  and  this  small  fractional 
change  in  specific  volume  in  the  solid  again  introduces  many  factors  of  10  to 
the  relative  rate  of  work  in  favor  of  the  gaseous  phase.  At  low  enough  pressures 
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Fig.  30  5.  Strong  shock  in  a  gas  channel. 


the  relative  rate  of  work  in  Eq.  (8)  approaches  the  ratio  of  the  familiar  "acous¬ 
tic  impedance^" 

^gas  _  (Po  co)  solid  ,  . 

W  (pc')  9 

solid  oygas 

which  is  a  ratio  like  10,  000  in  favor  of  propagation  in  the  gas. 

Such  considerations  show  the  enormous  tendency  for  the  shock  energy  to 
literally  spill  out  into  the  gaseous  phase  —  a  factor  of  1, 000  at  most  pressure 
levels  —  with  a  corresponding  depletion  of  the  shock  energy  in  the  solid  or 
liquid  phases. 

VIII.  Controlled  Excavations 


Despite  the  w*ell-known  condition  that  the  sound  velocity  in  solid  and  liquid 
media  may  be  five  to  ten  times  the  sound  velocity  in  air,  a  reverse  situation 
occurs  at  nearly  all  pressures  great  enough  to  be  of  practical  interest  for  ex¬ 
cavations. 

A  strong  tendency  for  the  shock  in  the  gaseous  media  to  outrun  the  shock 
in  the  dense  media  is  indicated  by  considering  the  equation  for  shock  velocity, 
given  as  Eq.  (5).  Assuming  we  have  equality  in  pressure  across  the  inter¬ 
face  as  in  Fig.  5,  the  ratio  of  the  shock  velocity  in  each  media  would  be  given 
by 


U 

gas 


^  solid 


(10) 


When  the  shock  strength  is  high  enough  to  achieve  appreciable  compression  of 
the  solid  (which  is  widely  significant  if  the  solid  contains  any  appreciable  air 
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fraction)  then  the  shock  velocity  in  air  outruns  the  shock  velocity  in  soil  by  the 
square  root  of  the  density  ratio*  A  factor  of  about  30  to  1  applies  for  many 
pressures  of  interest*  At  low  enough  pressures  the  shock  velocity  ratio  must 
obviously  approach  the  ratio  of  sound  velocities  in  each  medium*  a  ratio  in  favor 
of  increased  shock  velocity  in  the  solid  medium.  However*  over  the  large 
range  of  pressures  in  which  the  compression  in  solid  is  small  enough  to  be 
considered  linear  (although  much  too  strong  for  the  linear  approximation  to 
apply  in  air),  the  ratio  of  shock  velocities  behaves  as 


U  .  Cn  ,  air 


U  Cn  ,  solid 

solid  0 


(11) 


From  this  equation  it  is  clear  that  a  shock  in  air  will  outrun  the  shock  in  solids 
at  pressures  above  30  atmospheres  or  roughly  500  psi*  Close  in*  the  ratio 
reaches  the  limit  like  30  to  1  suggested  by  the  square  root  of  density  ratios 
alone  in  Eq*  (10)o  This  result  means  that  the  angle  of  inclination  of  the  ground 
shock  in  Fig*  5  is  very  flat;  the  shock  in  ground  moves  l/30  the  distance  out¬ 
ward  from  the  channel  as  the  front  moves  forward  in  the  channel*  For  prac¬ 
tical  purposes*  the  result  is  a  nearly  cylindrical  jet* 

As  a  result  of  such  comparisons  of  shock  velocity  and  with  considerations 
of  the  detailed  material  velocity  in  each  of  the  solid  and  gaseous  phases,  it  is 
possible  to  derive  some  general  criteria  for  the  conditions  under  which  such 
a  channel  may  be  stemmed  or  plugged*  The  details  are  more  complex  than 
can  be  shown  here  but  the  results  are  shown  qualitatively  in  Figs*  3*  6a  and 
3.  6b*  Figure  3*  6a  represents  the  befdre  and  after  pictures  for  an  open  chan¬ 
nel*  Categorically*  any  open  channel  in  a  radial  direction  to  the  device  will 
always  open  up*  The  criteria  for  stemming  depend  upon  the  physical  proper¬ 
ties  of  a  plug  shown  in  Fig*  3*  6b*  But  a  general  rule-of-thumb  appears  that 
a  plug  with  the  same  length  as  the  diameter  of  the  channel  will  usually  be  suf¬ 
ficient  to  stem  the  channel  as  shown*  The  hydrodynamic  reasons  are  as  fol¬ 
lows. 

When  the  shock  first  emerges  from  the  far  side  of  the  plug,  an  enormous 
rarefaction  in  shock  strength  occurs,  reducing  the  transmitted  shock  in  the 
channel  by  a  factor  roughly  equal  to  the  density  ratio  between  the  gas  and  the 
solid  phase.  A  corresponding  rarefaction  moves  backwards  through  the  plug* 
The  pressures  in  the  solid  phase  are  not  similarly  reduced  and  the  walls  of 
the  channel  begin  to  move  in  toward  the  center,  closing  the  channel*  If  the 
plug  is  too  thin,  the  subsequent  rarefaction  back  to  the  device  chamber  occurs 
with  sufficient  rapidity  that  the  device  energy  breaks  through*  But  if  the  plug 
is  long  enough  the  channel  may  already  be  closed  before  the  breakthrough  can 
occur*  The  net  result  is  a  continual  closure  of  the  channel*  At  the  same  time 
it  will  be  clear  from  the  figure  that  a  pinch  jet  will  almost  always  move  down 
the  channel,  ahead  of  the  ground  shock  because  of  the  continual  closure  of  the 
air  channel  behind  it*  The  process  is  somewhat  analogous  to  the  boyhood  trick 
of  shooting  watermelon  seeds  by  squeezing  them* 

Surprisingly  enough,  despite  the  very  much  lower  efficiency  of  a  nuclear 
explosion  with  spherical  symmetry  in  a  homogeneous  soil,  the  channeling  con¬ 
cept  provides  the  means  to  make  nuclear  explosions  more  efficient  than  TNT 
for  many  excavation  purposes*  The  argument  is  as  follows*  Because  of  the 
close  similarity  of  density  between  TNT  and  most  rock-like  materials,  the 
initial  shock  transmitted  outward  from  a  tamped  TNT  charge  is  of  the  same 
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Fig,  3.  6a.  Open  channel  (upper  figures). 

Fig.  3.  6b.  Stemmed  channel  (lower  figures). 


order  as  the  internal  pressure  initially  in  the  high  explosive,  about  300,  000 
atmospheres.  This  results  in  the  shattering  by  TNT  with  a  large  loss  in 
energy  due  to  crushing  and  waste  heat  within  the  rock,  although,  of  course, 
not  comparable  to  the  waste  heat  of  a  tamped  nuclear  explosion.  It  is  difficult 
to  exploit  channeling  with  TNT  even  with  an  open  channel  because  the  large 
mass  effect  of  the  TNT  is  itself  almost  equivalent  to  the  stemming  plug  de¬ 
scribed  above.  No  such  mass  effects  exist  for  nuclear  explosions;  by  the 
channeling  process  it  appears  possible  to  channel  the  energy  into  very  large 
volumes  of  air,  perhaps  reducing  the  pressures  down  to  levels  like  500  psi 
quoted  above,  at  least  down  to  the  crushing  strength  of  the  rock,  or  its  yield 
strength  in  situ.  Thus,  because  the  pressure  can  be  degraded  quickly,  the 
energy  from  a  nuclear  explosion  can  be  made  to  lift  the  overburden  with 
greater  efficiency  than  is  possible  with  TNT. 

The  area  of  a  cylindrical  shock  front  in  soil  (2trrf)  relative  to  the  area  of 
the  channel  itself  (trr^),  together  with  the  ratio  i/r  =  30,  is  sufficient  to  indicate 
that  the  area  of  the  ground  shock  is  about  60  times  the  area  of  the  shock  at  the 
front  of  the  jet.  This  means  that  despite  the  larger  rate  of  work  per  unit  area 
in  air  a  large  fraction  of  the  energy  from  an  explosion  can  be  directed  through 
the  jet  to  work  the  ground:  the  explosion  behaves  kinetically  like  a  line  charge. 

The  argument  has  often  been  advanced  that  an  open  radial  channel  from 
an  explosion  will  close  itself  because  of  the  greater  sound  velocity  in  dense 
media  than  in  air.  That  this  argument  is  invalid  at  pressures  of  interest  is 
clear  from  Eq.  (10)  and  (11).  Furthermore,  in  a  porous  rock  or  soil,  the 
sound  velocities  may  actually  be  comparable.  In  any  case,  Fig.  3.4a  shows 
that  the  shock  velocity  in  a  crushable  dense  medium  can  be  quite  low  even 
where  the  pressure  exceeds  the  crushing  strength.  Hence,  the  argument  for 
automatic  closure  applies  only  at  pressure  levels  which  are  too  low  for  the 
channel  to  crush  down.  The  argument  for  closure  of  an  open  radial  channel 
is  tantamount  to  tjie  incredulous  statement  that  a  pipe  leading  from  a  suddenly 
pressurized  boiler  will  implode  because  of  the  high  sound  velocity  in  steel. 
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4.  SUBSURFACE  EARTH  MOTION 
William  R.  Perret 
Sandia  Corporation 

Subsurface  earth  motion  is  the  time  history  of  transient  mechanical  effects 
produced  by  large  contained  explosions  and  observed  within  the  containing 
medium  in  the  region  extending  from  just  outside  the  explosion  cavity  to  an 
adjacent  free  surface  or  to  the  region  of  essentially  pure  elastic  reactions- 
The  studies  performed  in  this  region  provide  direct  knowledge  of  reaction  of 
the  medium  to  explosive  energy.  Such  knowledge  is  fundamental  to  under¬ 
standing  the  mechanics  of  containment  and  its  antithesis,  cratering. 

The  instruments  for  observing  such  transient  effects  consist  primarily  of 
accelerometers,  supplemented  by  several  specially  developed  velocity  gauges 
and  stress  or  strain  gauges.  Chief  dependence  has  been  placed  on  accelerom¬ 
eters,  not  because  acceleration  is  the  most  desirable  parameter  for  analysis, 
but  because  gauges  suitable  for  measuring  velocity  or  displacement  of  the  mag¬ 
nitudes  and  durations  anticipated  are  only  now  emerging  from  proof  test  phases 
of  development. 

Programs  designed  to  explore  motion  in  the  rock  surrounding  the  Rainier, 
Tamalpais,  and  Evans  shots  of  1957  and  1958  were  conducted  by  several  agen¬ 
cies.  The  Rainier  studies,  undertaken  by  Stanford  Research  Institute  (SRI), 
Sandia,  and  Engineer  Research  and  Development  Laboratories  (ERDL),  Fort 
Belvoir,  included  acceleration  measurements  on  vertical  and  horizontal  radii. 
Figure  4.  1  shows  the  relative  positions  of  the  instrumentation  for  the  Rainier 
explosion.  SRI  and  Sandia  instrumented  borings  adjacent  to  surface  zero. 

ERDL  and  Sandia  instrumentation  was  installed  in  the  access  tunnel.  These 
projects  met  with  only  qualified  success,  primarily  because  the  instrument 
ranges  chosen  were  in  some  cases  too  low,  and  recording  circuitry  limited 
recorded  peak  signals. 

Evans  studies  were  designed  to  overcome  some  of  the  difficulties  en¬ 
countered  by  the  Rainier  projects.  Two  vertical  instrumentation  systems  ex¬ 
tending  from  about  200  ft  from  the  zero  point  to  the  surface  were  installed  by 
SRI  and  Sandia.  Again  dependence  was  placed  primarily  on  accelerometers, 
but  SRI  included  several  newly  developed  velocity  gauges  in  their  array. 
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Fig.  4.  1.  Plan  and  elevation,  Rainier  tunnel. 
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Fig.  4.  2.  Distribution  of  vertical  instrumentation  for  Evans. 
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Figure  4 .  2  shows  the  distribution  of  vertical  instrumentation.  Horizontal  ra¬ 
dius  instrumentation  was  installed  adjacent  to  the  access  tunnel  by  ERDL  and 
Ballistic  Research  Laboratory  (BRL)  for  the  Tamalpais  and  Evans  event. 

Success  of  these  programs  was  again  seriously  limited,  this  time  principally 
by  the  fact  that  the  energy  derived  from  the  Evans  event  was  only  about  1%  of 
that  anticipated  and  measuring  and  recording  systems  were  operating  in  many 
instances  close  to  the  limit  of  resolution.  Precision  of  the  results  is  conse¬ 
quently  low. 

All  acceleration  data  taken  within  the  rock  during  Rainier  and  Evans  events 
were  radial.  These  data  from  both  vertical  and  horizontal  radii  are  shown 
plotted  against  range  in  Fig.  4.  3.  Two  plots  are  shown,  one  of  raw  data,  the 
second  of  scaled  data.  Dimensional  analysis  indicates  that  accelerations 
should  follow  a  relationship  of  the  form 

AW1/3  afRW'1/3)  n  , 

where  AW1/3  is  scaled  acceleration  and  RW'1/3  is  scaled  range.  The  data 
suggest  that  the  exponent  should  be  -4  out  to  a  range  of  about  700  ft,  beyond 
which  it  changes  toward  -2,  There  is  an  obvious  discrepancy  in  the  scaled 
plot  between  the  Rainier  and  Evans  data,  due  in  part  to  the  fact  that  most 
Rainier  peaks  were  low  by  a  factor  greater  than- 2,.  and  the  widths  of  Evans 
records  were  relatively  large  compared  with  deflections. 

The  inverse  fourth  power  relationship  between  accelerations  and  range 
near  an  explosion  appears  valid  for  both  plots  although  it  is  not  evident  that 
scaling  laws  are  valid.  Point- source  theory  implies  that  pressure  varies 
inversely  with  the  cube  of  radial  range,  but  the  radial  derivative  of  pressure 
must  be  proportional  to  acceleration.  Hence,  inverse  variation  of  accelera¬ 
tion  with  the  fourth  power  of  range  indicates  that,  out  to  about  700  feet  from 
an  explosion,  rock  reacts  as  to  a  point  source. 

There  is  a  slope  inversion  in  each  vertical  radius  data  curve  at  ranges 
which  encompass  a  region  within  300  feet  of  the  free  surface.  This  region 
includes  the  welded-tuff  caprock  of  the  Rainier  mesa.  Peaks  decrease  with 
range  up  to  a  point  about  250  feet  below  the  surface.  Above  this  point  accel¬ 
erations  increase.  Time  relationships  are  such  that  surface  reflection  effects 
cannot  cause  the  increase.  In  fact,  one  record,  AVS-5  Fig,  4.  4,  shows  a  re- 
fleeted  peak  following  the  major  one;  reflected  signal  arrivals  indicated  on  other 
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Fig.  4.  4.  Acceleration  —  time  records,  vertical  radius 
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records  are  derived  times  based  on  travel  times  to  surface  and  return,,  Ex¬ 
amination  of  propagation  velocity  data  and  rock  densities  shows  that  the  velo¬ 
city  decreases  through  a  factor  of  3  more  or  less  continuously  from  the  base 
of  the  caprock  toward  the  free  surface.  Densities  vary  by  not  more  than  20%. 
Thus  seismic  impedance,  the  product  of  propagation  velocity  by  density,  de¬ 
creases  toward  the  surface  from  the  base  of  the  caprock  and  particle  displace 
ments  and  accelerations  increase  accordingly  as  a  result  of  the  amplitude  re¬ 
lationship 


X 


2 


ZpjCj 

PlCl  +  p2c2 


9 


where  p!s  are  densities,  C?s  propagation  velocities,  and  X*s  amplitudes  or 
displacements  transmitted  from  medium  1  into  medium  2.  p^C-^  grea^er 

than  p2^2° 

This  enhancement  effect  is  more  pronounced  in  the  Rainier  data  than  in 
the  Evans  results,  but  it  appears  in  both  sets  of  records  at  about  the  same 
depth.  The  Evans  results  involve  small  signals  near  the  limit  of  resolution 
of  the  gauges,  and  energies  too  small  to  produce  gross  motion  and  partings 
of  the  caprock.  Therefore,  the  enhancement  in  the  Evans  data  does  not  fol¬ 
low  as  consistent  a  pattern  as  do  the  Rainier  data.  Only  at  the  free  surface 
were  accelerations  of  markedly  increased  magnitude  noted  in  Evans  records. 

Motion  of  the  caprock,  as  described  by  the  four  uppermost  acceleration 
records  from  Rainier  in  Fig.  4.  4,  indicate  upward  motion  followed  by  a  period 
of  free  fall.  These  records  are  a  composite  of  Sandia  and  SRI  data.  Free 
fall  ends  for  each  gauge  in  a  rebound  pulse,  implying  that  the  floating  rock 
came  back  to  earth.  Comparison  of  intervals  between  first  arrivals  and  ar¬ 
rivals  of  rebound  indicate,  by  differences,  at  least  three  partings  between  a 
depth  of  250  feet  and  the  surface. 

Particle  velocities  have  been  derived  from  both  Rainier  and  Evans  ver¬ 
tical  data.  Those  from  Rainier  were  derived  by  integration  of  the  accelera¬ 
tion-time  records,  but  those  from  Evans  were  measured  directly  by  the  SRI 
gauges  and  were  checked  by  integration  of  corresponding  acceleration  data. 

The  plot  of  peak  velocities  versus  scaled  range,  shown  in  Fig.  4.  5,  indicates 
that  below  the  caprock  discontinuity  velocities  vary  with  the  inverse  square 
of  scaled  range.  Scaling  relations  for  velocity  and  range  should  have  the  form 

V  a  (RW-1/3)11  . 
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The  graph  shows  that  the  exponent,  n,  is  -2,  Slope  of  the  velocity  curves  also 
undergo  inversion  within  the  caprock. 

Finally,  a  few  observations  of.  stresses  and  strains  are  available  from 
the  BRL  and  ERDL  Tamalpais  instrumentation.  The  five  observations  of 
radial  stress  by  BRL  ranged  from  about  1100  psi  at  100  feet  downward  to  18 
psi  at  1100  feet  from  the  57-ton  shot.  Tangential  stresses  were  about  60% 
of  radial  stresses.  The  few  data  available  indicate  that  both  components  of 
stress  observed  vary  inversely  with  the  cube  of  the  radial  range,  again  im¬ 
plying  that  within  this  region  the  rock  reacts  as  to  a  point  source. 

ERDL  strain  data  indicate  that  peak  transient  strains  are  compressions 
which  vary  inversely  with  the  square  of  scaled  radius.  Compressive  strains 
observed  ranged  downward  from  about  0.  1%  at  100  ft. 

The  data  available  concerning  motion  within  the  earth  near  a  contained 
nuclear  explosion  imply  that  the  rock  reacts  as  to  a  point  source  in  a  homo¬ 
geneous  medium  over  the  range  of  distances  beyond  the  gross  deformation  of 
the  cavity  out  to  at  least  700  ft.  Scaling  laws  do  not  appear  to  apply  to  peak 
acceleration  data,  but  circumstances  which  limit  precision  of  the  data  may 
account  for  many  of  the  discrepancies  in  scaling. 
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5.  SURFACE  EFFECTS 
Lu  M.  Swift 

Stanford  Research  Institute 


The  surface  effects  of  a  deep  underground  explosion  are  not  merely  a  con¬ 
tinuation  of  the  subsurface  effects.  Even  in  the  case  of  the  ideal  but  never  at¬ 
tained  ''infinite,  homogeneous,  elastic,  etc.  "  solid,  the  stress  wave  arriving 
at  surface  ground  zero  is  reflected  as  a  negative  of  its  former  self.  That  is, 
a  compressional  wave  is  reflected  as  a  tensile  wave,  and  vice  versa.  What 
is  not  so  obvious  is  that  this  reversal  in  both  sign  and  direction  of  travel  of 
the  stress  wave  requires  that  the  particle  motion  is  not  reversed.  Thus  the 
meeting  of  the  remains  of  the  incident  wave  with  its  reflections  near  the  sur¬ 
face  causes  the  stresses  to  cancel,  but  causes  an  amplification  of  the  particle 
motions. 

At  points  other  than  ground  zero,  where  the  incidence  of  the  wave  with 
the  surface  is  other  than  perpendicular,  these  cancellations  and  amplifications 
are  not  direct,  and  changes  occur  in  the  directions  of  stresses  and  movements. 
The  general  effect  is  to  translate  more  and  more  of  the  energy  into  a  horizon¬ 
tal  direction,  that  is,  parallel  with  the  free  surface.  When  the  medium  is 
layered,  as  real  conditions  usually  are,  this  tendency  may  be  offset  by  the 
tendency  of  the  energy  to  follow  a  refracted  path  in  reaching  the  surface,  so 
that,  beyond  a  critical  angle,  the  angle  of  arrival  is  constant,  for  high  fre¬ 
quency  waves  at  least.  Further  complications  are  added  by  the  presence  of 
faults  and  joints,  and  by  reflections  from  moderately  deep  layers. 

At  great  distances,  most  of  the  high  frequency  energy  is  lost,  and  the  re¬ 
maining  signals  are  low  frequency  and  consequently  of  long  wavelength.  These 
long  waves  simply  cannot  "see"  the  smaller  discontinuities  in  the  earth,  and 
are  affected  primarily  by  the  major  features  of  the  earth® s  structure.  That 
region  belongs  to  the  seismologist,  however,  and  I  am  sure  he  appreciates 
my  oversimplification  of  his  problems.  With  apologies,  then,  we  will  with¬ 
draw  to  the  close-in  surface  effects,  defined  as  those  at  distances  no  greater 
than  a  few  times  the  depth  of  burial  of  the  shot.  * 

On  shot  Rainier,  records  were  obtained  of  the  vertical  acceleration  at 
ground  zero,  and  of  vertical  and  horizontal  accelerations  at  six  other  ranges 


scussion  applies  to  burial  sufficient  for  containment. 
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out  to  about  a  half  mile.  At  most  of  these  points,  records  were  obtained  of 
the  strain  at  the  surface  as  well.  A  point  of  perhaps  maximum  interest  on  the 
surface  is  surface  ground  zero.  It  is  usually  the  closest  surface  point  to  the 
shot,  and  hence  the  one  where  maximum  effects  are  to  be  expected.  Also,  as 
mentioned  previously,  it  is  a  region  where  the  mechanisms  of  reflection  are 
simple.  Mr.  Perret  has  described  how  the  original  violent,  crushing  shock 
wave  was  modified  by  dissipation  and  geometric  spreading  to  a  more  complex 
wave  of  longer  duration  and  lower  peak  magnitude.  Beyond  a  few  hundred 
feet,  it.  differed  little  from  a  purely  elastic  wave.  As  this  wave  approached 
the  surface,  the  motions  were  amplified  by  reflections  from  the  surface,  so 
that  they  were  greater  at  the  surface  than  at  a  depth  of  250  feet,  in  spite  of 
the  greater  distance  from  the  shot.  Figure  5.  1  shows  the  motion  of  the  sur¬ 
face,  expressed  in  physical  quantities.  The  vertical  acceleration,  shown 
here  as  a  function  of  time,  is  seen  to  exhibit  a  sharp  pulse,  with  a  peak  of 
about  6  g  and  a  duration  of  50  msec,  followed  by  a  flat  region  of  -1  g  magni¬ 
tude  and  about  400  msec  duration.  (This  -1  g  is  with  reference  to  a  normal 
zero  established  with  the  accelerometer  being  supported  by  the  earth  -  -1  g 
actually  represents  a  free  fall  condition.  )  Following  this  is  another  pulse 
similar  to  the  first.  These  data  can  be  converted  into  terms  of  velocity  by 
integration  of  the  acceleration  with  respect  to  time.  The  first  acceleration 
pulse  produced  a  velocity  upward  of  about  6  fps.  During  the  -1  g  period,  this 
velocity  gradually  dropped  to  zero  and  on  to  about  the  same  velocity  in  a  down¬ 
ward  direction.  At  the  time  of  the  second  bump,  most  of  the  velocity  was  lost, 
and  the  surface  stopped  moving,  except  for  a  moderate  and  more  or  less  ran¬ 
dom  vibration. 

Integrating  once  more,  to  obtain  displacement,  we  see  a  parabolic  tra¬ 
jectory  which  begins  to  look  familiarj  it  shows  that  what  happened  was  simply 
that  the  surface  was  tossed  into  the  air  about  nine  inches,  and  landed  again 
about  where  it  started.  The  dotted  curve  of  displacement  is  one  obtained  from 
motion  picture  photography  of  the  surface,  taken  with  telephoto  lens  from  a 
point  about  4000  feet  away.  The  agreement  is  good. 

These  records  alone  do  not  demonstrate  it  clearly,  but  in  combination 
with  subsurface  data  they  show  that  it  was  not  only  the  very  surface  which  was 
in  free  fall,  but  a  layer  some  hundred  feet  or  more  thick.  In  fact  there  appear 
to  have  been  several  horizontal  parting  planes  at  depths  down  to  250  feet  or  more. 
This  effect  is  familiar  and  similar  to  that  o'f  "spalling,  "  frequently  observed 
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in  small  explosions,  but  here  it  is  complicated  by  the  fact  that  the  formation 
is  held  together  more  by  the  force  of  gravity  than  by  any  tensile  strength  of 
the  material. 

Most  of  us  have  had  the  experience  of  unexpectedly  stepping  off  a  curb  in 
the  dark,  It's  startling,  frequently  painful,  but  seldom  fatal.  That  is  about 
the  sensation  that  would  have  been  experienced  by  an  observer  standing  at 
ground  zero  when  Rainier  was  fired,  except  that  he  would  have  felt  two  blimps, 
each  of  only  about  5  or  6  g.  But  that  was  Rainier,  I  assure  you  that  my  hy¬ 
pothetical  observer  would  have  been  much  more  embarrassed  had  he  attempted 
to  sit  on  some  of  those  other  underground  shots. 

One  of  those  shots  was  Blanca,  This  was  a  larger  shot  at  about  the  same 
depth  as  Rainier,  so  one  would  expect  more  surface  effects.  Unfortunately, 
this  shot  was  under  the  slope;  for  this  and  other  reasons  close-in  measure¬ 
ments  were  not  feasible,  but  information  is  available  on  the  movement  of  a 
couple  of  points  on  top  of  the  mesa  but  near  ground  zero.  Figure  5.  2  gives 
the  vertical  displacement  as  a  function  of  time,  of  a  station  (1205,  18)  dis¬ 
placed  horizontally  about  750  feet  from  ground  zero.  The  dots  are  the  actual 
data  points  as  read,  and  the  solid  curve  is  the  smoothed  interpretation  applied 
by  the  analysts,  showing  a  peak  upward  displacement  of  some  33.  4  inches, 
occurring  about  436  msec  after  the  first  motion  was  observed.  I  have  drawn 
the  dotted  curve,  which  is  the  calculated  curve  for  an  object  in  free  fall,  start¬ 
ing  with  an  upward  velocity  of  about  160  inches  per  second.  I  think  this  agree¬ 
ment  is  about  as  good  as  that  in  Fig.  5.  1,  so  I  suggest  that  the  phenomena  are 
similar.  In  any  case,  a  velocity  of  160-180  ips  was  attained  in  about  20  msec, 
calling  for  an  acceleration  of  at  least  25-g  peak. 

Figure  5.  3  represents  similar  data  from  a  station  (1205.  17)  about  900  feet 
from  ground  zero  and  somewhat  farther  from  the  mesa  edge.  Its  peak  upward 
displacement  is  about  20  inches,  and  the  apparent  curve  does  not  fit  a  free-fall 
parabola  well,  so  at  this  range  the  formation  must  have  held  together. 

On  this  curve,  like  the  others,  the  photographic  data  stop  when  the  shock 
reached  the  camera.  In  this  case,  the  surface  has  returned  practically  to 
zero,  and  still  has  a  downward  velocity  of  near  40  ips.  For  this  one  station 
on  this  shot  we  have  pre-shot  and  post- shot  survey  data  which  show  that  it 
finally  stopped  about  19  inches  below  its  original  position.  This  is  partly  in 
accord  with  data  from  permanent  displacement  surveys  on  Rainier,  where  the 
permanent  displacements  were  found  to  be  downward  and  inward  toward  ground 
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zero.  On  Rainier,  though,  the  displacements  were  small,  none  over  2  inches 
downward.  The  permanent  displacement  of  the  station  1205.  18  on  Blanca  was 
large  but  not  meaningful.,  It  was  near  the  edge  of  the  cliff,  and  fell  off. 

Returning  to  Rainier,  we  have  data  on  the  ground  motion  as  a  function  of 
time  at  a  number  of  ranges.  As  might  be  expected,  these  fall  off  rapidly  with 
distance.  Figure  5.  4  shows  a  profile,  with  the  vertical  scale  greatly  exag¬ 
gerated,  of  the  vertical  displacements  at  different  times.  Superimposed  is  a 
curve  which  represents  the  calculated  shape  of  the  hump  produced  in  an  elastic 
medium  by  a  pressure  center  at  the  depth  of  shot  Rainier.  The  difference, 
though  marked,  is  not  surprising  in  view  of  the  inelastic  and  layered  condition 
of  the  medium  and  the  fact  that  we  are  comparing  a  dynamic  result  with  a  static 
calculation. 

In  forming  the  hump  shown  here,  an  appreciable  tensile  strain  was  pro¬ 
duced  in  the  near- surface  formations.  Rocks  of  the  general  types  found  at 
the  Nevada  Test  Site  simply  have  no  tensile  strength  in  the  gross  form.  An 
individual  small  rock  may  have  appreciable  strength,  but  in  the  rock  in  place 
there  are  faults  and  joints  at  frequent  intervals,  so  that  when  a  tensile  strain 
is  applied,  cracks  form  every  few  feet.  On  Rainier,  these  cracks  were  small 
—  they  probably  did  not  open  more  than  two  or  three  inches  at  their  widest,  and 
they  all  re-closed  almost  immediately. 

Figure  5.  5  shows  the  strain  at  the  surface  at  ground  zero,  as  measured 
by  a  gauge  with  a  total  span  of  50  feet.  The  maximum  observed  was  about 
four  parts  per  thousand,  or  0.4%.  In  the  total  span  of  50  ft,  this  represents 
an  elongation  of  about  2-1/2  inches.  Post- shot  observations  showed  that 
several  cracks  were  covered  by  this  span,  so  that  at  their  widest  these  cracks 
could  not  have  been  open  more  than  a  couple  of  inches. 

On  the  other  hand,  some  of  the  cracks  formed  in  the  surface  by  shot 
Blanca  must  have  been  quite  large,  even  when  we  do  not  include  those  quite 
near  the  edge  of  the  cliff.  A  peak-recording  strain  gauge  with  a  50-foot  span, 
located  at  station  1205.  17,  showed  a  peak  strain  of  about  nine  parts  per  thou¬ 
sand,  about  twice  the  maximum  recorded  on  Rainier,  even  though  this  station 
was  about  900  feet  from  ground  zero.  Figure  5.  6  is  a  photograph  of  the  post¬ 
shot  surface  conditions  near  that  station.  Note  that  there  is  considerable  ver¬ 
tical  displacement  across  these  faults.  There  is  reason  to  believe  that  at- 
least  some  of  the  cracks  opened  quite  wide,  but  most  of  them  re-closed  al¬ 
most  immediately. 
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In  discussing  the  surface  phenomena,  I  have  avoided  much  discussion  of 
movements  on  the  slope  of  the  mesa,  since  it  would  seem  that  many  such  move¬ 
ments  can  be  utterly  disproportionate  to  the  excitation  that  caused  them.  The 
material  of  this  slope  is  quite  naturally  lying  at  the  maximum  slope  which  is 
stable-  Occasional  slides  occur  naturally,  from  weathering,  etc-  It  seems 
inappropriate,  then,  to  pay  much  attention  to  the  quantitative  displacements 
of  portions  of  the  cliff  which  are  displaced  by  the  shock  but  which  fall  or  roll 
great  distances  down  the  slope-  On  shot  Rainier,  an  effort  was  made  to  meas¬ 
ure  the  transient  motion  of  a  series  of  outcrops  along  the  slope,  but  the  gauges 
were  either  destroyed  or  disturbed  by  the  rock  slides,  and  the  data  obtained 
were  not  very  helpful- 

In  summary,  then,  we  may  describe  the  close-in  surface  effects  of  a 
deep  underground  explosion  as  follows: 

a-  The  first  arrival  is  a  relatively  sharp  pulse  of  upward  acceleration 
which. imparts  an  upward  velocity  to  the  surface  formations- 

be  Under  the  influence  of  this  velocity,  the  surface  continues  upward, 
being  decelerated  only  by  the  force  of  gravity  or  a  little  more,  so  that  there 
is  an  appearance  of  a  slow  upward  stheave-  n 

ce  During  this  upward  movement,  cracks  tend  to  form  in  the  surface, 
but  there  is  little  tendency  toward  lateral  or  vertical  displacements  across  the 
cracks,  unless  a  major  discontinuity  of  the  formation  intersects  the  surface- 
dn  After  reaching  its  peak,  the  nhumpn  falls  back  to  its  original  position, 
landing  with  a  second  shock-  At  this  time,  most  of  the  cracks  close,  except 
for  those  which  are  somehow  wedged  apart  by  loose  rocks-  In  some  cases, 
there  are  distortions  producing  vertical  or  lateral  displacements  across  the 
cracks- 

e-  If  venting  occurs,  it  is  usually  in  the  late  stages  of  c  or  d,  unless 
the  burial  is  really  to  be  classed  as  shallow-  In  case  of  full  containment, 
d  is  followed  by  more  or  less  irregula.r  vibrations,  largely  caused  by  reflec¬ 
tions  or  refractions  from  irregularities  in  the  formation- 

f-  The  final  position  of  points  on  the  surface  near  ground  zero  has  been 
found  lower  than  their  original  position-  This  is  not  necessarily  a  general 
case  —  I  have  observed  permanent  humps  from  small  contained  shots,  but 
this  may  not  be  possible  with  large  shots- 

g-  The  effects  fall  off  rapidly  with  distance  from  ground  zero-  The  peak 
accelerations  tend  to  be  proportional  to  the  inverse  fourth  power  close-in, 


-70- 


UCRL-5675 


shifting  rapidly  to  the  inverse  cube  of  the  range;  the  peak  displacement  does 
not  fall  off  so  rapidly,  since  the  frequency  changes  as  well.  At  quite  great 
distances,  the  rate  of  decline  decreases,  since  indirect  propagation  effects 
become  prominent. 
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6.  SEISMIC  MEASUREMENTS 
Carl  Romney 

Headquarters,  U.  S.  Air  Force 
I.  INTRODUCTION 

This  report  presents  preliminary  technical  conclusions  based  on  the 
analysis  of  seismograms  recorded  at  more  than  40  stations  at  the  time  of  the 
Rainier  explosion  (September  1957)  and  seismograms  recorded  at  more  than 
90  stations  at  the  time  of  the  Hardtack  underground  nuclear  explosions  (Oc¬ 
tober  1958).  Recordings  are  available  for  five  underground  nuclear  explosions 
ranging  in  yield  from  about  50  tons  to  20  kilotons.  The  basic  information  on 
the  shot  times,  locations  (including  depth  of  burial)  and  yields  of  explosions 
has  been  made  available  to  seismologists  by  the  Atomic  Energy  Commission. 

II.  FIELD  PROGRAM 

When  it  was  recognized  that  the  underground  nuclear  explosions  planned 
for  Operation  Hardtack  could  provide  important  basic  scientific  information 
of  value  to  seismologists,  the  Air  Force  undertook  a  program  for  the  meas¬ 
urement  of  seismic  effects  from  these  shots.  With  the  assistance  of  a  con¬ 
tractor  (The  Geotechnical  Corporation),  the  Air  Force  was  able  to  install  17 
temporary  seismic  stations. 

The  results  to  be  reported  here  depend  to  a  large  extent  on  the  measure¬ 
ments  from  these  temporary,  stations  supplemented,  where  possible,  by  data 
from  the  U.  S.  Coast  and  Geodetic  Survey  stations  and  from  many  stations 
operated  by  universities.  Because  of  the  vast  amount  of  information  available, 
the  present  report  will  be  limited  to  a  summary  of  the  chief  conclusions  of 
interest  to  seismologists. 

A.  Location  of  Temporary  Stations 

The  stations  were  disposed  at  distances  extending  from  60  kilometers  from 
the  shot  points  to  a  distance  of  slightly  more  than4000  km.  Many  of  the  temporary 
stations  were  equipped  with  vehicles  enabling  them  to  move  to  various  locations 
as  required.  The  temporary  stations  were  situated  generally  along  a  line  ex¬ 
tending  eastward' from  the  Nevada  Proving  Ground  to  Arkansas  and  thence 

*AEC  release  on  Hardtack,  March  10,  1959,  No.  B-39.  Also  refer  to 
Table  I  on  page  81. 
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northeastward  to  Maine  (Fig.  6.  1).  Each  operating  site  and  alternate  sites 
had  been  preselected  by  a  team  of  geologists  who  located  suitable  outcrops  of 
hard  rock  remote  from  obvious  sources  of  man-made  noise.  It  was  planned 
that  five  temporary  stations  would  remain  at  fixed  locations  to  provide  control 
data  on  all  shots,  and  it  was  hoped  that  the  remaining  ten  mobile  stations 
could  be  moved  between  shots  in  order  to  increase  the  amount  of  recorded 
data.  These  mobile  teams  were  in  fact  able  to  move  between  shots  Logan  and 
Blanca.  In  all,  recordings  were  made  by  the  Air  Force  teams  at  28  different 
distances  from  the  shot  point.  Successful  recordings  were  thus  made  at  in¬ 
tervals  of  approximately  100  km  out  to  2300  km,  and  at  intervals  of  about 
250  km  thereafter  to  a  distance  of  4000. 

B.  Equipment 

The  temporary  stations  were  all  equipped  with  Benioff  short-period  ver¬ 
tical  seismographs;  twelve  of  them  were  additionally  equipped  with  two  hori¬ 
zontal  components  which  were  oriented  along  the  transverse  to  the  path  from 
the  shot  point  to  the  station.  Each  station  was  provided  with  a  good  timing 
device  and  means  for  recording  time  signals  from  WWV  on  a  continuous  basis. 
Each  station  also  was  provided  with  equipment  for  calibrating  the  seismo¬ 
graphs  and  for  determining  the  polarity  of  each  instrument  under  operating 
conditions.  In  addition,  special  efforts  were  made  to  check  both  the  polarity 
and  the  calibration  by  independent  methods  at  one  time  or  another  during  the 
course  of  the  observations. 

III.  TECHNICAL  RESULTS 

A.  Longitudinal  Waves 

1 .  Amplitude- Yield  Relation  for  P  Waves 

The  Hardtack  underground  explosions  produced  almost  identical  wave¬ 
forms  at  any  given  station.  It  was  therefore  possible  to  measure  the  ampli¬ 
tude  of  each  shot  at  several  corresponding  points  within  the  wave  train.  The 
separate  measurements  were  then  divided  by  the  corresponding  amplitudes 
recorded  from  Logan  so  that  all  results  are  normalized  to  Logan  amplitude 
equals  1. 0.  The  results  of  these  observations  are  shown  in  Fig.  6.  2.  It  may 
be  seen  that  the  points  described  in  approximate  first  power  relation  between 
amplitude  and  yield.  For  reasons  which  are  not  understood,  the  amplitudes 
from  shot  Evans  were  less  than  one-tenth  of  those  from  Tamalpais,  although 
the  two  shots  released  approximately  the  same  amount  of  energy  (55  and  70 
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tons,  respectively).  For  a  comparison,  the  amplitudes  recorded  on  the  Wood- 
Anderson  torsion  seismograph  at  Tinemaha  are  plotted  on  the  same  graph. 
These  measurements  are  of  the  maximum  amplitudes  in  the  shear  waves  rath¬ 
er  than  in  the  P  wave  group.  It  may  be  seen  that  the  Tinemaha  shear  wave 
amplitudes  have  about  the  same  dependence  on  yield  as  the  amplitudes  of  P 
wave  s .  - 

2.  Amplitude  of  P  Waves  as  a  Function  of  Distance 

The  maximum  amplitude  in  the  first  three  cycles  of  P  waves  is  shown  as 
a  function  of  distance  in  Fig.  6.  3.  The  data  indicated  by  the  solid  circles  are 
from  Blanca  (20  kt).  Logan  data,  indicated  by  open  circles,  were  used  to  help 
define  this  curve  by  multiplying  each  observed  Logan  amplitude  by  2.  36,  which 
was  the  average  amplitude  ratio  of  Blanca  to  Logan  observed  at  fixed  stations. 
Curves  for  other  yields  may  be  constructed  by  using  the  first-power  scaling  law 
mentioned  above. 

At  distances  between  about  200  and  1000  km  the  amplitudes  of  Pn  are  in¬ 
versely  proportional  to  the  cube  of  the  distance,  confirming  measurements 
made  from  earthquakes  and  previously  reported  in  the  literature  of  seismology. 

From  about  1200  to  2000  km  a  strong  P-wave  arrival  was  observed  in  the 
nshadow  zone"  discussed  by  Gutenberg  and  Richter  in  several  papers  on  earth¬ 
quake  magnitudes.  These  large  amplitudes  are  associated  with  longer  period 
waves  and  late  arrivals.  The  abrupt  change  in  the  amplitude -distance  rela¬ 
tionship,  the  travel -time  anomalies,  and  the  change  in  period  imply  that  a 
new  kind  of  wave  begins  to  arrive  at  about  1200  km.  Pn  waves  may  well  pre¬ 
cede  the  new  wave,  but  are  too  small  to  be  detectable  at  these  distances. 

Beyond  2000  km,  the  amplitudes  fall  off  slowly  and  regularly,  except  that 
there  is  a  low  amplitude  region  at  3300  and  3500  km.  Logan  was  not  detected 
at  3500  km  and  Blanca  was  not  detected  at  3300  km.  In  each  case,  measure¬ 
ments  of  microseismic  noise  showed  that  the  P-wave  amplitudes  must  be  less 
than  10  millimicrons.  Figure  6.4  shows  a  selection  of  the  vertical  component 
seismograms  obtained  at  distances  from  200  through  2650  km;  the  figure  in¬ 
cludes  recordings  of  both  shots.  The  seismograms  are  arranged  in  order  of 
distance  and  aligned  in  time  so  that  the  arrivals  should  all  lie  along  a  straight 
line  parallel  to  the  distance  axis  if  the  speed  of  the  P  waves  is  8.  1  km/ sec. 

The  marks  on  the  seismograms  are  at  10- sec  intervals;  upward  motion  on  the 
seismograms  is  compressional  motion  of  the  earth.  It  may  be  seen  that  the 
direction  of  first  motion  is  observed  as  compressional  at  some  stations  and 
rarefactional  at  others.  There  seems  to  be  no  systematic  dependence  upon 
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the  distance.  At  some  stations  (e.  g.  ,  the  station  at  2300  km)  there  is  appar¬ 
ently  clear  rarefactional  first  motion  in  spite  of  a  signal- to-noise  ratio  of  at 
least  10, 

3.  Travel-Times  of  P-Waves 

Figure  6.  5  shpws  the  observed  times  of  P  as  departures  from  a  velocity 
of  8,  1  km/sec.  Last  year  at  the  annual  meeting  of  the  Seismological  Society, 

I  reported  that  the  travel  times  of  P  from  the  Rainier  explosion  fit  the  equation: 

t  -  6.  9  +a/7.9 

between  180  and  1023  km  from  the  source.  The  new  data  may  be  seen  to  agree 
with  that  interpretation.  Beyond  1200  km,  most  of  the  data  show  higher  ve¬ 
locity  arrivals,  which  initially  are  late  with  respect  to  the  formula  mentioned. 
Between  1500  and  2500  km,  the  travel  times  and  the  pulse  shapes  show  mul¬ 
tiple  arrivals  which  may  be  interpreted  to  show  one  or  more  cusps  in  the 
travel-time  curve.  I  am  unable  to  explain  the  arrivals  in  detail  at  the  present 
time.  It  seems  probable  that  new  data  will  be  required  to  understand  the 
propagation  paths.  A  more  detailed  study,  with  more  data,  should  provide 
great  insight  into  the  velocity  structure  just  below  the  earth's  crust.  Beyond 
2500  kilometers,  the  times  approximately  fit  the  Jeffreys-Bullen  and  the 
Gutenberg  1953  travel  time  curves, 

B.  Shear  Waves 

Shear  waves,  arriving  at  a  speed  of  about  3.  5  km/sec,  were  observed  to  a 
distance  of  about  2000  km.  These  waves  have  not  been  analyzed  in  detail  at 
the  present  time.  As  a  general  observation,  however,  it  can  be  stated  that 
the  radial-  and  transverse-horizontal  components  of  motion  are  approximately 
equal  in  magnitude,  and  that  these,  in  turn,  are  approximately  equal  in  size 
to  the  vertical  component  of  motion.  At  some  distances  there  is  the  appear¬ 
ance  of  several  discrete  arrivals  which  might  help  to  confirm  that  these 
waves  are  high-mode  Love  waves  propagated  in  accordance  with  a  method 
suggested  by  Oliver.  Amplitudes  are  shown  in  Fig.  6.  6.  The  fact  that  these 
amplitudes  vary  as  the  inverse  cube  of  the  distance  is  not  consistent  with  the 
interpretation  that  these  are  surface  waves,  however. 

Co  Surface  Waves 

Rayleigh  waves  were  detected  by  long-period  seismographs  at  Berkeley, 
California;  Laramie,  Wyoming;  Fort  Sill,  Oklahoma;  and  Palisades,  New 
York.  In  addition,  there  are  clear  indications  of  Love  waves  at  Berkeley  and 
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Palisades.  These  latter  waves  were  completely  unexpected  from  blasts,  and 
their  presence  on  the  seismogram  is  believed  to  be  due  to  mode  conversion 
near  the  source,  by  processes  which  are  not  understood  at  the  present  time. 

D.  Magnitudes 

Magnitudes  were  first  determined  for  all  shots  recorded  by  Wood-Anderson 
torsion  seismographs.  It  should  be  emphaiszed  here  that  these  magnitudes 
merely  serve  to  compare  the  size  of  the  seismic  waves  from  blasts  with  those 
from  earthquakes;  it  is  not  believed  that  the  energies  released  can  be  compared 
by  this  means.  The  complete  results  are  shown  in  Table  I. 


Table  I.  Magnitude  Summary. 


Station 

Blanca 

Logan 

Rainier 

Palo  Alto 

5.  2 

4.  8 

4.4 

Berkeley 

4.9 

4.  5 

4.  1 

Mineral 

00 

4.4 

4.  3 

Mt.  Hamilton 

5.4 

5.  0 

4.  7 

USF 

4.  9 

4.  5 

4.  2 

Tinemaha 

5.  1 

4.  8 

4.  2 

Pasadena 

4.  8 

4.  4 

4.  0 

Barrett 

3.  9 

3.  7 

Riverside 

4.  7 

4.4 

Woody 

4.  2 

3.  9 

The  magnitude  of  Blanca  was  found  to  be  4.  8  +  0.  4*  based  on  measurements 
from  ten  stations.  The  magnitude  of  Logan  was  found  to  be  4.  4  +  .0.  4#  based 
on  measurements  from  the  same  ten  stations.  Only  seven  stations  recorded 
measurable  amplitude  for  Rainier;  these  gave  a  magnitude  of  4.  25  0.  2*.  The 

magnitudes  for  the  smaller  shots  were  observable  only  at  one  station  and  re¬ 
liance  should  not  be  placed  on  these  single  estimates.  It  is  perhaps  important 
to  note  that  if  the  magnitudes  of  Blanca  and  Logan  are  computed  from  the  same 
seven  stations  that  recorded  Rainier,  the  result  will  be  0.  2  of  a  unit  higher 
than  that  found  from  the  ten  stations.  The  implication  is  Strong. that  had  the 
amplitudes  been  measurable  at  the  other  three  stations  (Barrett,  Riverside, 
and  Woody)  the  magnitude  would  have  been  estimated  at  4.  05  rather  than  4.  25. 
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Magnitudes  were  next  computed  for  Blanca  and  Logan  from  nine  and  six 
stations,  respectively,  at  greater  distances.  The  results  are:  for  Blanca, 

M  =  4a  8  j-  0.  4*;  and  for  Logan,  M  =4.  4  -fO.  5*.  Thus,  magnitudes  estimated 
by  the  two  methods  are  in  agreement.  Based  on  earthquake  studies  it  might 
be  expected  that  the  magnitudes  based  on  P- wave  data  from  distant  stations 
should  exceed  those  based  on  shear  wave  data  from  nearby  stations  by  about 
one-half  magnitude  unit. 

A  magnitude-yield  relationship  can  be  derived  from  the  data  so  far  pre° 
sented.  The  basic  definition  of  magnitude  is 

!  M  =  log  Ad  -  log  Bd  ,  (1) 

where  Ad  is  the  amplitude  at  distance  d  and  B.d  is  the  amplitude  of  a  zero  mag¬ 
nitude  shock  at  the  same  distance.  The  observed  relation  between  amplitude 
and  yield  (Y)  is  given  by 

A=KdY  ,  (2) 

where  is  the  amplitude  observed  at  a  distance  d  from  an  explosion  of  yield 
Y  kilotons.  From  (1)  and  (2),  it  may  be  seen  that 

M  =  C  +  log  Y  ,  (3) 

where  C  is  a  constant  [  equal  to  log  (K/B)]  which  can  be  evaluated  if  the  magni¬ 
tude  is  known  for  a  given  yield.  Estimates  for  C  can  be  made,  based  on  the 
known  yields  and  magnitudes  of  Rainier,  Logan  and  Blanca.  The  relationship 
is  then  found  to  be 

M  =  3.  64  +  log  Y  . 

This  relationship  may  not  apply  for  shots  in  other  localities  or  in  other  media. 

IV.  SUMMARY 

Seismic  waves  from  recent  underground  nuclear  explosions  in  Nevada  were 
observed  at  a  number  of  temporary  stations  along  a  line  extending  eastward  to 
Maine.  Studies  of  seismograms  from  these  stations  and  from  a  large  number 
of  permanent  stations  have  given  the  following  tentative  results: 

(a)  Between  150  and  1100-kilometers  the  times  of  P  approximately  fit  a 
line  of  slope  of  1/7.  9  sec/km. 

Standard  deviation  of  a  single  observation. 
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(b)  At  1200  km  a  strong  wave  of  higher  apparent  surface  speed  arrives 
about  four  seconds  late  with  respect  to  the  straight-line  curve. 

(c)  Between  1200  and  2500  km  the  times  and  waveforms  show  multiple 
arrivals  which  may  be  associated  with  cusps  in  the  travel-time  curve. 

(d)  Short-period  shear  waves  of  speed  near  3.  5  km/ sec  were  observed  to 
about  2000  km. 

(e)  Between  200  and  1100  km  the  amplitude  of  P  varies  inversely  as  the 
cube  of  the  distance. 

(f )  Between  200  and  2000  km  the  amplitude  of  S  varies  inversely  as  the 
cube  of  the  distance.  The  vertical,  radial,  and  transverse  components  are  of 
approximately  equal  size. 

(g)  At  large  distances,  the  amplitudes  of  body  waves  are  proportional  to 
the  first  power  of  the  explosive  yield  (Y). 

(h)  The  explosions  produced  seismic  waves  equivalent  in  size  to  those 
from  natural  earthquakes  of  magnitude: 

M  =  3.  6  +  log  Ykt  . 
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7.  TEMPERATURE  AND  RADIATION  DISTRIBUTIONS 
IN  THE  RAINIER  SHOT  ZONE 

Jo  Lo  Olsen,  W0  P0  Bennett,  and  D.  E.  Nielsen 
Lawrence  Radiation  Laboratory,  Nevada 

In  an  attempt  to  understand  better  the  phenomenology  of  underground  nuclear 
explosives  and  their  applications  to  the  Plowshare  Program,  extensive  radiation 
and  temperature  measurements  were  made  in  the  Rainier  shot  zone.  Measure¬ 
ments  are  also  being  made  at  the  present  time  in  the  shot  zones  of  Blanca, 

Logan,  Neptune  and  Tamalpais. 

Prior  to  the  Rainier  event,  some  background  temperature  measurements 
were  made  in  deep  vertical  holes  at  the  Rainier  site  by  the  USGS.  Broadview 
Research  Corporation  made  pre-Rainier  temperature  studies  and  did  early 
planning  for  the  post- shot  temperature  measurements  for  Rainier,  The  post¬ 
shot  temperature  and  radiation  measurements  on  this  event  and  similar  meas¬ 
urements  on  all  subsequent  underground  nuclear  detonations  were  conducted  by 
LRLo 

In  order  to  make  the  measurements,  it  was  necessary  to  drill  holes  through 
the  zone  of  effect  using  standard  underground  drilling  techniques  with  modifica¬ 
tions  for  remote  operations  and  radiation  control.  The  cores  taken  from  these 
holes  were  used  for  radiochemical  analyses  and  physical  properties  studies. 

The  radiation  distribution  in  the  shot  zones  was  determined  by  logging  the 
drill  holes  for  gamma  radiation  using  a  standard  Jordon  RAMS  unit  having  a 
range  of  1  mr/hr  to  1000  r/hr„  This  unit  utilizes  as  the  detector  a  thick  wall 
ionization  chamber.  For  areas  of  low  gamma  intensity,  a  standard  "Drill  Hole11 
gamma  probe  was  employed.  This  unit  uses  a  Geiger  tube  mounted  in  a  thick 
wall  probe  and  has  a  range  of  0.  01  mr/hr  to  20  mr/hr.  Both  units  were  cali¬ 
brated  with  a  Co^  source  before  and  after  each  logging.  The  radiation  values 

V»- 

recorded  are  accurate  to  +  15%. 

Temperature  distributions  were  measured  using  platinum  resistance  ther¬ 
mometers  (100  ohms  at  20°C)  and  thermistors  (2000  to  18000  ohms  at  20°C). 
These  units  were  also  calibrated  before  and  after  each  logging.  Calibration  in 

^T.  C.  Goodale  et  al.  »  "Temperatures  from  Underground  Detonation,  Shot 
Rainier,"  Operation  Plumbbob  Report  WT-1527,  July  1958. 
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this  case  was  done  in  a  Fisher  temperature  bath  using  National  Bureau  of 
Standards  mercury  thermometers  accurate  to  +  0,  05  °C,  The  overall  calibra¬ 
tions  were  reproducible  to  0,  5°C,  The  position  accuracy  for  the  probes  in 
the  holes  was  +  6  inches. 

Figure  70  1  is  an  elevation  and  Fige  7,  2  a  plan  view  of  the  Rainier  shot  zone 
showing  the  locations  of  the  drill  holes  and  an  exploratory  drift  driven  through 
the  zone. 

The  first  hole  drilled  into  the  zone  of  effect  was  a  vertical  hole  (No,  R-6) 
drilled  from  the  surface  ground  zero  to  a  depth  of  912  feet.  This  hole  was 
completed  approximately  two  months  after  the  shot.  No  radiation  above  back¬ 
ground  (0,  05  mr/hr)  was  detected  in  this  hole  to  a  depth  of  848  feet.  Measure¬ 
ments  could  not  be  made  below  848  feet  because  of  drilling  tools  lost  in  the  hole. 
Gas  and  soil  samples  were  obtained  from  the  cavity  (Fig,  7,  1)  386  feet  above  the 
shot  point.  Radiochemical  analysis  of  these  samples  showed  a  very  small  amount 
of  Kr  "  and  traces  of  Nd  and  Y.  Temperatures  recorded  in  this  hole  by  USGS 
(Fig,  7,  3)  were  lower  than  the  background  temperatures  measured  by  them  in 
an  adjacent  pre-shot  vertical  hole  (Fig0  7,4),  This  difference  can  be  explained 
by  the  68s  000  gallons  of  drilling  water  (10  to  15°C)  which  were  lost  to  the  me¬ 
dium  while  drilling. 

The  B}  CP  and  D  holes  (Fig,  7,  1)  were  completed  from  3  to  5  months  after 
the  detonation.  Figures  7,  5  through  7,  10  show  the  results  of  the  radiation  and 
temperature  measurements  made  in  these  holes  both  5  and  18  months  after  the 
shot.  In  December  1958S  the  investigation  of  the  Rainier  zones  which  had  pre¬ 
viously  been  suspended  because  of  other  activities  at  the  test  site3  was  resumed. 
Holes  GP  Hp  and  I  were  drilled  and  an  exploratory  drift  was  driven  as  shown  in 
Fig,  7,  1,  The  results  of  the  radiation  and  temperature  measurements  made  in 
these  holes  are  shown  in  Figs,  7,  11  through  7,  15,  No  radiation  above  background 
(0,  05  mr/hr)  was  observed  in  the  G  hole.  Figure  7,  16  shows  the  radiation 
levels  measured  along  the  wall  of  the  exploratory  drift. 

All  of  the  radioactivity  was  found  within  a  radius  of  68  +  9  feet  from  ground 
zero  (Fig,  7,  17)  except  for  the  small  quantities  of  Kr,  Nd  and  Y  mentioned  ear¬ 
lier. 

Figure  7,  18  is  an  isothermal  plot  of  the  temperatures  observed  in  the  B, 

Cp  and  D  holes  approximately  5  months  after  detonation,  ,  Figure  7,  19  is  an 
isothermal  plot  of  the  temperatures  observed  in  the  B,  Cp  D,  G  and  H  holes 
approximately  18  months  after  the  shot.  Temperatures  above  background  were 
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observed  to  a  radius  of  approximately  180  feet  from  ground  zero.  However, 
since  the  temperature  profiles  approach  background  asymptotically,  tempera¬ 
ture  differences  of  the  order  of  the  accuracy  of  the  measuring  equipment  cor¬ 
responded  to  large  changes  in  the  volume.  Therefore,  it  was  necessary  to 
select  a  well-defined  point  for  thermal  energy  calculations.  The  knee  of  the 
temperature  curve,  that  is,  the  point  at  which  the  temperature  begins  to  rise 
sharply  above  background,  was  selected. 

On  Rainier,  the  knees  of  the  curves  from  the  B,  C,  and  D  holes  as  meas 
ured  approximately  5  months  after  the  shot  all  correspond  to  a  temperature  of 
20°C  and  a  distance  of  96  4-  1  ft  from  ground  zero.  A  rough  integration  of  the 
temperature  distribution  given  in  Fig.  7.  18  shows  that  the  heat  added  to  the 
medium  within  a  radius  of  96  feet  is  equal  to  8.4  X  10^  calories  or  approxi¬ 
mately  1/2  the  total  energy  of  the  device  (1.  7  X  10^  calories).  An  integration 
over  the  same  volume  using  the  temperature  distribution  shown  in  Fig.  7.  19  as 
measured  18  months  after  detonation  shows  that  the  heat  remaining  is  7.  2  X  10 
calories,  a  decrease  of  14%  in  13.  5  months.  In  this  same  period,  the  zone 
heated  above  20 °C  increased  in  volume  by  57%. 
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8.  THE  DISPOSITION  OF  RADIOACTIVITY 
R.  H.  Goeckermann 

Lawrence  Radiation  Laboratory,  Livermore 

The  material  of  this  paper  is  based  on  radiochemical  analyses  by  the  forty 
members  of  the  Radiochemistry  group  at  LRL,  Livermore,  of  samples  obtained 
from  underground  explosions  of  0.  07  to  19  koloton  yield  range  in  tuff  at  the  NTS 
during  1957-1958.  The  sampling  included  rock  cores  and  assorted  samples  of 
the  vented  material  from  totally-  or  partially-contained  explosions.  The  radio¬ 
nuclides  studied  included  tracer  activities,  fission  products,  and  induced  ac¬ 
tivities  of  observable  half-life  —  weeks  to  years. 

It  is  necessary  to  review  briefly  the  early  history  of  an  underground  nuclear 
explosion  in  tuff.  In  explosions  of  a  few  kilotons,  the  maximum  bubble  radius 

1  /3 

is  reached  in  5-10  milliseconds.  This  time  scales  roughly  like  W  .  The  vapor 
phase  in  tuff  is  at  high  temperature  (1200-1  500° C)  and  pressure  (~40  atmos¬ 
pheres)  and  consists  largely  of  water  vapor.  The  viscous  fused  tuff  acts  to 
seal  the  cavity.  Collapse  begins  at  times  ranging  from  immediately  to  tens  of 
seconds  later  depending  on  the  nature  of  the  medium  and  degree  of  containment. 

If  the  cavity  persists,  flow  and  dripping  of  the  fused  material  will  Occur.  With 
collapse  comes  a  sudden  release  of  pressure  and  the  escape  of  the  vapor  phase 
from  the  immediate  vicinity  of  the  cavity. 

The  tracer  species  that  have  been  studied  in  several  explosions  are  plu¬ 
tonium,  uranium,  the  actinide  earths,  protactinium,  arsenic,  cesium,  hafnium, 
and  tantalum.  Almost  all  of  the  plutonium  has  been  found  in  the  glass  resulting 
from  solidification  of  the  fused  tuff.  The  concentration  indicates  200-500  tons  of 
fused  rock  mixed  with  debris  per  kiloton  yield  results  from  explosions  in  the 
1  to  20  kiloton  range.  Spongy-looking  material  containing  a  factor  of  20  lower 
concentration  of  Pu  has  been  found  in  fissures  in  the  rock  at  distances  of  four 
cavity  radii.  The  disposition  of  uranium  is  such  as  to  indicate  its  behavior  is 
considerably  more  volatile  than  plutonium  as  a  reference  species.  Glassy  shell 
samples  have  been  found  which  have  only  one-third  the  expected  uranium  content 
and  samples  away  from  the  cavity  such  as  the  fissure  material  are  enriched 
in  uranium  by  a  factor  of  approximately  2.  When  collapse  occurs  immediately, 
as  much  as  half  of  the  uranium  possibly  escapes.  The  actinide  earths  behave 
more  like  plutonium,  as  do  protactinium,  hafnium,  and  tantalum.  Arsenic 
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behaves  like  a  very  volatile  species  in  these  circumstances.,  The  majority  of 
the  fused  rock  retains  less  than  1%  of  the  arsenic,  although  isolated  droplets 
have  been  found  which  have  a  normal  concentration,,  When  venting  to  the  ground 
surface  occurred,  the  vented  debris  was  enriched  in  arsenic  by  an  order  of 
magnitude.  Studies  with  cesium  tracer  have  shown  that  approximately  half  of 
the  cesium  is  retained  in  the  glass  while  vented  material  is  enriched  by  a  factor 
of  2  in  cesium. 

The  detonation  of  a  nuclear  explosive  produces  fission  products  in  propor¬ 
tion  to  the  fission  yield  of  the  explosive.  Twenty  kilotons  of  fission  energy  pro¬ 
duces  approximately  1  kg  of  fission  products,  including  20  g  of  strontium-90, 
which  is  comparable  to  the  waste  resulting  from  one  day’s  operation  of  a  1000- 
Mw  reactor.  The  distribution  of  fission  product  activity  follows  the  familiar 
fission  yield  curve.  The  behavior  of  the  observable  species  can  be  understood 
only  after  examination  of  the  decay  chains  leading  to  them  (these  chains  are 
listed  in  Table  I). 

Table  L  Fission  Product  Decay  Chains, 


A/T  nr  A  0,43  40 

Mass  85:  As - ^Se  — 


s 

15 


■Br..  3‘  °>Kr 


-S>  Sr 


89:  Kr-lli — ^Rb- 

m  m 


33  2  7 

90:  Kr  _ii - >  Rb - 111  >  Sr 


m 


91:  Kr 
95:  Rb 
99:  Zr 


9.  8 


-5>Rb 


2  and  14 


sh 


25% 

30 


Sr 
^  Nb 


m 

sh 


->  Y 


2.  5 


137:  I  —Bl - 3*Xe 


m 
3.  8 


Mo 


140:  Xe 


16 


->  Cs 


m 

66 


Cs 


141:  Xe  -  7  >  Cs 


s 

sh 


->  Ba 


m 


W  Sr. 


9.  7^ 
IT 

10 


m 


>Zr 


„  .  18  .  3.  7  _ 

-$►  Ba - >La — r^-Ce 

m  h 
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The  history  of  the  products  of  these  chains  is  given  below: 


Krypton- 8  5 

Strontium- 89 

Strontium- 90 

Yttrium-91 

Zirconium- 9  5 
Molybdenum- 99 
Cesium- 137 

Barium-140 

Cerium- 141 

Cerium-144 
Neodymium- 147 


less  than  1%  appears  to  be  trapped  in  the  glass  in  a  2-kt 

9 

explosion,  the  majority  was  diluted  into  ~10  cu  ft  of  air 
in  the  chimney. 

3-10%  in  the  major  part  of  the  glass,  enriched  in  the  rub¬ 
ble  by  several-fold  and  in  the  vented  debris  by  nearly  an 
order  of  magnitude. 

20-30%  in  the  glass,  vented  debris  enriched  by  5-fold, 
local  fallout  less  enriched. 

30-50%  in  the  glass,  vented  debris  and  rubble  enriched 
several-fold. 

.70-100%  in  the  glass,  rubble  enriched, 
no  data  on  glass,  vented  debris  enriched  by  30-60%, 
20-40%  in  the  glass,  vented  material  enriched  by  several¬ 
fold. 

30-60%  in  glass  decreasing  as  yield  increases,  vented 
debris  enriched  several-fold. 

approximately  60%  in  the  glass,  escaping  debris  enriched 
few-fold. 

~1 00%  in  the  glass,  vented  debris  has  same  proportion  as 
plutonium. 

precursors  are  rare  earths,  behavior  like  plutonium. 


The  species  calcium-45,  produced  by  neutron  capture  in  the  small  per¬ 
centage  of  calcium  in  tuff,  has  been  measured  in  several  samples.  The  activa- 

-4  45 

tion  has  varied  from  (1  to  20)  X  10  atoms  Ca  per  fission,  depending  on  the 

amount  of  neutron  absorbing  material  around  the  nuclear  explosive.  In  a  given 

45 

explosion,  the  ratio  of  Ca  to  plutonium  varies  by  a  factor  of  2  in  different 

,.45 

glass  samples,  Isolated  debris,  such  as  the  fissure  sample,  is  low  in  Ca 

In  conclusion,  the  disposition  of  the  radioactivity  is  a  function  of  the  yield 
as  it  affects  the  time-temperature  history,  the  degree  of  containment  and  as¬ 
sociated  collapse  phenomena,  and  certainly  the  nature  of  the  medium  in  which 
the  explosion  takes  place,  Subsequent  measurements  will  be  made  in  explo¬ 
sions  which  approximate  more  closely  typical  Plowshare  applications. 
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9.  THE  LEACHING  OF  RAINIER  DEBRIS 
W.  D.  Bond  and  W.  E.  Clark 
Oak  Ridge  National  Laboratory 

The  possibility  of  serious  ground  water  contamination  is  an  important 
consideration  in  the  peaceful  application  of  contained  nuclear  explosions.  A 
series  of  leaching  tests  were  performed  on  Rainier  debris  to  establish  the 
extent  of  ground  water  contamination  which  would  result  from  such  debris. 
Leaching  conditions  were  chosen  so  that  the  results  would  represent  the  max¬ 
imum  extractability  of  radioactivity. 

Samples  of  the  most  radioactive  Rainier  debris  available  were  crushed 
and  sieved  to  divide  the  particles  into  five  size  classifications  dfefined  by 
diameters  as  follows: 

(1)  less  than  53  microns  (<  270  mesh) 

(2)  53  to  74  microns  (200-270  mesh) 

(3)  74  to  149  microns  (100-200  mesh) 

(4)  149  to  297  microns  (48-100  mesh) 

(5)  297  to  1 190  microns  (14-48  mesh). 

Specific  activities  were  determined  for:  the  size  ranges  (1),  (2),  and  (4),  re¬ 
spectively,  and  one-gram  samples  of  these  size  ranges  were  digested  over¬ 
night  with  100  ml  of  the  various  reagents  of  interest  as  indicated  in  Tables 
Ia-Ic.  Solids  were  then  centrifuged  off,  the  activity  of  the  supernatant  liquid 
measured,  and  the  percentage  of  each  type  of  activity  extracted  was  calcu¬ 
lated. 

The  results  are  listed  in  Tables  Ia-Ic.  Though  there  is  considerable 
scatter  in  the  data  it  is  obvious  that  extraction  of  any  material  represented 
by  the  measured  activities  would  be  completely  impractical  since  the  highest 
recoveries  obtained  are  about  27%.  The  addition  of  MnC^  to  give  the  leach 
oxidizing  power  shows  no  clear  trend  so  far  as  extraction  efficiency  is  con¬ 
cerned.  It  is  possible  that  any  increase  in  extraction  due  to  oxidation  is  off¬ 
set  by  adsorption  on  the  suspended  particles  of  MnC^. 

The  extraction  of  activity  by  ground  water  from  the  Nevada  Test  Site 
(Hagestad  Hole)  is  of  interest  because  it  gives  a  measure  of  the  amount-of 
activity  which  can  be  expected  to  be  leached  from  the  residue  after  a  test  shot. 
The  room  temperature  data  are  presented  in  Table  II  and  Fig.  9.  1.  Since 
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Table  la.  Leaching  Experiments  on  Rhyolite  Debris  from  Rainier. 
Conditions:  100  ml  of  solution,  1  g  of  debris,  temperature 


boiling.  <  53  (i  diam  fraction  of  the  debris. 
No  stirring. _ 


Leaching 

agent 

MnC>2 

Leach  time 

Activity  extracted  (%) 

(hr) 

Alpha3" 

Beta 

Gamma 

H  SO. 

2  4 

0.  5  M 

Excess 

15 

9.  12 

0.  85 

0.  28 

0.  5  M 

Exc  ess 

16 

4.  84 

0.  98 

0.  29 

0.  5  M 

None 

16.  5 

20.  2 

2.  02 

0.  90 

3.  0  M 

Excess 

16 

4.  64 

0.  81 

0.43 

3.  0  M 

None 

16.  5 

26.  8 

2.  00 

1. 26 

hno3 

6  M 

None 

16 

2.  69 

1.  74 

0.  82 

HC1 

6  M 

Excess 

16.  5 

6.  87 

2.  99 

1. 29 

6  M 

Excess 

16 

6.  26 

2.  65 

1. 30 

6  M 

None 

16.  5 

16.4 

1. 29 

0.  69 

NaOH 

1  M 

None 

16.5 

20.  1 

1. 22 

3.  20 

Na2CO  ,  NaHCO 

0.5M-0.5M 

Excess 

16.  5 

12.  2 

0.  28 

0.  10 

0.  5  M  -  0.  5  M 

Excess 

16 

6.  34 

0.  12 

0.  03 

0.  5  M  -  0.  5  M 

None 

16.  5 

14.  5 

0.  66 

0.  28 

Distilled  H^O 

None 

21. 5 

0 

0.  31 

0.  28 

Distilled  H^O 

None 

21.  5 

0 

0.  32 

0.  38 

Ground  H^O^ 

None 

16.  5 

3.  21 

0.  55 

0.41 

cl 

Includes  all  alpha  activities  except  uranium. 


From  the  Hagestad  Hole  on  the  NTS.  pH  =  10.  2. 
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Table  lb.  Leaching  Experiments  on  Rhyolite  Debris  from  Rainier. 

Conditions:  149-297  (i  diam  fraction  of  the  debris. 

_ Solutions  stirred. _ 

,  .  _  .  »  Activity  extracted  (%) 

Leaching  .  .  n  Leach  time  - - - 

agent _ MnU2 _ (hr) _ Alpha3- _ Beta _ Gamma 


H2S°4 


0.  5  M 

Excess 

16.  5 

2.47 

0.  38 

0.  25 

0.  5  M 

None 

16.  5 

2.  01 

1. 08 

0.  70 

3.  0  M 

Exc  ess 

16.  5 

2.48 

0.  56 

0.45 

HC1 

6  M 

Excess 

16.  5 

1.  50 

1.  16 

0.  65 

6  M 

None 

16.  5 

1.  32 

0.  89 

0.  71 

Na2C03,  NaHC03 

0.  5  M  -  0.  5  M 

Excess 

16.  5 

2.49 

0.  092 

0.  044 

0.  5  M  -  0.  5  M 

None 

16.  5 

1. 29 

0.  31 

0.  35 

Ground 

None 

16.  5 

1.44 

0.  52 

0.44 

Ground 

Exc  ess 

16.  5 

0.  69 

0.  30 

0.  17 

hno3 

6  M 

None 

16.  5 

2.  32 

1.  18 

0.  73 

aIncludes  all  alpha  activitie 

s  except  uranium. 

From  the  Hage 

stad  Hole  on  the  NTS. 

pH  =  10.  2 

Table  Ic.  Leaching  Experiments  on 

Rhyolite  Debri 

s  from  Rainier. 

Conditions: 

297-1190 

|i  diam  fraction  of  the  debris. 

Solutions 

stirred. 

Leaching 

lVT-nO 

_  .  . .  Activity  extracted  (%) 

Leach  time  - - - 

agent 

iviiiv-/  2 

(hr) 

Alpha 

Beta 

Gamma 

H2S°4 

0.  5  M 

Excess 

16.  5 

1. 64 

0.  19 

0.  01 

0.  5  M 

None 

16.  5 

1. 31 

1. 27 

0.  75 

HC1 

6  M 

Exc  ess 

16.  5 

1. 05 

0.  56 

0.  035 

6  M 

None 

16.  5 

5..  61 

0.  68 

0.40 

hno3 

6  M 

None 

16.5 

1. 87 

0.  55 

0.  029 

Na?COv  NaHC03 

0.  5  M  -  0.  5  M 

Exc  ess 

16.  5 

1.  31 

0 

0.  005 

11 

includes  all  alpha  activities  except  uranium. 
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Table  II.  Leaching  of  Rainier  Debris  by  Ground  Water  at  Room 

T  emperature  • 

1-g  samples  ground  to  <  53  p  diameter,  stirred  magnetically 
with  100  ml  H^O  (Hagestad  Hole). 

5-ml  specimens  removed  periodically  for  analysis. 

Activity  of  solid  specimen: 

Gross  alpha  (except  U)  =150  counts/min/mg 
Beta  =  1. 82  X  103  counts/min/mg 

Gamma  (scintillation)  =  9.  05  X  103  counts/min/mg 


Specimen 

Total  leach 

Percent  activity  extracted 

No. 

time  (hr) 

Alpha 

Beta 

Gamma 

1 

2 

0.  20 

0.  11 

0.  065 

1 

4 

0.  20 

0.  14 

0.  050 

1 

6 

0.  20 

0.  20 

0.  11 

1 

8 

0.  39 

0.  26 

0.  12 

2 

16 

0.  20 

0.  29 

0.  15 

1 

24 

0.  39 

0.  40 

0.  28 

2 

24 

0.  20 

0.  34 

0.  18 

1 

32 

0.  39 

0.  53 

0.  39 

2 

40 

0.  59 

0.  48 

0.  34 

2 

48 

0.  52 

0.  54 

0.  41 

Includes  all  alpha  activities  except  uranium. 


these  tests  were  run  on  the  finest  size  classification,  under  stirred  conditions 
and  with  a  100/1  weight  ratio  of  water  to  specimen,  it  is  believed  that  this 
should  represent  the  most  favorable  condition  for  removal  of  activity  from 
the  rock.  Since  the  maximum  activity  leached  is  less  than  0.  6%  it  appears 
that  there  will  be  relatively  little  contamination  of  the  ground  water.  It 
should,  however,  be  pointed  out  that  at  the  end  of  the  48 -hr  contact  period  the 
rate  of  leaching  appeared  to  be  relatively  constant  at  about  0.  05%/hr;  it  is, 
therefore,  possible  that  over  longer  residence  times  very  sizeable  activities 
could  be  built  up.  Long  residence  time,  however,  implies  a  slow  rate  of 
movement  for  the  water  through  the  ground  water  table  with  a  correspondingly 
good  chance  of  removal  of  activity  by  the  inert  medium  traversed.  It  should 
be  noted  (Table  la)  that  less  activity  is  removed  from  the  rhyolite  by  distilled 
water  than  by  the  ground  water  which  had  a  pH  of  approximately  10.  2. 
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C.  The  Theory  of  Underground 
Nuclear  Explosions 


1.  CAVITY  FORMATION 
G.  T.  Pelsor 

Lawrence  Radiation  Laboratory,  Livermore 

In  Fig.  1.  1  the  solid  curve  approximates  the  Hugoniot  for  water  satu¬ 
rated  tuff  which  surrounded  the  Rainier  nuclear  detonation.  A  shock  of  1.1  Mb 
pressure  would  increase  the  specific  internal  energy  of  this  material  by  an 
amount  equal  to  the  area  of  the  triangle  below  the  straight  line  joining  the  ini¬ 
tial  and  final  states.  As  the  pressure  of  this  same  material  falls,  the  material 
expands  nearly  adiabatically.  Using  the  Hugoniot  as  an  approximation  to  the 
adiabat,  the  hatched  areq  represents  the  net  gain  in  internal  energy  when  the 
pressure  has  fallen  to,  say,  100  bars.  The  3200  cal/g  was  assumed  to  be 
sufficient  to  vaporize  the  rock  under  this  pressure. 

A  one -dimensional  hydrodynamic  calculation  was  made  on  a  rather  crude 
machine  code  (as  few  as  4  zones  were  used  occasionally)  which  did,  however, 
provide  for  two  different  behaviors  of  the  material  after  passage  of  the  shock. 
That  is,  if  the  shock  pressure  exceeded  1.  1  Mb,  the  expansion  was  that  of  an 
ideal  gas  of  y  =  1.4;  on  the  other  hand,  material  subjected  to  weaker  shocks 
than  this  was  allowed  to  expand  along  the  Hugoniot.  Figure  1.  2  shows  two 
pictures  taken  from  this  calculation.  In  the  one  at  the  top,  the  shock  pressure 
has  just  gotten  down  to  1.  1  Mb.  Material  which  was  on  the  wall  of  the  cavity 
initially  at  1.  2  m  has  now  moved  out  to  about  2  m  and  has  a  larger  specific 
volume  than  it  had  initially.  The  material  now  at  the  front  at  2.  65  m  will 
move,  as  shown  in  the  later  picture  below,  to  a  point  where  there  is  a  jog  in 
the  specific -volume  curve,  indicating  the  difference  in  behavior  of  the  two 
phases.  Many  milliseconds  later,  the  material  outside  this  jog  forms  the 
liquid  wall  of  the  cavity  at  a  radius  of  about  55  ft. 

But  the  situation  in  the  liquified  rock  is  complicated  by  the  fact  that  one 
component,  namely,  water,  is  much  more  volatile  than  the  rest.  Most  of  this 
water  will  eventually  be  found  as  part  of  the  cavity  vapor.  Its  emergence  is 
slowed  by  the  viscosity  and  surface  tension  of  the  liquid  rock.  In  fact,  the 
water ’in  a  thin  outer  layer  of  the  fused  rock  is  trapped  in  the  viscous  glass. 

It  may  be  that  most  of  the  water  vapor  from  the  molten  rock  gets  into  the  cavity 
so  late  in  the  process  that  it  takes  little  part  in  determining  the  cavity  size. 
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When  this  water  does  arrive,  having  started  from  a  region  of  lower  tempera¬ 
ture  than  the  cavity  vapor  and  having  cooled  still  further  by  expansion,  it  chills 
the  rock  vapor  and  leads  to  some  condensation  as  has  been  discussed  by  Kennedy 
and  Higgins .  They  deduce  a  water  vapor  pressure  of  40  bars  during  the  growth 
and  freezing  of  the  small  beads  which  have  been  recovered  in  drilling  through 
the  Rainier  site  after  the  shot.  The  total  pressure  in  the  cavity  including  that 
exerted  by  the  rock  vapor  was  probably  considerably  higher,  perhaps  even 
higher  than  the  overburden  pressure,  since  the  thick  cushion  of  crushed  rock 
between  55  ft  and  130  ft  would  allow  a  higher  equilibrium  pressure  in  the  cavity 
than  the  radial  stress  at  the  outer  edge  of  the  crushed  zone0  The  fact  that  the 
roof  took  a  minute  or  two  to  collapse  is  further  indication  of  cavity  pressure 
higher  than  overburden  pressure. 

Radioactive  products  found  beyond  the  cavity  wall  lead  one  to  suspect  the 
formation  of  fissures  (see  Fig.  10  3)c  Although  the  location  was  such  that  fis¬ 
sure  propagation  could  have  been  influenced  by  proximity  of  the  tunnel  or  a 
weak  faulted  zone  in  the  rock,  there  may  have  been  other  fissures  which  have 
been  missed  in  the  limited  drilling  and  in  the  exploratory  drift.  One  might 
imagine  these  fissures  to  be  formed  somewhat  as  follows:  In  the  early  stages 
of  the  explosion,  the  material  just  outside  the  liquid  is  deformed  plastically 
with  radial  compression  and  lateral  spreading  as  it  moves  to  larger  radii. 
Eventually  as  the  pressure  decreases  this  material  is  able  to  resist  further 
distortion.  This  point  may  be  reached  while  it  is  still  moving  outward  so  that, 
having  low  tensile  strength,  it  tea.rs  along  radial  cracksQ  These  cracks  grow 
rather  slowly  at  first  because  they  are  moving  into  a  region  where  the  pressure 
is  higher  than  the  cavity  pressure.  But  if  the  cavity  pressure  is  higher  than 
the  overburden  pressure,  as  the  pressure  wave  goes  out  and  so  becomes  weaker, 
the  pressure  gradient  is  soon  reversed.  Then  the  cracks  grow  longer  rapidly 
because  of  stress  concentration  at  the  end  of  the  crack.  This  gives  a  large  area 
of  cold  rock  on  which  vapor  can  condense  and  leads  to  the  sudden  fall  in  pressure 
inferred  by  Kennedy  and  Higgins.  The  cracks  would  form  preferentially  in  the 
roof  of  the  cavity  where  overburden  pressure  is  least  and  the  unstable,  water- 
saturated,  crushed  rock,  no  longer  supported  by  cavity  pressure,  would  fall 
into  the  cavity  in  a  pile  of  rubble  (Fig.  1.  3).  The  displaced  vapor  would  then  be 
even  more  intimately  mixed  with  cold  rock  and  would  further  condense.  The 
crushed  rock  beyond  the  cracks,  and  even  the  unaltered  rock  above,  left  without 
sup^~t,  would  proceed  to  collapse  —  perhaps  in  a  succession  of  falls  which  in  a 
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couple  of  months  had  reached  a  height  of  386  ft  above  the  detonation  point.  The 
condensed  water  percolated  back  down  to  the  hot  region  where  it  evaporated  and 
started  another  cycle  of  diffusion  upward,  condensation  in  the  cold  chimney,  and 
so  forth.  This  gradually  decreased  the  maximum  temperature  in  the  system  to 
the  boiling  point  of  water,  after  which  the  chief  mechanism  of  cooling  is  the  con¬ 
duction  of  heat  through  the  surrounding  rock. 

The  compressed  rock,  in  place,  is  relatively  impermeable  while  that  which 
was  broken  in  the  cave-in  is  open  and  permeable.  In  a  recent  hole  inclined  up¬ 
ward  through  the  Rainier  chimney,  circulation  of  drilling  mud  was  lost  at  a  point 
about  34  ft  horizontally  outward  and  74  ft  above  the  similar  point  on  the  horizontal 
hole  through  the  detonation  point.  Thus  the  caving  region  flares  out  about  25° 
from  the  vertical. 

If  the  process  description  we  have  given  is  correct,  it  depends  at  a  number 
of  points  and  in  a  rather  sensitive  way  on  the  properties  of  the  medium  in  which 
the  detonation  takes  place.  For  example,  the  strength  of  the  medium  for  dynamic, 
one-dimensional  loading  may  be  an  important  factor  in  the  formation  of  fissures 
and  thus  in  the  final  shape  and  state  of  the  cavity.  We  hope  eventually  to  have 
data  on  cavities  formed  in  media  other  than  tuff. 
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2„  A  COMPUTER  CALCULATION  OF  RAINIER 
(The  First  100  Milliseconds) 

John  H0  Nuckolls 

Lawrence  Radiation  Laboratory,  Livermore 

The  purpose  of  these  calculations  is  to  test  and  develop  our  capability 
for  predicting  the  phenomena  of  underground  nuclear  explosions  from  the  laws 
of  physics  and  the  properties  of  materials,,  We  also  seek  to  understand  the 
experimental  observations  and  measurements.  Analytic  methods  are  extremely 
difficult  because  of  the  nonlinearity  of  the  partial  differential  equations,  the 
complexity  of  the  stress -strain  relations,  and  the  transitions  between  gaseous, 
liquid,  plastic,  fractured,  and  elastic  states „ 

A  special  code,  UNEC  (Underground  Nuclear  Explosion  Code),  was  de¬ 
veloped  for  the  IBM  709  Electronic  Computer  to  make  these  calculations,,  They 
extend  in  time  from  a  few  microseconds  to  approximately  100  milliseconds 
after  the  Rainier  explosion,  and  in  space  to  a  radius  of  several  hundred  feet 
from  the  site  of  detonation,,  Results  are  generally  in  good  agreement  with 
experiment,, 


METHOD,  ASSUMPTIONS 

Technically,  UNEC  is  a  spherically  symmetric  Lagrangian  hydro-elasto- 
dynamic  code.  The  basic  difference  equations  are  derived  from  the  partial 
differential  equations  for  conservation  of  momentum  and  energy: 


3r 


(t  —  q)  +  •—  (t  - 
rr  ^  r  rr 


w  = 


oo 

pr 


Conservation  of  momentum 


elasto-dynamic  I 


Conservation  of 
energy 


hydrodynamic  ) 
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where  r  =  radius 

S  =  radial  displacement 
V  =  volume 
p  =  density 

e  =  internal  energy  per  initial  volume 
P  =  pressure 

Trr*  T00  =  radial  and  tangential  components,  in  spherical  coordinates, 

of  the  stress  tensor 

q  =  viscosity. 

The  system  of  integration  automatically  conserves  mass. 

Outside  the  elastodynamic  region  (where  the  elastic  limits  of  the  tuff  have 
been  exceeded),  a  generalized  Richtmyer-von  Neumann  artificial  viscosity  is 
used.  It  spreads  the  shock  front  over  several  zones,  and  increases  entropy  as 
required  by  the  second  law  of  thermodynamics.  In  the  elastodynamic  region,  the 
viscosity  is  proportional  to  the  rate  of  strain  and  produces  Kelvin-Voigt  type 
damping  of  small  amplitude  stress  waves.  ^ 

The  tuff  is  represented  by  the  generalized  Hooke  !s  law  when  elastic,  by  a 
bulk  modulus  type  equation  when  fractured,  by  the  initial  shock  Hugoniot  in  the 
plastic  and  liquid  states,  and  as  a  Thomas-F ermi-Dirac  gas  when  vaporized: 


Trr  =  (X  +  2„)|f +  2  X.  § 


3r 


Tee  =  T^  =  x|f  +  2(\  +  n)§ 


elastic  (generalized  Hooke's  law) 


P  =  -k 


AY 


fractured 


P  =  F  (V) 


Hugoniot 


crushed,  plastic 


P  =  P(V,e) 


TFD’ 


gas 


J.  von  Neumann  and  R.  D.  Richtmyer,  "A  Method  for  the  Numerical  Calculation 
of  Hydrodynamic  Shocks",  J.  Appl.  Phys.  _21,  232,  (1950). 

L.  Knopoff  and  G.  J.  F.  MacDonald,  "Attenuation  of  Small  Amplitude  Stress 
Waves  in  Solids",  Revs.  Modern  Phys.  _30,  1178  (1958). 
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where  \,  u  =  Lame  constants 

k  =  bulk  modulus. 

In  "vertical"  calculations,  a  uniform  radial  gravitational  field  is  superimposed 
on  the  system. 

MODEL 

Initially,  the  energy  of  the  nuclear  explosion  is  distributed  uniformly  as 
the  internal  energy  of  a  gas  contained  in  a  sphere  of  the  same  volume  as  the 
original  room  (rectangular,  6x6x7  feet).  It  is  assumed  that  the  initial  den¬ 
sity  of  this  gas  is  approximately  equal  to  the  mass  of  material  vaporized  during 
the  first  few  microseconds  divided  by  the  volume  of  the  room. 

It  is  assumed  that  the  tuff  has  negligible  large-scale  tensile  strength.  As 
long  as  the  tension  does  not  exceed  the  compressive  stress  due  to  the  overbur¬ 
den,  and  the  shear  strength  is  not  exceeded,  the  tuff  is  treated  as  a  linear  elas- 

3 

tic  solid.  The  elastic  constants  were  determined  by  static  laboratory  tests  to  be: 

IP  P 

bulk  modulus,  k  =  0„  036  megabars  (Mb)  (1Mb  =  1.01  X  10  dynes/cm  ) 

shear  modulus,  |jt  =  0.  055  Mb 

sound  speed,  c  =  0.  Z41  cm/fji-sec. 

After  the  tension  in  a  spherical  shell  exceeds  the  lithostatic  stress  the  shear 
modulus  there  is  set  to  zero,  (i.  e.  ,  the  diagonal  components  of  the  stress  tensor 
are  set  equal  to  a  pressure  which  is  related  to  the  volume  by  a  bulk  modulus  type 
equation  of  state). 

When  the  shear  stress  exceeds  the  dynamic  confined  shear  strength,  it  is 
assumed  that  the  tuff  flows  plastically.  This  plastic  material  resembles  mud  and 
is  relatively  impermeable  to  water*  Laboratory  tests  indicate  the  static  unconfined 

4 

shear  strength  of  tuff  is  175-350  bars.  The  dynamic  strength  is  not  known.  The 

assumption  that  it  is  twice  the  maximum  static  strength  (i0e0,  0.7  kilobars)  leads 

to  a  crushing  radius  of  130  ft  and  is  consistent  with  experimental  ratios  of  dynamic 

5 

to  static  strengths  of  concrete.  The  plastic  tuff  is  assumed  to  obey  an  equation 
of  state  identical  with  the  initial  shock  Hugoniot.^  At  pressures  less  than  10  kilobars, 

3 

G0  W.  Johnson  and  C.  Ec  Violet,  "Phenomenology  of  Contained  Nuclear  Explosions,  M 
UCRL-5124,  Rev.  I. 

4 

G.  W.  Johnson  and  C*  E.  Violet,  op.  cit. 

5 

Da  Watstein,  J.  Am.  Concrete  Inst.,  April  1953,  pp.  7Z9-440 

^Measured  experimentally  at  LRL  with  high  explosives*  Method  is  described  by 
J0  M*  Walsh  and  Re  H.  Christian,  Phys,Rev0  97,  1544  ( 1 95 5 )0 
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the  slope  of  the  Hugoniot  is  faired  into  static  measurements  of  the  bulk  modulus. 
Although  use  of  the  Hugoniot  equation  produces  the  correct  conditions  at  the 
shock  front,  it  results  in  the  deposition  of  somewhat  too  much  waste  heat  during 
the  subsequent  expansion. 

When  the  shock  is  strong  enough  to  vaporize  the  tuff  upon  rarefaction,  an 

irreversible  transition  is  made  to  a  gaseous  equation  of  state  calculated  by  the 

7 

Thomas-Fermi-Dirac  theory  and  normalized  to  the  Hugoniot  at  the  transition 

point.  It  is  estimated  that  the  total  heat  of  fusion  is  715  cal/g  and  the  total  heat 

8 

of  vaporization  is  3200  cal/g.  These  energies  will  be  deposited  by  shocks  with 
peak  pressures  of  0.  4  and  1.0  Mb,  respectively. 

Summarizing,  the  calculations  with  this  model  will  predict  a  vaporized 
region  determined  by  a  1.  0-Mb  peak  shock  pressure,  a  melted  region  correspond¬ 
ing  to  0.  4  Mb,  a  plastic  region  where  the  shear  strength  is  exceeded,  a  fractured 
region  where  the  tension  exceeds  the  overburden,  and  beyond  this  an  elastic  region. 


RESULTS 

Shock 

The  calculated  shock  arrival  times  are  shown  in  Fig.  2.1.  At  approximately 

60  feet  the  shock  velocity  becomes  subsonic.  Thereafter,  the  curve  is  the  arrival 

time  of  the  dilatational  elastic  wave.  The  results  of  Dickens  and  Porzel’s  meas- 
9 

urements  are  represented  by  circles.  Disagreement  at  small  radii  is  probably 
mainly  due  to  the  initial  nonspherical  geometry.  At  larger  radii,  the  shock 
quickly  approaches  sound  speed  and  the  arrival  time  becomes  relatively  insensitive 
to  yield. 

Figure  2.  2  shows  the  peak  shock  pressure  as  a  function  of  distance.  To  10 

-2.  35 

meters,  it  decreases  as  r  .  Tuff  is  vaporized  to  a  radius  of  7  feet  (peak 

pressure  1.  0  Mb)  and  melted  by  the  shock  to  10  feet  (0.  4  Mb).  Enough  energy  is 
deposited  within  a  radius  of  14  ft  to  melt  all  the  tuff  to  this  distance,  or  660  metric 
tons.  Radiochemical  analysis  indicates  that  600  ±  250  tons  were  melted.  ^  At 
130  ft  the  peak  pressure  is  0.7  kilobars,  twice  the  static  shear  strength. 


7 

See  R.  D.  Cowan  and  J.  Ashkin,  Phys.  Rev.  105,  144  (1957). 

g  - 

G.  T.  Pelsor,  private  communication. 

9 

R.  G.  Dickens  and  F.  B.  Porzel,  uClose-In  Time  of  Arrival  Measurements  for 


Yield  of  Underground  Shot,  M  ARF-D-133  (PR- 13). 


R.  H.  Goeckermann, 


private  communication. 


RADIUS  (METERS) 


-124- 


UCRL-5675 


TIME  OF  ARRIVAL  (MILLISECONDS)  MUL-10255 


Fig. 2.  1. 


PEAK  PRESSURE  (KILOBARS) 
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RADIUS  (METERS)  10 


Fig.  2. 2. 


70 

MJL-10256 
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Figure  2.  3  (pressure  vs  time  at  60  feet)  is  a  typical  supersonic  shock.  A 
subsonic  pressure  pulse  is  fundamentally  different  as  shown  in  Fig.  2.4  (pres¬ 
sure  vs  time  at  130  ft).  The  elastic  precursor  travels  at  dilatational  sonic  speed, 

1  /z 

[(X.  +  2 (jl) / p3  ,  and  the  main  pressure  pulse  moves  approximately  40%  as  fast 

between  130  and  200  ft. 

As  it  propagates,  the  shock  becomes  broader.  Figure  2.5  shows  the  times 

for  the  pressure  to  decay  to  1  /Z  and  l/lO  peak  values  as  a  function  of  distance. 

For  example,  at  60  ft  (peak  pressure  9  kilobars)  the  pressure  is  greater  than 

4.5  kilobars  for  ~  1.5  msec,  and  greater  than  0.  9  kilobar  for  ~12  msec.  To 

- 1 . 65 

about  20  meters,  the  lines  in  Fig.  2.5  are  approximately  parallel  to  r 

CAVITY 

Figure  2.6  shows  the  growth  of  the  central  cavity.  The  maximum  radius 
is  about  63  ft,  compared  to  62  ft  determined  by  post-shot  investigation.^  The 
average  cavity  pressure  as  a  function  of  time  is  plotted  in  Fig.  2.7.  It  decreases 
from  4.5  megabars  to  1  kilobar  in  less  than  10  msec,  and  to  about  100  bars  in  35 
msec. 


ENERGY 


At  90  msec,  less  than  10%  of  the  energy  is  kinetic.  The  total  energy  is 
distributed  as  follows: 


State 

Gas 
Liquid 
Crushed 
F  ractured 
Elastic 


Radii 

0-62' 
62'-62'3" 
62 • 3"- 1 30 1 
1 30 '  —  2  8  0  * 
>280' 


%  of  total  energy 

8.  2 
19.  1 
47.  0 
21.  2 


At  this  time,  the  temperature  distribution  is  as  in  Fig.  2.8. 


11 


G.  H.  Higgins,  private  communication. 
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Fig.  2.  5. 
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Fig.  2.  8. 
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CONCLUSIONS 

Much  of  the  phenomena  that  occur  during  the  first  100  msec  after  a  con¬ 
tained  underground  nuclear  explosion  in  tuff  can  be  predicted  with  considerable 
accuracy.  Calculations  are  now  being  made  of  other  proposed  Plowshare  ex¬ 
periments. 


NOTE  ADDED  DECEMBER  1959 

The  elasto-hydrodynamic  model  used  for  these  calculations  is  being  gen¬ 
eralized  to  a  reversibly  cracked  elasto-plasto-hydrodynamic  system:  Radial 
and  tangential  cracks  are  opened  and  closed  reversibly.  At  initial  time  the  diag¬ 
onal  components  of  the  elastic  stress  are  compressive  by  amount  pgh  (p  =  density, 
g  =  acceleration  of  gravity,  h  =  depth  of  burial).  When  one  of  these  compenents 
(radial  or  tangential)  becomes  tensile  (due  to  the  deformation  resulting  from  the 
explosion),  it  is  set  to  zero  during  the  current  time  cycle  and  the  other  com¬ 
ponent  is  calculated  from  its  strain.  When  the  stress  becomes  compressive 
again,  the  material  is  treated  as  a  linear  elastic  solid  with  no  cracks. 

Under  high  confining  pressures  the  rock  flows  plastically  and  supports 
large  shearing  stresses.  The  magnitude  of  the  shear  is  calculated  from  a  function 
of  the  pressure  (where  pressure  is  the  average  of  the  stress  components),  and 
the  sign  is  determined  by  the  sign  of  the  radial  velocity.  The  pressure  is  com¬ 
puted  from  the  initial  shock  Hugoniot.  Upon  liquifaction  a  transition  is  made  to  the 
hydrodynamic  phase  (no  shear).  The  main  effect  of  the  plasto-dynamic  theory  is 
to  damp  the  system,  i.  e.  ,  slow  the  expansion  of  the  cavity  and  increase  the  rise 
time  of  the  disturbances  generated  by  the  explosion.  Experiments  are  being 
planned  to  determine  the  constants  needed  for  the  theory. 

The  various  states  and  transitions  are  indicated  in  the  following  diagram: 
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UNEC  Material  Behavior 


r  >  0  (tension) 


Vaporize 


-135- 


UCRL-5675 


3.  SPALLING  AND  LARGE  BLASTS 
John  S.  Rinehart 

Professor  of  Mining  Engineering,  Colorado  School  of  Mines 

INTRODUCTION 

Large  contained  nuclear  blasts  impart  tremendous  forces  to  their  imme¬ 
diate  surroundings  and  these  forces,  becoming  distributed  in  accordance  with 
well  known  laws  of  physics,  produce  damage  in  areas  both  near  and  remote 
from  the  blast-  Many  qualities  of  damage  manifest  themselves  and  of  these 
one  of  the  most  esoteric  is  spalling,  i.  e.  ,  separation  of  rock,  frequently  close 
to  a  free  surface,  as  a  consequence  of  preferential  partitioning  of  the  momen¬ 
tum  of  the  stress  wave,  The  elemental  principles  of  the  dynamics  and  me¬ 
chanics  of  spalling  in  soils  and  rocks  are  discussed  here.  These  principles 
are  then  used  to  appreciate  more  fully  the  significance  of  the  observations  of 
ground  movement  at  the  Rainier  blast- 

PROPAGATION  OF  TRANSIENT  DISTURBANCES 

When  a  buried  nuclear  device  explodes,  the  action  of  the  blast  is  not  in¬ 
stantaneously  transmitted  to  regions  remote  from  the  blast-  Changing  stress 
situations  will  be  communicated  through  most  materials  at  velocities  of 
several  thousand  feet  per  second,  the  exact  value  depending  upon  the  material, 
the  type  of  stress  distribution,  the  state  of  stress  within  the  medium,  and  the 
various  boundaries  involved-  In  large  earth  masses,  velocities  of  two  general 
types  are  encountered: 

The  dilatational  velocity,  c^,  the  velocity  of  propagation  of  longitudinal 
waves  in  a  medium  of  infinite  extent,  is  given  by 

-1/2 

» 

where  K  is  the  bulk  modulus;  p  is  the  density  of  the  material;  and  v  is  its 
Poissonfs  ratio-  Variations  of  K,  p,  and  v  with  stress  affect  the  value  of  c^. 

Shearing  displacements  advance  through  a  material  with  a  velocity  given 

by  ct 

ct  =  (G/p)l/2  , 
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where  G  is  the  rigidity  modulus. 

If  the  stress  level  is  exceedingly  high,  the  material  will  behave  as  a  fluid, 
having  no  rigidity,  and  transverse  or  shear  waves  cannot  exist  within  it.  Un¬ 
der  these  conditions,  stress  changes  will  be  propagated  with  a  velocity,  c^, 
given  by 

cf=(K/p)l/2 

The  state  of  stress  will  affect  both  K  and  p.  Indeed,  at  the  exceedingly  high 
pressures  produced  by  nuclear  blasts,  changes  in  values  of  K  and  p  can  amount 
to  30  or  40%. 

A  stress  disturbance  passes  through  a  body  in  the  form  of  a  transient 
wave  of  particle  motions.  A  point  or  infinitesimally  small  particle  within 
the  body  will  be  subjected  to  forces  which  impart  transitory  motions  to  it. 

The  instantaneous  velocity,  v,  of  such  a  particle  is  directly  related  to  the 
instantaneous  stress  cr  at  that  point  by  the  equation 

cr  =  pcv  , 

where  c  is  the  velocity  of  propagation,  viz.  ,  c^,  c^,  or  of  the  disturbance. 
Particle  velocity  in  a  longitudinal  wave  is  parallel  to  the  direction  of  propaga¬ 
tion  of  the  wave,  moving  along  with  an  advancing  compressional  wave  and  op¬ 
posite  to  an  advancing  tension  wave.  Particle  motion  within  a  transverse 
wave  is  normal  to  the  direction  of  propagation  of  the  wave.  In  a  typical  granite 
a  particle  velocity  of  3  fps  corresponds  roughly  to  a  stress  of  1000  psi. 

A  typical  disturbance,  such  as  might  have  been  generated  by  an  under¬ 
ground  nuclear  blast,  is  shown  in  Fig.  3.  1.  The  particular  parameters  used 
in  plotting  the  curve  are  selected  for  convenience.  Note  that  there  is  a ■  Ill- 
correspondence  between  stress  and  particle  velocity  and  that  whether  time  or 
distance  is  used  as  the  abscissa  will  depend  upon  whether  the  phenomenon  is 
viewed  at  all  points  in  space  simultaneously  or  at  a  single  point  through  which 
the  disturbance  is  passing.  The  instantaneous  velocity  of  propagation  of  any 
point  on  the  wave  depends  upon  the  level  of  stress  at  that  point,  being  given 
by  the  expression 

C  =  [(d<r/d€)/pl  1/2, 

where  (dtr/dc )  is  the  slope  of  the  stress-strain  curve  at  that  stress  level.  Un¬ 
less  the  stress-strain  curve  is  linear,  which  is  seldom  the  case,  the  different 
points  on  the  disturbance  will  range  in  velocities.  Thus  the  shape  of  the 
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disturbance  will  be  continually  changing  as  it  moves  along.  In  the  particular 
case  shown  in  Fig.  3.  1,  the  wave  would  tend  to  steepen  its  front  since  the 
higher  stresses  would  be  traveling  at  greater  velocities  than  the  lower 
stresses. 

Rocks  exhibit  a  wide  variety  of  stress- strain  curves.  Typical  curves  for 
limestone  and  sandstone  are  shown  in  Fig.  3.  2.  Note  that  the  limestone  curve 
is  concave  downward  which  means  high  stresses  will  travel  at  relatively  low 
velocities  while  the  sandstone  is  concave  upward,  implying  that  the  high 
stresses  will  travel  at  high  velocities.  In  limestone,  the  waves  will  tend  to 
elongate  and  flatten  out,  with  just  the  reverse  being  true  in  sandstone. 

DYNAMICS  OF  SPALLING 

Viewed  simply  a  material  is  said  to  spall  when  one  segment  of  its  parts 
and  moves  off  from  another  segment.  This  will  occur  if  the  velocity  of  the 
one  segment  is  greater  than  the  velocity  of  the  other,  as  illustrated  in  Fig. 

3.  3.  Stress  waves,  on  reaching  a  free  surface  and  reflecting,  are  prone  to 
give  rise  to  this  type  of  partition  of  velocity  or  momentum.  Consider  the  semi¬ 
infinite  body  of  Fig.  3.4,  in  which  a  stress  wave  of  the  shape  shown  is  advanc¬ 
ing  toward  a  free  surface.  Assume,  further,  that  the  stress- strain  curve  is 
linear.  When  a  disturbance  strikes  the  free  surface,  it  will  be  reflected  as 
a  tension  wave  without  change  in  form.  The  incident  compression  wave  and 
its  reflected  tension  counterpart  will  interfere  with  each  other  as  indicated. 

The  resultant  distribution  of  stress  within  the  body  at  a  slightly  later  instant 
is  shown  in  the  right  hand  drawing,  the  tension,  AB,  increasing  as  the  re¬ 
flected  wave  moves  to  the  left.  At  some  point,  the  material  will  no  longer  be 
able  to  support  this  tension;  it  will  fracture;  and  a  spall  will  fly  off,  trapping 
much  of  the  momentum  of  the  wave.  The  stress  required  to  rupture  the  ma¬ 
terial  is  called  its  critical  normal  fracture  stress,  cr  .  In  metals,  this  frac- 

c 

ture  stress  is  exceedingly  high,  corresponding  to  differential  particle  veloc¬ 
ities  of  100  to  200  ft/ sec.  In  rocks,  the  stress  is  quite  low,  being  in  the 
neighborhood  of  500  psi  (1  to  2  fps  particle  velocity  for  granite).  The  critical 
normal  fracture  stress  is  an  extremely  important  mechanical  property  of  the 
material  and  is  the  factor  that  governs  parting  of  the  material.  Note  that  two 
factors  are  important  in  describing  spalling:  the  shape  of  the  stress  wave, 
and  the  critical  normal  fracture  stress.  Determination  of  the  thickness  of 
the  spall  is  a  simple  geometrical  problem,  provided  cr^  and  the  shape  of  the 
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wave  are  known.  Thickness  of  spall  will  be  equal  to  1/2  the  distance  within 
the  wave  that  corresponds  to  a  decrease  in  stress  equal  to  ov. 

Multiple  spalling  occurs  when  the  stress  level  within  the  wave  is  more 
than  double  the  critical  normal  fracture  stress  of  the  material.  In  such  cases 
several  juxtaposed  and  parallel  spalls  are  generated  as  illustrated  in  Fig.  3.  5. 
These  arise  in  the  following  way:  Assume  that  the  first  spall  is  generated  in 
the  manner  just  described.  The  remainder  of  the  wave  will  suddenly  find  it¬ 
self  impinging  on  a  freshly  created,  free  boundary  surface.  The  result  will 
be  the  formation  of  another  spall.  This  process  will  repeat  itself  until  the 

stress  level  of  the  wave  has  been  reduced  to  a  value  less  than  <r  .  In  the  case 

c 

illustrated  the  thickness  of  the  two  spalls  are  vastly  different.  The  first  spall 
is  quite  thick  because  the  leading  portion  of  the  wave  is  relatively  flat.  The 
subsequent  abrupt  decrease  in  stress  thins  down  the  second  spall. 

As  pointed  out  earlier,  underground  blasts  can  give  rise  to  disturbances 
of  a  wide  variety  of  shapes  which  continually  change  as  they  propagate.  Four 
idealized  shapes  are  illustrated  in  Fig.  3.  6.  In  view  of  the  strong  dependence 
that  the  thickness  of  the  spall  has  upon  shape  of  disturbance,  it  is  apparent 
that  each  particular  shape  of  disturbance  generates  its  own  characteristic 
spalls.  Given  the  shape  of  the  wave  and  the  fracture  stress  of  the  material, 
it  is,  of  course,  possible  to  describe  precisely  thes  character  of  these  spalls. 
The  simplest  shape  to  consider  is  the  sawtooth  curve  in  the  upper  left  hand 
corner  of  Fig.  3.  6.  This  shape  does  not  differ  appreciably  from  that  of  the 
disturbance  produced  in  the  consolidated  tuff  900  feet  from  the  Rainier  shot 
(see  Fig.  3.7). 

INFLUENCE  OF  MATERIAL  PROPERTIES  ON  SPALLS 

The  tensile  strength,  compressive  strengths,  competencies,-  and  aniso¬ 
tropies  in  mechanical  properties  range’  widely  in  rocks  and  soils.  The  fore¬ 
going  section  described  spalling  for  the  very  particular  situation  where  the 
rock  is  perfectly  elastic  and  anisotropic.  It  was  also  tacitly  implied  that  the 
material  could  support  a  compressive  stress  greater  in  magnitude  than  the 
tensile  stress  it  could  support.  In  real  situations  such  an  ideal  circumstance 
will  not  be  obtained.  •  Compare  now  the  three  differing  situations  shown  in 
Figs.  3.  8,  3.  9,  and  3.  10. 

The  first  of  these.  Fig.  3.  8,  illustrates  the  spalling  which  will  occur 
when  a  sawtooth  wave,  of  length  \  and  maximum,  stress  cr0  ,  strikes  the  free 
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surface  of  an  elastic  material  having  a  critical  normal  fracture  stress  <r  , 
which  is  between  1/2  and  l/3  of  <r0  .  Two  spalls  are  thrown  off.  Both  are  of 
equal  thickness,  given  by 

<r=  C<rc/2<ro)X- 

but  of  different  velocities,  the  piece  to  the  right  having  the  highest  velocity. 
Specifically,  the  two  velocities,  v^  and  v^,  are 

V1  =  2vo  ■  (V^c)  ;  v2  =  2vo  -  [^g/Pc)]' 

As  the  detached  spalls  move  to  the  right  they  will  become  separated  by  an 
amount  which  is  easily  calculable  using  the  above  equations. 

Figure  3.  9  illustrates  what  will  happen  under  similar  circumstances  when 
the  material  is  incompetent,  incapable  of  supporting  any  tension  whatsoever. 

It  is  assumed,  however,  that  a  compres sional  wave  will  be  transmitted  with¬ 
out  change  in  form.  As  soon  as  the  wave  reaches  the  free  surface,  the  mate¬ 
rial  will  begin  to  flake  off.  The  velocities  of  these  first  flakes  will  be  equal 
to  2v0  where  v0  is  given  by 

vo  =  °-o/pc  • 

Each  bit  of  flaking  absorbs  a  small  portion  of  the  momentum  of  the  incident 

wave.  The  remaining  portion  of  the  incident  wave  will  continue  to  advance 

against  the  freshly  created  free  surface,  flaking  continuing  so  long  as  any 

portion  of  the  wave  remains.  The  velocity,  v  ,  of  any  small  particle  will  be 

P 

given  by 


where  <ris  the  value  of  the  stress  at  the  front  of  the  incident  wave  at  the  time 
the  particle  leaves  the  surface.  The  initial  thickness  of  the  flaked  off  mate¬ 
rial  will  be  (l/2)X..  The  flakes  themselves  will  appear  as  a  cloud  of  debris, 
the  rear  portion  of  which  is  stationary  and  the  front  moving  forward  with  a 
velocity  2v0  .  Thus  the  extent  of  the  cloud  of  debris  is  continually  increasing, 
its  forward  portion  becoming  less  dense. 

The  situation  for  incompetent  rock  or  soil  is  examined  in  a  little  more 
detail  in  Fig.  3.  10.  Here  are  compared  the  displacements  of  points  within  a 


-1  50- 


UCRL-5675 


competent*  elastic  solid  which  does  not  spall  and  an  incompetent  mass  which 
flakes  off,  Note  that  in  the  competent  solid*  the  two  points*  A  and  B,  even¬ 
tually  undergo  the  same  displacement  but  do  not  do  so  simultaneously.  The 
reason  their  total  displacements  are  the  same  is  simply  that  the  material  does 
not  fly  apart*  and  such  must  therefore  be  the  case.  The  incompetent  material 
behaves  differently.  Point  A  initially  acquires  a  velocity  which  is  maintained 
so  that  A  moves  farther  and  farther  away  from  C,  B  moves  along  at  a  some¬ 
what  lower  velocity.  If  the  free  surface  of  such  a  mass  were  horizontal*  the 
restraining  force  of  gravity  would  affect  the  material  so  that  it  would  be 
brought  to  rest.  These  two  drawings  illustrate  graphically  the  primary  rea¬ 
son  why  damage  is  so  much  greater  in  unconsolidated  soils  and  unconsolidated 
material  than  in  solid  rock. 

The  action  of  the  wave  on  a  series  of  juxtaposed  slabs  is  illustrated  in 
Fig,  3,  11,  The  joints  between  the  slabs  are  assumed  to  be  perfect  transmit¬ 
ters  of  compression  but  possessing  no  strength  in  tension.  Each  slab  is  of 
thickness  L,  On  striking  the  free  surface*  the  incident  wave  will  be  reflected, 
with  the  reflected  portion  proceeding  backward  (into  the  left)  into  the  right¬ 
most  slab.  When  the  wave  reaches  the  first  joint*  the  rightmost  slab  will  fly 
off*  trapping  in  it  the  momentum  in  the  wave  back  to  a  distance  2L,  The  ve¬ 
locity  of  this  slab*  vg^*  will  be 

v.l  =  vo  +  V1  • 

The  second  slab  will  then  be  carried  off  in  a  similar  fashion  with  a  velocity 
given  by 

vs2  =  V1  +  v2 

and  finally,  the  third  slab  with  a  velocity 

v  o  =  +  v0  o 

s3  2  3 

THE  RAINIER  SHOT 

It  was  observed  in  the  Rainier  shot  that  one  or  more  large  slabs  of  con¬ 
solidated  tuff,  or  rhyolite,  situated  about  900  feet  above  the  blast,  as  illus¬ 
trated  in  Fig,  3,  12,  rose  to  a  height  of  9  in,  and  then  fell.  From  these  data 
and  accelerometers  placed  in  a  vertical  hole  situated  directly  above  the  blast, 
it  w*  c  possible  to  make  a  reasonable  guess  of  the  shape  of  the  stress  wave 
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which,  reached  this  top  layer  of  rock0  The  most  probable  curve  is  shown  in. 

Fig.  3.  7.  The  leading  portion  of  this  curve  approximates  the  sawtooth  waves 
whose  actions  have  been  discussed  above.  In  Fig.  3.  13,  the  wave  is  redrawn 
as  a  stress  versus  distance  plot  and  shown  approaching  the  free  surface  of 
the  rhyolite.  The  specific  gravity  of  the  rhyolite  was  taken  as  2.  7.  The  max¬ 
imum  stress  in  the  wave  is  1500  psi,  corresponding  to  a  maximum  velocity  of 
7  fps,  taking  an  average  velocity  for  propagation  of  a  wave  of  6000  fps,  close 
to  the  measured  value.  Such  a  wave  would  produce  three  spalls  of  equal  thick¬ 
nesses  in  a  material  having  a  tensile  strength  of  500  psi,  a  value  reasonably 
close  to  that  which  might  be  expected  for  rock  of  this  type.  The  three  spalls 
would  move  away  with  initial  velocities,  5.  8,  3.  5,  and  1.  5  fps.  Gravity 
would  slow  these  slabs  down  so  that  the  total  excursion  of  the  top  slab  would 
be  7  to  9  in.  Taking  a  somewhat  lower  tensile  strength,  say  250  psi,  several 
more,  but  thinner,  slabs  would  be  generated.  The  thickness  of  the  consolidated 
tuff  layer  was  230  ft.  It  seems  likely  that  the  entire  slab  probably  moved  as 
a  unit  with  additional  spalling  taking  place  in  the  material  below.  All  of  the 
observations  made  are  entirely  consistent  with  what  can  be  predicted  from  an 
appreciation  of  the  elements  of  spalling.  The  two  most  important  parameters 
needed  are  the  shape  of  the  wave  and  the  fracture  strength  of  the  rock.  In 
most  situations,  the  rock  will  be  stratified  so  that  it  will  come  off  in  layers 
of  thicknesses  preselected  by  nature. 

In  considering  the  possibility  of  damage,  it  should  be  noted  that  the  effect 
of  gravity  restricts  substantially  the  movement  of  rock,  primarily  because 
the  velocities  involved  are  so  low.  Spalling  produced  on  the  walls  and  roofs 
of  tunnels  and  other  underground  openings  will  not  have  the  mollifying  effects 
of  gravity. 
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4.  A  SEISMIC  SCALING  LAW  FOR  UNDERGROUND  EXPLOSIONS 
A.  L.  Latter  and  E.  A*  Martinelli,  RAND  Corporation 
Edward  Teller,  Lawrence  Radiation  Laboratory 

Seismic  measurements  have  been  made  on  five  underground  nuclear  ex¬ 
plosions  ranging  in  energy  release  from  about  70  tons  to  19  kilotons  TNT- 
equivalento  All  these  explosions  took  place  in  the  same  general  area  of  the 
Nevada  Atomic  Test  Site,  The  table  lists  the  names,  dates,  depths  of  burial 
and  sizes  of  the  five  shots* 


Depth  of  burial  Energy  release 


Name 

Date 

(feet) 

(kt  of  TNT -equivalent) 

Blanca 

10/30/58 

835 

19  ±  1.5 

Logan 

10/15/58 

830 

5  ±  0.4 

Rainier 

9/19/57 

790 

1.  7  ±  0. 1 

Neptune 

10/14/58 

99 

0.  09  ±  0.  02 

Tamalpais 

10/8/58 

330 

0.  07  ±  0.  01 

Relative  amplitude  of  horizontal  earth  motion  as  a  function  of  the  explo- 

sion  energy  is  shown  in  Fig.  4*  1.  These  data  were  obtained  at  the  Tinemaha 
Station  of  the  Cal-Tech  Seismological  Network  using  a  Wood- Anderson  torsion 
seismograph*  The  distance  between  the  test  site  and  the  point  .of  observation 
is  approximately  180  kilometers.  The  relative  amplitudes  shown  in  the  figure 
refer  to  the  amplitude  maxima,  which  were  shear  waves,  but  within  the  ac¬ 
curacy  of  the  measurements,  the  same  relative  values  appear  to  hold  for  the 
P-wave  motion  as  well*  (Frequency  is  not  a  function  of  yield,  but  constant*  ) 
The  observed  periods  were  similar  for  all  the  shots  and  roughly  equal  to  one 
second.  Although  the  precise  values  of  the  energy  release  are  uncertain,  a 
best  power- fit  indicates  that  the  amplitude  is  roughly  proportional  to  the  first 
power  of  the  energy.  This  dependence  is  somewhat  unexpected  since,  for  an 
elastic  wave  of  a  given  period,  the  amplitude  scales  as  the  square  root  of  the 
energy. 

An  explanation  of  the  observed  dependence  is  proposed  bn  the  basis  of  the 
following  simple  model.  Close  to  the  point  of  the  explosion  there  is  a  shocked 
region  in  which  energy  is  propagated  and  dissipated  by  nonlinear  processes. 
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Fig.  4.  1.  Relative  amplitude  of  horizontal  earth  motion  as  a  function  of  the  explosion  energy. 
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This  region  is  assumed  to  be  similar  for  all  explosions  so  that  its  radius  ja. 
scales  like  the  cube  root  of  the  energy*  consistent  with  the  assumption  that 
the  nonlinear  processes  stop  at  a  definite  shock  pressure  p0  .  In  the  region 
r  >  a  the  medium  is  assumed  to  be  elastic  and  linear  so  that  a  harmonic  anal¬ 
ysis  of  the  pressure  wave  is  permissible  in  the  sense  that  the  modes  do  not 
mix.  It  is  further  assumed  that  energy  is  being  dissipated  in  this  region  in  a 
way  which  depends  on  the  frequency  but  not  on  the  amplitude.  The  effect  of 
this  dissipation  is  to  damp  the  high-frequency  components  of  the  pulse,  allow¬ 
ing  only  the  low-frequency  components  to  be  transmitted  to  large  distances. 

In  the  elastic  region  r  >a  the  seismic  propagation  is  complicated  by  the 
free  surface,  inhomogeneities,  and  dissipation,,  If  we  were  interested  in  the 
detailed  wave  shape  of  the  distant  signal,  we  would  have  to  deal  with  these 
complications o  However,  we  are  interested  only  in  the  way  in  which  the  mag¬ 
nitude  of  the  earth  displacement  depends  on  the  radius  a  and  hence  on  thecen- 
ergy  of  the  explosion.  For  this  purpose  the  free  surface  and  inhomogeneities 
can  be  neglected,  and  the  dissipation  can  be  treated  in  a  simplified  manner. 

We  assume  only  that  the  medium  is  homogeneous  within  a  fixed  radius  r  =  R 
which  is  at  least  a  few  times  greater  than  any  of  the  values  of  a  considered. 

At  this  radius  the  displacement  is  determined  as  a  function  of  time  for  each 
value  of  a.  This  function  is  then  regarded  as  the  boundary  condition  for  the 
seismic  propagation  in  the  region  r  >  R.  It  is  necessary  of  course  that  Jl  be 
much  less  than  the  distance  between  the  explosion  point  and  any  surface  of  dis¬ 
continuity  (in  particular  the  free  surface)  in  order  that  interference  between 
the  outgoing  wave  and  reflected  and  refracted  waves  at  r  =  R  be  negligible. 

Since  R  is  chosen  independent  of  a,  all  the  ^-dependence  of  the  seismic  signal 
is  contained  in  the  boundary  condition  at  R  . 

With  this  model  a  simple  argument  can  be  given  to  determine  the  depen¬ 
dence  of  the  amplitude  of  the  distant  seismic  signal  on  the  energy  of  the  ex¬ 
plosion.  It  follows  from  the  assumption  that  the  pressure  p0  at  the  elastic- 
inelastic  interface  a  is  independent  of  the  yield,  that  the  amplitude  of  the  dis¬ 
placement  produced  at  a  is  proportional  to  ^a  .  Since  ^R  is  taken  reasonably 
large  compared  to  a  ,  only  the  radiated  wave  need  be  considered  and  from  a 
to  R  the  amplitude  will  decrease  by  the  factor  a/R.  The  amplitude  of  the  dis- 
placement  at  R  is  therefore  proportional  to  a  /R. 

If  the  pressure  at  a  remains  high  for  a  long  time  compared  to  the  period 
a/c  in  which  the  sound  velocity  c  traverses  the  distance  _a ,  then  the  displacement 
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at  a  will  reach  its  maximum  value  in  a  time  approximately  equal  to  a/c  and 
will  thereafter  be  rapidly  damped  to  the  static  value.  It  follows  that  at  the 
point  R  a  disturbance  will  arrive  whose  amplitude  is  proportional  to  a  /R  and 
whose  duration  is  proportional  to  a/c.  In  the  model  we  are  now  considering 
this  disturbance  will  have  essentially  only  outwardly  directed  components.  The 
Fourier  amplitude  of  a  wave  of  sufficiently  long  wavelength  will  then  be  propor¬ 
tional  to  the  product  of  amplitude  and  duration  at  the  point  R.  Thus  we  obtain  an 
amplitude  which  is  proportional  to  a  and  therefore  proportional  to  the  energy 
of  the  explosion. 

Another  limiting  case  of  interest  is  one  in  which  the  pressure  is  not  sus¬ 
tained  but  disappears  in  a  time  proportional  to  a/c.  In  this  case  the  outward 
motion  is  followed  by  a  similar  inward  displacement  in  a  time  which  is  pro¬ 
portional  to  a/c.  The  result  at  R_  is  the  superposition  of  a  positive  and  neg- 
ative  pulse,  each  of  an  amplitude  proportional  to  a  /R  and  of  a  duration  pro¬ 
portional  to  a/c,  and  separated  by  a  time  interval  again  proportional  to  a/c. 

The  Fourier  component  of  a  long  wave  will  then  be  proportional  as  before  to 
the  product  of  a  /R  and  a/c,  but  will  furthermore  be  proportional  to  the  time 
interval  between  the  positive  and  negative  pulses.  This  introduces  a  further 
factor  a/c  and  the  resulting  Fourier  amplitude  is  therefore  proportional  to 
a  or  the  -power  of  the  energy  of  the  explosion. 

These  simple  arguments  will  now  be  verified  by  a  more  rigorous  calcula¬ 
tion.  Since  the  explosion  is  spherically  symmetrical  in  the  region  r  <  R,  the 
equation  of  elasticity  in  this  region  reduces  to  the  one-dimensional  wave  equa¬ 
tion  in  terms  of  a  displacement  potential  which  is  related  to  the  displace¬ 
ment  by 

C  =3F~(t)  • 

The  one-dimensional  wave  equation  is  satisfied  if  i| ;  =^(t),  where  r  =t  -  ^  , 

the  retarded  time.  The  form  of  t  has  been  chosen  to  represent  an  outgoing 
spherical  wave. 

We  shall  analyze  the  problem  by  relating  the  Fourier  analysis  of  the 
source  pressure  p  at  point  ^a  to  the  displacement  £,  at  the  greater  distance 
R  which  is  in  turn  related  to  the  displacements  at  great  distances.  The  pres¬ 
sure  p(t)  at  r  =  a  is  connected  to  the  displacement  and  its  derivative  at  r  =  a 
through  Hooke's  law: 

p(t)  =  -  (X  +  2\i)  -  2\ 


r=a 


(2) 
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where  X  and  p  are  the  Lame  constants.  In  terms  of  the  displacement  potential, 

Pit)  =  rw  -  ^  +(t>  .  (3) 

ac  a  c  a 

Denoting  the  Fourier  transform  of  p(t)  by  p*  and  the  Fourier  transform  of 
Z,  by  it  follows  from  Eqa(l)  and  (3)  that  ^  at  r  =  R  is 


To  treat  the  dissipation  we  assume  that  the  earth  acts  as  a  filter  charac¬ 
terized  by  a  response  function  S{oa)0  This  means  that  if  the  displacement  at 
r  =  R  is  given  by 

/OOO 

t,  (r)  =  \  C  (w)  e1WTdw  , 

oo 

then  the  distant  signal  will  contain  Fourier  components  given  by 

tdistM  =  f(*)  S(o>)  .  (5) 

r  is  time  measured  after  the  arrival  of  the  wave  front,  i.  e.  ,  retarded  time. 

The  function  S(w)  depends  on  the  nature  of  the  physical  mechanism  which 
causes  the  energy  loss.  It  will  not  be  necessary  for  us  to  consider  this  mech¬ 
anism  in  detail  nor  will  we  need  to  know  the  precise  form  of  S(w).  We  shall 
assume  only  that  S(oo)  has  the  property  of  admitting  low  frequencies  and  sup¬ 
pressing  high  frequencies.  The  cut-off  frequency  of  S(w)  will  further  be 
assumed  to  be  much  less  than  the  natural  elastic  frequency  o)0  of  the  inner 
region  (r  <  a). 

In  this  case  w  can  be  neglected  compared  to  ojq  and 
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The  constant  of  proportionality  in  this  equation  is  independent  of  a,  If  we 
assume  that  p(t)  is  a  jump  function,  equal  to  zero  for  t<0  and  equal  to  p0  e 
for  t>0,  then 


9 


Po 

a  +  ito  s 


(7) 


where  p0  is  a  constant  independent  of  a. 

There  are  several  cases  to  be  considered  depending  upon  the  value  of  a. 

If  a  is  independent  of  ji*  'p'will  be  independent  of  a  and  the  amplitude  of  the 

distant  signal  will  be  proportional  to  a  and  therefore  proportional  to  the  first 

power  of  the  explosion  energy.  If  a  depends  on  a  but  is  much  smaller  than 

cocs  i.  e.  ,  the  pressure  falls  off  very  slowly  compared  to  the  period  of  the 

transmitted  waves,  'p'is  again  essentially  independent  of  a3  and  again  the 

distant  signal  goes  like  the  energy  to  the  first  power.  If9  however,  a 

and  a  is  also  proportional  to  go  =  c/a*  then  it  follows  that  the  amplitude  of  the 
. .  .  °  4  4 

distant  signal  is  proportional  to  a  or  energy  to  the  -power.  To  fix  the  de¬ 
pendence  more  precisely  requires  an  analysis  of  the  nonlinear  region  which 
we  shall  not  undertake  here. 

The  data  shown  in  Fig .4.1  are  not  adequate  to  distinguish  between  the  first 

4 

power  and  the  -powen  for  the  amplitude  yield  relationship.  In  fact  the  data 
are  intrinsically  uncertain  due  to  lack  of  knowledge  of  the  precise  energy  re¬ 
lease  of  the  shots*  particularly  for  Blanca  and  Tamalpais.  Furthermore*  as 
can  be  seen  from  the  table*  the  depth  of  burial  for  the  five  shots  varies  over 
a  considerable  range  and  this  may  introduce  effects  not  described  by  the 
model. 

We  have  seen  how  the  displacement  depends  on  the  energy  of  the  explosion. 

It  is  also  of  interest  to  see  how  the  acceleration  depends  on  the  energy  of  the 

explosion.  This  can  be  done  readily  from  the  preceding  arguments.  The 

2 

Fourier  analysis  of  the  acceleration  is  co  times  the  Fourier  analysis  of  the 
displacement.  For  the  distant  seismic  signal  the  Fourier  analysis  of  the  dis- 
placement  cuts  off  in  the  neighborhood  of  gocj>  and  the  factor  oj  does  not  intro¬ 
duce  any  further  a-depend ence.  At  large  distances  therefore  the  acceleration 
and  the  displacement  depend  in  exactly  the  same  way  on  the  energy  of  the  ex¬ 
plosion.  Close-in*  however*  where  the  high  frequency  waves  have  not  been 

dissipated*  the  Fourier  analysis  of  the  displacement  cuts  off  in  the  neighbor- 
«  2 

nood  of  io  instead  of  go  ,  and  the  factor  gj  does  introduce  some  additional 
o  c 
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a-dependence.  In  this  case,  since  w  =c/a,  the  Fourier  analysis  of  the  accel- 

—  2  ° 

eration  is  proportional  to  (l/a  )  times  the  Fourier  analysis  of  the  displace- 
ment.  From  Eq„  (4)  it  can  be  seen  that  the  Fourier  analysis  of  the  displace¬ 
ment  at  close-in  distances  (but  still  in  the  radiation  zone)  is  proportional  to 
2 

a  .  It  follows  that  the  acceleration  at  close-in  distances  is  independent  of  ^a, 
and  therefore  independent  of  the  energy  of  the  explosion.  Experiments  should 
be  made  to  check  this  amusing  consequence  of  the  model. 

We  are  grateful  to  Dr0  Frank  Press  for  providing  us  with  the  seismic 
data  contained  in  Fig.  4.  1. 


/ro 


- 1 63- 


UCRL-5675 


LIST  OF  PREVIOUS  PLOWSHARE  AND/OR  RELATED  REPORTS 


Report  No. 


Title 


UCRL-4659 
UCRL-5026 
UCRL-5124  Rev.  I 
UCRL-5253 
UCRL-5257  Rev. 
UCRL-5281 

UCRL-5457 

UCRL-5458 

UCRL-5538 

UCRL-5542  Rev. 

UCRL-5623 

UCRL-5626 

UCRL-5709 

UCRL-5757 

UCRL-5766 

UCRL-5882 

UCRL-5932 

UCRL-5949 

UCRL-601 3 


Deep  Underground  Test  Shots. 

Non-Military  Uses  of  Nuclear  Explosions 

Phenomenology  of  Contained  Nuclear  Explosions. 

Industrial  Uses  of  Nuclear  Explosives. 

Peaceful  Uses  of  Fusion. 

Temperatures  and  Pressures  Associated  with  the 
Cavity  Produced  by  the  Rainier  Event. 

Large  Scale  Excavation  with  Nuclear  Explosives. 

Mineral  Resource  Development  by  the  Use  of  Nuclear 
Explosives. 

Evaluation  of  the  Ground  Water  Contamination  Hazard 
from  Underground  Nuclear  Explosions. 

Properties  of  the  Environment  of  Underground  Nuclear 
Detonations  at  Nevada  Test  Site.  Rainier  Event. 

Radioactivity  Associated  with  Underground  Nuclear 
Explosions. 

Underground  Nuclear  Detonations. 

Hydroclimatology  and  Surface  Hydrology  of  San 
Clemente  Island. 

Geologic  Studies  of  Underground  Nuclear  Explosions 
Rainier  and  Neptune. 

The  Neptune  Event.  A  Nuclear  Explosive  Cratering 
Experiment. 

Chemical  Reactions  Induced  by  Underground  Nuclear 
Explosions. 

The  Soviet  Program  for  Industrial  Applications  of 
Explosions. 

An  Application  of  Nuclear  Explosives  to  Block  Caving 
Mining. 

Probing  the  Earth  with  Nuclear  Explosions. 


164- 


UCRL-5675 


LIST  OF  PREVIOUS  PLOWSHARE  AND/OR  RELATED  REPORTS 


Report  No. 


UCRL-5675 

UCRL-5676 
UCRL-  567*7 


UCRL-5678 


UCRL-5679 


Title 

Plowshare  Series:  Report  No.  2.  Proceedings  of  the 
Second  Plowshare  Symposium  (Held  at  San  Franciscos 
May  13-15,  1959): 

Parti  :  Phenomenology  of  Underground  Nuclear 
Explosions. 

Part  II  :  Excavation. 

Part  III:  Recovery  of  Power  and  Isotopes  from 

Contained  Underground  Nuclear  Explosions. 

Part  IV:  Industrial  Uses  of  Nuclear  Explosives  in 

the  Fields  of  Water  Resources,  Mining,  Chemical 
Production,  and  Petroleum  Recovery. 

Part  V  :  Scientific  Applications  of  Nuclear  Explosives 
in  the  Fields  of  Nuclear  Physics,  Seismology, 
Meteorology  and  Space. 


LEGAL  NOTICE 


This  report  was  prepared  as  an  account  of  Government  sponsored  work. 
Neither  the  United  States,  nor  the  Commission,  nor  any  person  acting  on 
behalf  of  the  Commission: 

Ac  Makes  any  warranty  or  representation,  expressed  or  implied,  with 
respect  to  the  accuracy,  completeness,  or  usefulness  of  the  information  con¬ 
tained  in  this  report,  or  that  the  use  of  any  information,  apparatus,  method, 
or  process  disclosed  in  this  report  may  not  infringe  privately  owned  rights;  or 

B*  Assumes  any  liabilities  with  respect,  to  the  use  of,  or  for  damages 
resulting  from  the  use  of  any  information,  apparatus,  method  or  process  dis¬ 
closed  in  this  report* 

As  used  in  the  above,  person  acting  on  behalf  of  the  Commission  fr 
includes  any  employee  or  contractor  of  the  commission,  or  employee  of  such 
contractor,  to  the  extent  that  such  employee  or  contractor  of  the  Commission, 
or  employee  of  such  contractor  prepares,  disseminates,  or  provides  access 
to,  any  information  pursuant  to  his  employment  or  contract  with  the  Commis¬ 
sion,  or  his  employment  with  such  contractor* 


Printed  in  USA.  Price  $3,  00  Available  from  the 
Office  of  Technical  Services 
U,  S.  Department  of  Commerce 
Washington  25,  D,  C. 


